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This article will give an overview on both theoretical and experimental developments concerning states with 
lattice symmetry breaking in the cuprate high-temperature superconductors. Recent experiments have provided 
evidence for states with broken rotation as well as translation symmetry, and will be discussed in terms of nematic 
and stripe physics. Of particular importance here are results obtained using the techniques of neutron and x-ray 
scattering and scanning tunneling spectroscopy. Ideas on the origin of lattice-symmetry-broken states will be 
reviewed, and effective models accounting for various experimentally observed phenomena will be summarized. 
These include both weak-coupling and strong-coupling approaches, with a discussion on their distinctions and 
connections. The collected experimental data indicate that the tendency toward uni-directional stripe-like ordering 
is common to underdoped cuprates, but becomes weaker with increasing number of adjacent Cu02 layers. 
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I. INTRODUCTION 

High-temperature superconductivity in the copper oxides 
constitutes one of the most fascinating and challenging prob- 
lems in modern condensed matter physics. Since its dis- 
covery in 1986 by Bednorz and Miiller ( 1) it has influenced 
and inspired a vast variety of both experimental and theo- 
retical developments, ranging from the tremendous improve- 
ments in techniques like photoemission and scanning tunnel- 
ing microscopy over the development of theoretical tools for 
strongly correlated and low-dimensional systems to the dis- 
covery of fundamentally new states of matter. 



A. Doped Mott insulators 

Superconductivity in the copper oxides arises from dop- 
ing of half-filled Mott insulators. 1 While a number of low- 
temperature properties of the superconducting state are com- 
patible with BCS-like pairing (2) of d-wave symmetry, the 
normal state is highly unconventional and appears to violate 
Fermi-liquid properties (except, perhaps, at strong overdop- 
ing), in contrast to most other superconductors where pairing 
derives from a Fermi-liquid metallic state. 

Understanding doped Mott insulators is at the heart of the 
high-T c problem (H;|3), however, here both our conceptual and 
methodological toolboxes are still rather limited. Phenomeno- 
logically, the pseudogap regime in underdoped compounds is 
not understood, as is the doping evolution of the normal-state 
Fermi surface and the fate of local magnetic moments. Further 
puzzles are connected to the linear-in-T resistivity around op- 
timal doping above T c and the asymmetry between electron- 
and hole-doped materials (H|;|^;|3). 

It has become clear that doped Mott insulators are charac- 
terized by a plethora of competing ordering tendencies, with 
commensurate antiferromagnetism and superconductivity be- 
ing the most prominent ones. Others are incommensurate spin 
and charge density waves, orbital magnetism with circulating 
currents, and exotic fractionalized states with topological or- 
der (Sii). 

Incommensurate uni-directional s pin and charge density 
waves, often dubbed "stripes" (H M [H El H 0), 
play an interesting role: Those states were first predicted in 
mean-field studies of Hubbard models d M m Q 

and 

later observed in La2- a; Ba 2 ;Cu04 and La2-a-j,Ndj,Sr 2; Cu04 
( fl9l I20I) . belonging to the so-called 214 family of cuprate 
compounds. Those stripe states break the discrete translation 
and rotation symmetries of the square lattice underlying the 
Q1O2 planes. A conceptually related state is an (Ising) "ne- 
matic" which only breaks the lattice rotation symmetry and 
which may occur as intermediate state upon melting of stripe 
order dUlH]). 



More precisely, in terms of a three-band model of Cu 3d x 2_ y 2 and O 
2p x , y orbitals, the C11O2 planes at half-filling are charge-transfer, rather 
than Mott, insulators. 



Stripes and nematics, together with the associated phase 
transitions and fluctuation phenomena, constitute the focus of 
this review article. While stripe states were originally thought 
to be a very special feature of the 214 compounds, this view 
has changed. The last few years have seen exciting experi- 
mental progress: Signatures of states with lattice symmetry 
breaking have been identified in a number of other cuprates 
and in a number of different probes, and the dependence of 
the ordering on external parameters has been mapped out in 
more detail. In parallel, stripe-like charge ordering has been 
identified and investigated in other correlated oxides, most 
importantly nickelates l22t I23I) and manganites {24]). While 
there are some important differences between those systems 
and cuprates, it is likely that the phenomena are related. 

Taken together, the observations suggest that stripe and ne- 
matic states are an integral part of the physics of doped Mott 
insulators, qualifying them an interesting subject of funda- 
mental research on their own right. However, their role for su- 
perconductivity in cuprates is unclear at present. On the one 
hand, interesting theoretical ideas of stripe-driven or stripe- 
enhanced superconductivity have been put forward. On the 
other hand, in the cuprates there seems to be an anticorrela- 
tion between T c and the strength of stripe ordering phenom- 
ena, and no (static) stripe signatures have been identified in 
the cuprates with the highest T c . I will return to this discus- 
sion towards the end of the article. 



B. Why another article on stripes? 

Both experimental and theoretical aspects of strip es have 
been subject of other reviews in the past (H [13 EE [H CI 
[Til) . In addition, stripes have been discussed in more gen- 
eral introductory or review articles on cuprate superconduc- 
tors, both from experimental (H; [2§ U% and theoretical 
•HI S S S H; HH HI) perspectives. The most recent and 
comprehensive review closely related to the present subject is 
the one by Kivelson et al. dllb on "How to detect fluctuating 
stripes in the high-temperature superconductors", published in 
2003. 

The main motivation for a new article is that a number of 
key experiments were performed in the last few years (i.e. af- 
ter 2003), such that those and subsequent theoretical devel- 
opments are not covered in the above articles. I feel that the 
recent developments have changed the perception of the com- 
munity with regard to states with broken lattice symmetries, 
making this old subject a very timely one. 



C. Focus 

The goal of this article is to review concepts, experimen- 
tal results, and theoretical ideas relevant to lattice-symmetry- 
broken states in doped cuprate superconductors. The em- 
phasis is on uni-directional spin and charge density waves 
(stripes) and on nematic states with broken rotational sym- 
metry. These types of order can co-exist with bulk supercon- 
ductivity, and modulated superconducting states will naturally 
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FIG. 1 Left: Schematic phase diagram of hole-doped cuprates in the temperature-doping plane. Shown are the phases with Neel-like 
antiferromagnetic order (AF) and superconductivity (SC) as well as the pseudogap, strange-metal, and overdoped Fermi-liquid (FL) regimes. 
The phases with lattice symmetry breaking, discussed in this article, occur primarily in the underdoped regime at low temperatures (shaded). 
Various features of the phase diagram are not settled, such as shape of the crossover lines near the superconducting dome, the character of 
the pseudogap line, and the character of the overdoped FL-like regime. Not shown are additional crossovers associated with the onset pairing 
fluctuations and with spin glass behavior. Right: Electronic structure of the CuC>2 planes, with a unit cell containing a Cu 3d and two O 2p 
orbitals. 



appear in the discussion. 

The insulating regime of zero and very small doping will 
be discussed only briefly, as the ordered states occurring there 
can likely be described by quasi-classical concepts and are of 
less relevance to cuprate superconductivity (with caveats to be 
noted later). 

An exciting upcoming topic is that of circulating-current 
states, of the rf-density-wave type d30f> and of the type pro- 
posed by Varma J3ll [32h . Here, experiments and interpreta- 
tions are partially controversial, such that I believe, at present, 
it is too early for an extensive review. Therefore, the cover- 
age of circulating-current states will be brief and restricted to 
a few important experimental and theoretical results. 

Even with these restrictions, the amount of published liter- 
ature is vast. Therefore I will focus most of the discussion to 
works which appeared over the last few years, and the cover- 
age of older results will be restricted to a few key references. 



D. Outline 

To set the stage, Sec. [TT] will introduce the phenomenol- 
ogy of ordered phases in the language of symmetries, order 
parameters, and Landau theory. An important aspect is the 
coupling of order parameters to structural distortions and to 
quenched disorder arising from dopant atoms. Sec. HLT1 will 
give on overview on experimental data, which either estab- 
lish the existence of symmetry-breaking order or else can be 
used to extract characteristics of such order. Important classes 
of results, providing direct evidence for order, are those from 
neutron and x-ray scattering and from scanning tunneling mi- 
croscopy and spectroscopy. 

The remainder of the article is devoted to theory: In Sec. HVl 
microscopic considerations in the language of Hubbard and t- 



J models (or variants thereof) will be reviewed, with the fo- 
cus on identifying mechanisms and conditions for symmetry 
breaking. Sec. [V] then describes concrete theoretical efforts 
in modeling experimental characteristics of such symmetry- 
broken phases - due to the difficulties with microscopic or 
ab-initio approaches, most published work in this area is based 
on effective phenomenological models which take some input 
from experiment. Finally, Sec.[Vl]contains attempts to discuss 
the "global picture", i.e., the implications of both experiments 
and effective theories for our understanding of cuprate high- 
temperature superconductivity. Conclusions and an outlook 
will wrap up the article. 

To shorten the notations, the following abbrevia- 
tions for cuprate compounds will be frequently used: 
LSCO for La 2 _ zSr^CuC^, LBCO for La^Ba^aiC^, 
LESCO for La 2 - 2; -yEuySr 2; Cu04 with y = 0.2, 
LNSCO for Laa-^-yNd^Sr^CuC^ with y = 0.4, 
YBCO for YBa 2 Cu 3 6+(5 , BSCCO or Bi-2212 for 
Bi 2 Sr2CaCu20 8+<5 , Bi-2201 for Bi 2 Sr 2 Cu0 6+< 5, and 
CCOC for Ca 2 _ 2; Na 2; Cu0 2 Cl 2 . The doping level will be 
specified as LSCO-0.12 or YBCO-6.6. Moreover YBCO-124 
denotes YBa 2 Cu 4 Os. 



II. ORDERED PHASES: PHENOMENOLOGY 

This section will introduce the order parameters (in the Lan- 
dau sense) of various ordering phenomena which have been 
observed or discussed in the cuprates. The orders break sym- 
metries of the underlying microscopic models, such as the 
translation symmetry, the C^h point group symmetry of the 
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two-dimensional (2d) C11O2 square lattice, 2 and charge con- 
servation. The coupling between order parameters will be 
discussed in terms of Landau theory, which will be used to 
derive global scenarios for phase diagrams and sequences of 
phase transitions. 

Importantly, the following discussion will only make refer- 
ence to symmetries, but not to mechanisms of ordering phe- 
nomena. Consequently, terms like "spin density wave" or 
"charge density wave" will be used (throughout the article) 
for states with spatially modulated spin or charge densities, 
independent of whether the underlying physics is of weak- 
coupling or strong-coupling nature. 3 

An order-parameter field, 4>(r, t), is often introduced such 
that, in an ordered phase or close to the ordering transition, 
it varies slowly in space and time, i.e., the typical length and 
time scales of fluctuations are large compared to microscopic 
scales. To this end, modulations of the relevant observable on 
microscopic scales are absorbed in the definition of <j>, such 
that cf> is associated with a finite lattice momentum. 4 

Close to an ordering transition, the effective theory for 
4>(r, t) can often be written in the form of a <f> A or Landau- 
Ginzburg-Wilson (LGW) model, 5 
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here for a real n-component field <fi a (r, r), and h(r, r) is the 
field conjugate to the order parameter, tq is the (bare) mass of 
the order parameter field <j), used to tune the system through 
the order-disorder transition, c a velocity, and uq the quartic 
self-interaction. Sd yn contains the <j> dynamics, which is often 
of the form (<9 T 0) 2 ; exceptions arise from fermionic Landau 
damping or field-induced magnetic precession dynamics. 



A. Charge and spin density waves 

A charge density wave (CDW) is a state where the charge 
density oscillates around its average value as 6 



(p(R, r)> = p avg + Re e^ n <f> e (A, r) 



(2) 



2 For most of the article, only ordering phenomena within the C11O2 planes 
will be discussed. Ordering perpendicular to the planes is often short- 
ranged, with exceptions to be noted below. 

3 An introduction to weak-coupling and strong-coupling approaches of 
stripe formation is given in Sec. |IV| A discussion of the question which 
of the two is better suited to describe actual cuprate experiments appears in 
Sec.lYHl 

4 Occasionally, we will also refer to an order parameter description on 
the lattice scale, where the order parameter simply represents a Hubbard- 
Stratonovich field. 

5 Throughout this article, Einstein's summation convention will be em- 
ployed, such that indices occurring twice are summed over, e.g., (<?V) 2 = 

4>a4>a = J2 a 4>a4>a- 

6 Higher harmonics may always be present in the modulation pattern, but 
depend on microscopic details. 



The scalar order parameter cj) c is associated with the finite 
momentum (or charge-ordering wavevector) Q c . 7 In general, 
(f> c is complex, with the phase describing a sliding degree of 
freedom of the density wave; for the sites of a square lat- 
tice, exceptions are simple modulations with Q cx , Q cy = ir, 
where <p c is real. For wavevectors Q c which are incommen- 
surate with the lattice periodicity, all values of the complex 
phase of <\> c are equivalent, whereas the phase will prefer dis- 
crete values in the commensurate case - the latter fact reflects 
lattice pinning of the density wave. In particular, depend- 
ing on the phase, the resulting modulation pattern can pre- 
serve a reflection symmetry at bonds or sites of the lattices, 
then dubbed bond-centered or site-centered density wave. For 
stripe order in cuprates we will be mainly concerned with 
Q cx = {2n/N, 0) and Q cy = (0, 2tt/N) where N is the real- 
space periodicity, and N = 4 appears to be particularly stable 
experimentally. 

The modulation describes a uni-directional density 
wave, with a single wavevector Q c , which also breaks the 
point group symmetry. 8 Alternatively, multiple modulations 
with different wavevectors may occur simultaneously. If a 
wavevector and its symmetry-equivalent partners occur with 
equal amplitude, then lattice point group symmetry is unbro- 
ken. This is the case in a checkerboard state with equal mod- 
ulations along Q cx and Q cy . 

In this article, we shall use the term "charge density wave" 
in its most general sense, referring to a periodic spatial mod- 
ulation in observables which are invariant under spin rota- 
tions and time reversal; apart from the charge density this in- 
cludes the local density of states and the magnetic exchange 
energy, the pairing amplitude, or the electron kinetic en- 
ergy, all defined on the links of the square lattice. Note that 
even in a strongly ordered state, the modulation in the total 
charge density may be unobservably small because longer- 
range Coulomb interactions tend to suppress charge imbal- 
ance. This applies e.g. to valence-bond-solid (VBS, or spin- 
Peierls) order in a Mott insulator. In the described sense, a 
columnar VBS state (3) is a bond-centered CDW with Q c = 
(tt,0). 

A spin density wave (SDW) is specified by a vector order 
parameter <j) sa (r, r), a = x,y, z, and the spin density modu- 
lation is given by 



(S a {R,r))=Re 



(3) 



Here Q s = describes a ferromagnet, while the parent anti- 
ferromagnet of the cuprates has Q s — (tt, tt). The field 4> sa 
transforms as angular momentum under 0(3) spin rotations 
and is odd under time reversal. 



7 The lattice spacing ao (of the Cu square lattice) is set to unity unless 
otherwise noted. 

8 In low-symmetry crystals, a single-Q modulation may not break the point 
group symmetry. This applies to stripes in planes with orthorhombic sym- 
metry, relevant for the LNSCO, LESCO, LBCO, and YBCO, see Sec. HLGl 
below. 
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Both collinear and non-collinear SDWs can be described by 
Eq. ©: 



collinear: 
spiral: 



= e ie n a with n a real, (4) 
= nia + iri2a with n\ 2a real, ni a ri2Q = 0. 



In cuprates, both spiral and collinear order have been dis- 
cussed. In the context of stripe order, our focus will be on 
collinear order with Q sx — 27r(0.5 ± 1/M, 0.5) and Q sy = 
2tt(0.5,0.5 ± 1/M), with M = 2N such that 2Q S = Q c 
modulo reciprocal lattice vectors. 

Spin anisotropies play an important role, as they can 
freeze out low-energy fluctuations in some directions. In the 
cuprates, a combination of spin-orbit and crystal-field effects 
leads to antiferromagnetic moments lying in the a-b plane, and 
the Dzyaloshinskii-Moriya (DM) interaction results in a small 
spin canting out of the plane in La 2 Cu04 d33h . 

Charge and spin ordering can be expressed via expectation 
values of particle-hole bilinears. Here, charge and spin order 
simply correspond to ordering in the spin singlet and triplet 
channels of the particle-hole pair. Assuming a one-band de- 
scription and neglecting the time dependence, we have for 
charge order: 



(ct(f)c(?)) = F mg (r-r>) 



F c (f-f)Re 



iQc-R 



MR) 



(5) 



where cj.(f) is an electron creation operator at coordinate r 
and spin er, and R = (r + r*)/2. F avg and F c are short-ranged 
functions, characterizing the unmodulated state and the mod- 
ulation, respectively. After Fourier transformation, k in F c (k) 
refers to the internal momentum of the particle-hole pairs (in 
contrast to their center-of-mass momentum Q c ), it determines 
the structure of the modulation within a unit cell. The sim- 
plest possible order corresponds to F c (r — r') ~ Spf/t, i.e., 
F c (k) is /c-independent, and we will loosely refer to this as "s- 
wave" order. However, proposals of order with non-trivial k 
dependence have also been made. This applies to the so-called 
"(i-density wave order" d30» and d-wave checkerboard order 
(HH [35h . where the modulations are on lattice bonds, while 
the modulation on sites vanishes identically. (Note that the 
rf-density wave state breaks time reversal whereas the rf-wave 
checkerboard does not.) d-wave order parameters have the in- 
teresting feature that they couple primarily to antinodal quasi- 
particles, leaving nodal quasiparticles along the Brillouin zone 
diagonal unaffected. In general, the symmetry classification 
of F c has to be done according to the lattice point group of the 
ordered state: For a square lattice and uni-directional CDW 
order with Q along the (1,0) or (0,1) direction, s-wave and 
d 3 .2_ y 2-wave representations always mix, because the order- 
ing wave vector Q c breaks the symmetry from C4 to C 2 . Nev- 
ertheless, stripes may be dominated by either the s-wave or 
the d-wave component of F c , the latter leading to valence- 
bond stripes (Bo) , with little modulation on the square-lattice 
sites and little coupling to nodal quasiparticles. 

Fig. [2] schematically shows the real-space structure of 
stripes with period 4 in the charge sector, illustrating bond- 
centered and site-centered stripes. The STM results of 



Kohsaka et al. are most consistent with bond-dominated 
stripe order as in Fig. [2jl. 

As noted above, the terms "charge density wave" and "spin 
density wave" will be used without reference to the underly- 
ing cause of the modulation, hence "stripe" is equivalent to 
"uni-directional charge density wave". CDW and SDW states 
can be insulating (like Wigner crystals), or metallic, or super- 
conducting. 



B. Superconducting pairing 

For singlet superconducting pairing, the order parameter is 
a charge-2 scalar field tp(R, r), where R and r are the center- 
of-mass coordinates of the Cooper pair. Most superconduct- 
ing states have a zero-momentum condensate, i.e., an ordering 
wavevector Q p = 0. However, modulated pairing is possible 
as well, and we may write in analogy to Eq. {5): 



(4(r>!(r<)) = F p (r-7) 



JQ P R 



-iQ p -R 



(6) 



Because the superconducting order parameter is complex, two 
fields ip and rp are required to implement the sliding degree of 
freedom of the "pair density wave". The form factor F p de- 
termines the internal structure of a Cooper pair: For homoge- 
neous pairing, Q p = 0, the angular dependence of F p (r) leads 
to the usual classification into s-wave, d x 2_ y 2-wave etc. In 
the cuprates, there is good evidence for the pairing symmetry 



being d x 2 _ y i (|3 

Finite-momentum pairing has first been proposed by Fulde 
and Ferrell (39) and Larkin and Ovchinnikov (1401) . as a state 
realized for large Zeeman splitting of the Fermi surfaces in an 
external field. More than 40 years later, evidence for such a 
FFLO state has indeed been found in the organic supercon- 
ductor k-(BEDT-TTF) 2 Cu(NCS) 2 ©. In cuprates, finite- 
momentum pairing in zero field has been proposed in the con- 
text of stripe phases (USUI), see Sec. |V.F.2| 



C. Nematics 

States which spontaneously break the (discrete or contin- 
uous) real-space rotation symmetry of the underlying Hamil- 
tonian are often called "nematic" states. The discussion here 
will be restricted to homogeneous (Q = 0), time-reversal in- 
variant, and spin-symmetric nematic order. 9 Moreover, the 
focus will be on electronically driven symmetry breaking: 
While lattice distortions invariably follow the broken sym- 
metry, transitions involving simple structural distortions only 
shall not be considered "nematic". 



9 Modulated nematic order is equivalent to density wave order. Spin- 
antisymmetric nematic order and "spin nematic" order, corresponding to 
spontaneously broken spin rotation symmetry, will not be of interest in this 
article. 
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FIG. 2 Schematic real-space structure of stripe states with period 4 in the charge sector, a) Bond-centered and b) site-centered stripes with 
period-8 magnetic order, showing the on-site hole densities and magnetizations, c), d) Bond-centered charge-only stripes, now showing the 
modulations in both the on-site hole density and the bond strength (i.e. kinetic energy). Depending on the form factor F c in Eq. {5]l, the 
modulation in the spin-singlet sector can be primarily of c) s-wave or d) d-wave type - the latter structure is dominated by bond modulations. 



The term "nematic" originates in the physics of liquid crys- 
tals where it refers to a phase with orientational, but without 
translational, order of molecules. An electronic nematic phase 
can arise as a Pomeranchuk instability of a Fermi liquid J44l) . 
in which case the ordered phase possesses a full Fermi sur- 
face, but can also be driven by strong correlations and even 
occur in an insulator, i.e., in the absence of a Fermi surface. In 
particular, a nematic phase may occur upon melting of stripe 
order, namely when the rotation and the translation symme- 
try are broken at separate transitions. Then, a nematic is the 
intermediate phase between the stripe and disordered phases 

dUl). 10 

In the case of a continuous rotation symmetry, the order 
can be characterized by an angular momentum I > 0; e.g. for 
/ = 2, the order parameter is a director, specifying an axis 
in real space. In contrast, on the square lattice with C4 rota- 
tion symmetry, I = 2 (or rf-wave) nematic order corresponds 
to breaking the rotation symmetry down to C2. One can dis- 
tinguish d xy and d x 2_ y 2 nematic order, here we shall focus 
on the latter. The order parameter is of Ising type; globally, 
it may be defined via the electronic momentum distribution 
according to 

K = XX( C L C fe<r)' d k = COsk * ~ COsk V ( ? ) 
fe 

A local order parameter can be defined from any spin-singlet 
bond observable which is even under time reversal, like the 



In the physics of liquid crystals, phases are termed "liquid", "nematic", 
"smectic", and "crystalline", according to their broken spatial symmetries. 
Here, a smectic breaks translation symmetry in one direction, whereas a 
crystal breaks it in all directions. Some workers have applied this termi- 
nology to cuprates, identifying a conducting stripe with a "smectic" and an 
insulating Wigner crystal with a "crystal". To avoid terminology problems 
with e.g. multiple-Q modulated states, this article will employ the terms 
"CDW" and "SDW" instead, as described above. 



bond kinetic energy, 

<t>n{r) = (4(r)c CT (r + x) - ci(r)c a (r + y)), (8) 

or the bond magnetic energy or pairing strength. 

The so-defined nematic order parameter (f> n measures ro- 
tation symmetry breaking. As such, it is finite in a uni- 
directional density wave state as well: The density wave will 
induce a <j> n pro porti onal to \<p cx \ 2 - \4> cy \ 2 or \(j) sax \ 2 - 
|<^sai/| 2 ' see Sec. III. E| Strictly speaking, the density wave is 
not a nematic state, as translation symmetry is broken as well. 
However, here we shall use the term "nematic order" synony- 
mously to lattice rotation symmetry being spontaneously bro- 
ken. 11 

In recent years, nematic phases have been proposed as ex- 
planations of some puzzling experimental observations, such 
as the enigmatic ordering transition in URu2Si2 d45l) and the 
low-temperature phase near the metamagnetic transition of 
Sr 3 Ru 2 7 (S3). 



D. Loop-current order 

While SDW, CDW, and nematic order have been experi- 
mentally identified at least in some cuprates with reasonable 
confidence, other symmetry-breaking order parameters (apart 
from superconductivity) have been suggested but are more 
controversial. Most prominent are various forms of loop- 
current order, proposed as origin for the pseudogap behav- 
ior of underdoped cuprates. Charge current loops break the 
time reversal symmetry and induce orbital moments. In par- 
ticular, orbital antiferromagnets are characterized by current 



Experimentally establishing the existence of a nematic instead of a smectic 
(in the liquid-crystal terminology) would require disproving the existence 
of translation symmetry breaking, which is rather difficult. 
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loops with directions alternating from plaquette to plaquette 
on the lattice, such that the total moment is zero and hence 
there is no a global edge current. 

Varma J3H : I32I) has proposed different patterns of current 
loops within a unit cell of the CuC>2 plane, such that the or- 
der exists at wavevector Q = 0, but has vanishing total mo- 
ment. The description of these states in general requires a 
three-band model of the CUO2 plane, with Cu and O orbitals, 
however, the symmetry of the state On of Ref. (132h ca n also 
be represented in a one-band model (case G in Ref. (1471) '). The 
magnetic moments are expected to point in a direction perpen- 
dicular to the CuC<2 layers, but in principle loop-current order 
could also involve oxygen orbitals outside layers, resulting in 
canted moment directions. 

A distinct type of state was proposed by Chakravarty et al. 
d30b as resolution of the pseudogap puzzle, namely an alter- 
nating current pattern with Q — (tt, it) with <i-wave form fac- 
tor fia; |49|) . dubbed d-density wave. In fact, symmetrywise 
this state is a generalization to finite doping of the so-called 
staggered flux phase which appeared in early mean-field stud- 
ies of Hubbard-Heisenberg models (50). 

Recent neutron scattering experiments reported magnetic 
order in the pseudogap regime at Q = 0, which has been in- 
terpreted as evidence for the loop current of On type, and we 
shall come back to these experiments in Sec. IIII.A.3I Some 
other neutron scattering experiments have also been argued 
to be consistent with d-density wave order, and will be men- 
tioned in Sec. IIH. Al We shall, however, refrain from a detailed 
discussion of loop-current states in this article, considering 
that some aspects of both theory and experiment are either 
controversial or not fully settled. 



condensate ip, if>, Eq. ©, according to 

X S (4>* C 4>1 + cc.) + Xi((p* c tpip* + c.c.) 
+ A 5 (CVoV>* + <f>cM* +c.c). (10) 

The couplings A3,4 are allowed only if the ordering wavevec- 
tors obey Q c = 2Q S and Q c — 2Q p , respectively, whereas A5 
is only allowed if Q c = Q p . 

The listed terms imply, e.g., that a uni-directional density 
wave induces nematic order via A1.2 and that a collinear SDW 
with Q s induces a CDW with 2Q S via A3. As a consequence, 
a transition into a stripe-ordered state may occur as a direct 
transition from a disordered into a CDW+SDW state, or via 
intermediate nematic and CDW phases, see Fig. [3] The phase 
diagram of the corresponding Landau theory has been worked 
out in detail by Zachar et al. (1511) . Whether the intermedi- 
ate phases are realized depends again in microscopic details; 
weak-coupling theories typically give a direct transition into a 
CDW+SDW state. 

Similar considerations can be applied to the coupling be- 
tween stripes and superconducting condensates arising from 
A4.5. In particular, the co-existence of stripes and supercon- 
ductivity induces a modulated component of the condensate. 
A detailed analysis can be found in Ref. ( |53I) . 

An external Zeeman magnetic field couples only quadrat- 
ically to any of the order parameters, as all have vanishing 
uniform magnetization. The orbital part of an external field 
will have a strong influence primarily on the superconducting 
order, by inducing vortices. These considerations will be im- 
portant in discussing the field tuning of ordering phenomena 
in Sec. lIII.E.3l 



E. Order parameter coupling and global phase diagrams 

In the presence of multiple ordering phenomena, the vari- 
ous order parameters are locally coupled. The general form 
of this coupling can be deduced from symmetry arguments. 
For any two order parameters <fii and (j>2, a density-density 
coupling term v\4>i(f, t)\ 2 \4>2(^, t )\ 2 m the LGW theory is 
generically present. The sign of v depends on microscopic pa- 
rameters and decides about repulsion or attraction between (f>i 
and 4>2. For instance, if (j>i $ represent horizontal and vertical 
CDW order parameters, then v > will lead to uni-directional 
(stripe) order where v < results in bi-directional (checker- 
board) order. 

More interesting are couplings involving one order parame- 
ter linearly. Those terms are strongly constrained by symme- 
try and momentum conservation. A nematic order parameter 
(\> n in a tetragonal environment couples to CDW and SDW 
according to 

Xl4>n(\4>cx\ - \4>cy\ ) + <Wn(|<^s ax I |0say r), (9) 

note that \(t) cx \ 2 etc. carry vanishing lattice momentum. A 
CDW order parameter couples to a spin density wave <f> s , 
a uniform superconducting condensate i/jq, and a modulated 



F. Concept of "fluctuating order" 

"Fluctuating order" describes a situation on the disordered 
side of, but close to, an ordering transition, such that precur- 
sor effects of the ordering phenomenon are visible in physical 
observables. Regarding symmetries or symmetry breaking, 
"fluctuating order" is equivalent to "no order". For metallic 
systems and in the absence of additional symmetry breaking, 
a regime of fluctuating order can then be adiabatically con- 
nected to the weakly interacting Fermi liquid. 

In weakly interacting and/or high-dimensional systems, the 
fluctuation regime is usually tiny and restricted to the imme- 
diate vicinity of the critical point. In contrast, large fluctua- 
tion regimes typically occur in strongly coupled and/or low- 
dimensional systems. 

Fluctuating order is characterized by a correlation length £ 
and a typical fluctuation energy (or frequency) A. The phys- 
ical properties at distances larger than £ and energies smaller 
than A will be that of the disordered phase. On length scales 
smaller than £, the system is critical (not ordered, as some- 
times implied) - this follows directly from standard scaling 
arguments. However, if the anomalous dimension r\ of the 
order parameter is small, then spectral features of the crit- 
ical regime are not very different from that of the ordered 
phase, e.g., the branch points in the critical spectrum follow a 
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FIG. 3 Schematic phase diagrams, illustrating possible transition scenarios into stripe-ordered states. The vertical axis is temperature, whereas 
the horizontal axis - where x may represent doping - tunes the interplay between superconductivity and spin/charge order, which is assumed to 
be competitive. The solid (dashed) lines are transition into lattice-symmetry-breaking (superconducting) states. Case a) corresponds to a weak- 
coupling scenario, with a single stripe ordering transition, while c) may be realized at strong coupling, with distinct transitions for nematic, 
charge, and spin order. Assuming that increasing x also moves the system from strong to weak coupling, the intermediate case b) is possible 
as well. Effects of lattice pinning and incommensurability are ignored, as is the interplay between stripe/checkerboard and horizontal/diagonal 
order. Quenched disorder will smear out both CDW and nematic transition and likely turn the SDW into a cluster spin glass ( 52). In addition, 
a structural anisotropy will smear out the nematic transition as well, see Sec. III. Gl 



dispersion similar to that of Goldstone modes of the ordered 
state. 12 Close to a quantum critical point, A is typically the 
gap in the collective-mode spectrum, and at low temperatures 
£ and A scale with the distance to criticality according to 
£ oc \r — r c \~ v ', A oc \r — r c \ vz , where r is the tuning pa- 
rameter (e.g. doping), and v and z are the correlation length 
and dynamical exponents of the transition at ?* = r c . 

Fluctuating order is invisible to static probes (unless there 
is some form of pinning, see below). Some probes, like elastic 
neutron scattering or fiSR, are quasi-static, i.e., average over 
a time scale which is large compared to electronic scales, but 
can be comparable to the time scale, 1/ A, of collective fluctu- 
ations. As a result, the "ordering" temperature as determined 
by quasi-static probes will not be unique, but instead depend 
on the type of probe - this is a typical signature of a large 
fluctuation regime. 

A direct probe of fluctuating order is given by the low- 
frequency dynamic susceptibility, u>). The imaginary 
part x" (k, u>) will be strongly peaked at the ordering wavevec- 
tor, k ~ Q, and at uj ~ A (provided sufficient energy resolu- 
tion of the experiment), and the real part x'{Qi w = 0) w iU 
diverge at the critical point. In the spin sector, \" can be 
measured by inelastic neutron scattering; in the charge sec- 
tor, electron energy-loss spectroscopy (EELS) is in principle 
the appropriate method, however, to date its energy resolu- 
tion (> 0.1 eV) is insufficient to detect fluctuating stripes in 
cuprates. Note that the static structure factor S{k), being an 
energy-integrated quantity, is not a suitable probe for fluc- 
tuating order near a quantum phase transition, because it is 
not directly sensitive to the low-energy part of the fluctuation 



The critical exponent 77 characterizes the propagator of order-parameter 
correlations G(k, ui) at criticality. At a quantum critical point with dynam- 
ical exponent z = 1, G takes the form G(k, ui) oc [k 2 — uj 2 } (^ 2 + r i)/ 2 . 



spectrum. 13 (This is different near a classical phase transition, 
where S(k) and cj = 0) contain the same information, 
because the temperature T is much larger than the relevant 
fluctuation frequencies.) 

While all statements here were for a clean system, the pres- 
ence of quenched disorder qualitatively modifies the picture. 
In particular, disorder can induce pinning of otherwise slowly 
fluctuating order, such that it is detectable by static probes. 
This will be discussed in more detail in the next subsection. 



G. Influence of structural effects and disorder 

So far, we have discussed the concepts of "order" and 
"phase transitions" without taking into account the effects of 
structurally broken lattice symmetry and of quenched disor- 
der. These effects often severely complicate (and also some- 
times simplify) the identification of ordered phases (fTlh . 

1. Uni-axial in-plane anisotropy 

For our discussion, the most important structural effect is 
that of a uni-axial lattice anisotropy of the CUO2 plane, which 
breaks the C4 rotation symmetry. 14 While such anisotropies 
are absent in the BSCCO 15 and CCOC compounds, they are 



13 The static structure factor is given by S (k) = J dujS(k, ui) / (2ir), and the 
dynamic structure factor S(k, ui) is related to x"{k, ui) via the fluctuation- 
dissipation theorem, x"(k, w) = (l - e-0" J )S(k,ui)/2, where /3 = 1/T 
is the inverse temperature. 

14 While such an anisotropy can be the result of an interaction-driven sym- 
metry breaking of the correlated electron system, we are here concerned 
with anisotropies of structural origin. 

15 BSCCO displays a structural supermodulation along the (1,1) direction 
of the Cu lattice, with a wavelength of 4.8 unit cells, whose origin and 
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important in YBCO and 214 materials. 

In YBCO, in addition to the CuC<2 planes there exist chain 
layers with CuO chains running parallel to the orthorhombic 
b axis. For doping 6 > 0.4, this results in an orthorhombic 
structure with inequivalent a and b axes (with lattice constants 
a < b). Single crystals commonly are "twinned", i.e., contain 
both orientations of chains, such that macroscopic measure- 
ments average over the inequivalent a and b directions. How- 
ever, it has become possible to produce de-twinned crystals, 
and relevant experimental results will be described below. 

In the order parameter language, the anisotropic structure 
results in a (small) global field coupling linearly to nematic 
and quadratically to CDW and SDW order. Hence a transition 
to a nematic state will be smeared; conceptually, the electronic 
nematic state remains well defined only for a large electronic 
anisotropy in the presence of a small structural anisotropy. In 
contrast, the transition to a stripe state remains sharp, with its 
critical temperature being enhanced. 

In the 214 compounds, various structural modifications 
occur as function of temperature. A transition from 
a high-temperature tetragonal (HTT) structure to a low- 
temperature orthorhombic structure (LTO) occurs at a doping- 
dependent transition temperature between 200 and 500 K. 
This HTT^LTO transition occurs as a result of the bond 
length mismatch between the CuC<2 planes and the La202 bi- 
layers. This mismatch is relieved by a buckling of the CuC<2 
plane and a rotation of the CuOe octahedra. In the LTO phase, 
the crystallographic axes are rotated by 45° w.r.t. those of the 
HTT phase; in this article, we shall use the coordinate nota- 
tion of the HTT phase unless otherwise noted. In terms of 
electronic parameters of the Cu lattice, the LTO phase is char- 
acterized by inequivalent diagonals. 

While LSCO remains in the LTO phase down to lowest 
temperatures, the compounds LBCO, LESCO, and LNSCO 
display an additional low-temperature tetragonal (LTT) phase. 
The LTO-^LTT transition occurs between 50 and 150 K and 
is driven by the smaller radius of e.g. the Nd and Eu ions com- 
pared to La. In the LTT phase, the CuOe octahedra are rotated 
such that now the a and b axes of Cu lattice are inequivalent. 
Therefore, in the LTT phase, there is again a field coupling 
to nematic order in each CUO2 plane, which is relevant for 
stabilizing stripe order. 16 However, the direction of the in- 
plane anisotropy alternates from plane to plane, rendering the 
global crystal symmetry tetragonal and macroscopic in-plane 
anisotropies absent. 

Concrete numbers for the in-plane anisotropies can be ex- 
tracted from first-principles calculations. For the LTT phase 
of 214 compounds, a simple estimate can be obtained from 
the octahedral tilt angles. For a tilt angle of a 4° ... 5° 
in LNSCO- 1/8, the relation t x /t y ~ | cos(7r — 2a) gives a 
hopping anisotropy of about At/t ~ 1% . . . 1.5% i 



In YBCO the situation is more complicated, as both the struc- 
tural distortion and the CuO chains contribute to the (effec- 
tive) hopping parameters in the planes. Recent LDA calcula- 
tions <|57|) indicate t a < t b and At/t ~ 3 . . . 4%. 

In principle, an electronic in-plane anisotropy can also 
be induced by applying uni-axial pressure to an otherwise 
isotropic sample. However, those experiments tend to be diffi- 
cult. A few experimental results are available, demonstrating 
the interplay of lattice distortions and stripes (58). 



2. Quenched disorder 

Most superconducting cuprates are "dirty" materials, in the 
sense that chemical doping in non-stoichiometric composi- 
tion inevitably introduces disorder (exceptions are e.g. the 
oxygen-ordered YBCO-6.5 and YBCO-124 compositions). 17 
The physics of disorder in cuprate superconductors is ex- 
tremely rich and only partially understood; the reader is re- 
ferred to a recent review of both experimental and theoretical 
aspects d59l) . 

For the discussion of ordering phenomena, it is important 
that the disorder potential couples to the charge sector of the 
CUO2 plane. From symmetry, this disorder then is of random- 
mass type for an SDW order parameter, i.e. an impurity at Xq 
acts as ((/^(xq), whereas it is of much stronger random-field 
type 18 for a CDW order parameter, (<fi c (xo)- For a nematic 
order parameter, the coupling is also of random-field type, 
except for pure site disorder which does not locally select a 
direction. 

Ordering in the presence of randomness is a difficult prob- 
lem in statistical mechanics, which is not fully understood 
even in the classical case. Some important questions are about 
(i) the existence of a true ordered phase, associated with a 
sharp ordering transition, (ii) the nature of the phase transi- 
tion, (iii) the existence of regimes with anomalous properties 
near the (putative) phase transition. In the following, I shall 
only touch upon a few aspects relevant to cuprates and the fo- 
cus of this article, and I refer the reader to literature, namely to 
Ref. d60l) for the famous Harris criterion about the stability of 
phase transitions in the presence of randomness, Ref. (16 II) for 
the Imry-Ma argument about the stability of ordered phases, 
Ref. (1621) for an overview on the random-field Ising model, 
and Ref. d63l) for an overview on rare regions near phase tran- 
sitions. 

In the random-mass case, the ordered phase is characterized 
by true symmetry breaking and survives as a distinct phase. 
The phase transition is expected to be sharp, albeit perhaps 
modified compared to the clean case according to the Harris 
criterion, with the exception of certain quantum phase transi- 
tions with order parameter damping (1631) . 



properties are not completely understood {E^tl . 

Consequently, electronic stripe ordering in the LTT phase of LBCO, 
LESCO, and LNSCO is only accompanied by spontaneous breaking of 
translation symmetry, i.e., those stripes do not possess electronic nematic 
order. 



The relevant dopants are typically located away from the CUO2 planes, 
therefore the disorder potential is often assumed to be smooth. 
The crucial difference between random mass and random field cases is that 
a random field breaks the order parameter symmetry whereas the random 
mass does not. 
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In contrast, in the random-field case the effects of quenched 
disorder are more drastic. In low dimensions, the ordered 
phase ceases to exist, because the system breaks up into do- 
mains which are pinned by the local fields ( |6ll) . Consequently, 
there is no sharp phase transition upon cooling, and the corre- 
lation length is finite even in the zero-temperature limit. This 
applies in particular to the random-field Ising model in two 
space dimensions, and is expected to hold for discrete Zm 
symmetries as well. 19 As the system is at its lower critical di- 
mension d~ = 2, the domains are exponentially large for weak 
disorder. In the case of continuous symmetry the lower criti- 
cal dimension is shifted to d~ = 4. Experimentally, random- 
field pinning implies that static probes will see a gradual, i.e. 
smeared, onset of order. Rather little is known about dynam- 
ical properties in the quantum case. On the disordered side 
of a quantum phase transition, randomness will induce in-gap 
spectral weight, and it is conceivable that the smeared phase 
transition is accompanied by slow glass-like order-parameter 
dynamics. 20 

For stripes, the conclusion is that the charge ordering tran- 
sition is generically smeared due to random-field effects (as- 
suming the inter-plane coupling to be weak). The spin or- 
dering transition could still be sharp, however, the coupling 
between spin and charge sectors may invalidate this guess: 
Short-range charge order can induce magnetic frustration in 
the spin sector, leading the spin-glass behavior. Whether this 
setting allows for a sharp phase transition is not known. Ex- 
perimentally, both spin and charge ordering transitions appear 
to be broad and glassy, see Sec. IIII.Fl 

Pinning of multiple equivalent collective modes, e.g. hori- 
zontal and vertical stripes, can lead to interesting phenomena 
d65l l66h . If the clean system would display uni-directional 
order, random-field pinning tends to globally restore the rota- 
tion symmetry, because domains of different orientation will 
occur. Locally, rotation symmetry will still be broken (un- 
less the repulsion between the two order parameters is much 
weaker than the disorder effect). In contrast, if the clean sys- 
tem has a tendency to bi-directional (checkerboard) order, pin- 
ning may locally induce some anisotropy. As a consequence, 
random-field pinning makes it difficult to deduce the nature 
of the clean order (stripes vs. checkerboard), in particular for 
weak order and strong pinning. 21 



III. EXPERIMENTAL EVIDENCE FOR LATTICE 
SYMMETRY BREAKING 

Let me turn to actual results of experiments on cuprate su- 
perconductors. While it seems that detecting a symmetry- 



Zff is the cyclic group of N elements, formed e.g. by rotations about a 
single axis with ./V-fold rotation symmetry. 

It has been demonstrated that stripes arising from frustrated interactions 
g display glass-like behavior even in the absence of quenched disorder 

A detailed analysis of the impurity-induced patterns in a magnet with a 
tendency towards valence-bond order is in Ref. i6T\) . 



breaking order should be a straightforward undertaking, both 
fundamental and practical problem complicate matters. 

The most important fundamental problem is related to the 
effect of quenched disorder, described in Sec. III.GI Disorder 
from chemical doping acts as a random field for the CDW or- 
der parameter, hence it is a relevant perturbation. Thus, the 
charge ordering transition is smeared, most likely eliminating 
thermodynamic singularities. A second fundamental problem 
is that macroscopic manifestations of broken rotation symme- 
try can only be expected when uni-directional stripes in all 
CuC<2 planes are commonly aligned in one direction. Such 
a single-domain situation cannot be expected to be realized, 
mainly because of the random-field effects, unless a global 
symmetry-breaking field exists which selects one stripe direc- 
tion, Sec. III.GI The best candidate here is the YBCO family. 

The direct observation of superstructures is nevertheless 
possible and will be discussed below. Here, practical issues 
such as the availability of sufficiently large single crystals, 
clean surfaces etc. become important, but those issues have 
been resolved at least partially. Scattering experiments which 
average over a large spatial area are again faced with the do- 
main problem: If e.g. domains of both horizontal and vertical 
stripes are simultaneously present, distinguishing this from a 
local superposition of both (i.e. checkerboard order) requires 
a careful analysis. 

As a result of those efforts, signatures of translation sym- 
metry breaking have been found in a variety of hole-doped 
cuprates, most notably in LNSCO, LESCO, and LBCO. The 
phase diagrams of the two latter are shown in Fig. 2] where the 
results from different measurement techniques have been col- 
lected. For the three materials, the order can be consistently 
interpreted in terms of uni-directional SDW and CDW order, 
i.e. stripes, over a wide doping range. Moreover, static incom- 
mensurate SDW order has been established in La2Cu04+a, in 
LSCO for x < 0.13, and for YBCO for 5 < 0.45 (but in the 
latter case the order is only short-ranged and of glassy charac- 
ter). An exciting recent development is the clear observation 
of stripes in STM experiments on BSCCO and CCOC com- 
pounds, although it has to be kept in mind that STM probes the 
surface layer only. As discussed in more detail in Sec. IVI.Al 
stripe signatures seem to weaken with increasing number of 
adjacent CUO2 layers in hole-doped compounds, and have not 
been reported in electron-doped cuprate materials. 

Direct observations of stripe order are corroborated by more 
indirect probes: NMR and /iSR experiments provide evidence 
for inhomogeneous magnetism in 214 compounds, with a 
temperature dependence similar to that seen in scattering ex- 
periments. Besides the temperature and doping dependence 
of the symmetry-breaking order, its dependence on an applied 
magnetic field has been studied as well, which allows to draw 
conclusions about the relation between superconductivity and 
stripe order. 

In principle, translational symmetry breaking should leave 
well-defined traces in the dispersion of all elementary exci- 
tations due to band backfolding and the opening of Bragg 
gaps. However, clear-cut experimental signatures are difficult 
to identify, due to a variety of complications: Disorder and the 
simultaneous presence of horizontal and vertical stripes tend 



11 




0.00 0.05 0.10 0.15 0.20 0.25 0.00 0.05 0.10 0.15 0.20 0.25 

Sr content x Ba content x 

FIG. 4 Phase diagrams of LESCO (left) and LBCO (right). For LESCO, the data points represent the superconducting T c (68), the charge- 
ordering temperature T c h from resonant soft x-ray scattering J69l) . the spin-ordering temperature T sp from neutron scattering d70l) and from 
/iSR i53i) . and the LTO-LTT transition temperature Tltt from x-ray scattering i53) . For LBCO, the data are T c d7 It 1721) (stars, squares), 
Tltt from x-ray scattering < F73l : [74h (crosses,circles) and T c h from x-ray scattering < f73h . and T ap from neutron scattering < r73l ; p74l) (filled 
squares,triangles) and from /iSR (75). The lines are guides to the eye only. The published data display a rather large spread, with transitions 
often being broad, which may be due to effects of disorder and/or sample inhomogeneities. Moreover, data obtained on polycrystals and single 
crystals may differ substantially: For LBCO- 1/8, (iSR results show T sp ~ 29 K for polycrystals d75l) and T sp ~ 40 K for single crystals ( 17a) . 



to smear the signal, and matrix element effects do not allow to 
observe all bands. The experimental results for magnetic ex- 
citations, phonons, and the single-electron spectrum as mea- 
sureed by ARPES will be discussed in subsections below. 

Signatures of rotation symmetry breaking have been most 
clearly identified in underdoped YBCO, both in transport 
and in neutron scattering. In particular, the magnetic exci- 
tation spectrum of YBCO-6.45 was found to develop a spatial 
anisotropy below about 150 K, see Sec. IIII.BT4l Whether these 
data should be interpreted in terms of a Pomeranchuk insta- 
bility of the Fermi surface or in terms of fluctuating stripes is 
open at present and will be discussed in Sec. IV.II Locally bro- 
ken rotation symmetry is clearly visible in STM data obtained 
from the surface of underdoped BSCCO and CCOC, where it 
is accompanied by stripe formation, see Sec. IIII.CI 



A. Static order in neutron and x-ray scattering 

Long-range order accompanied by breaking of lattice trans- 
lation symmetry leads to sharp superlattice Bragg peaks in 
diffraction experiments. Those can be detected by neutrons or 
by x-rays. Experiments require sufficiently large single crys- 
tals, which are by now available for many cuprate families. 

1 . Spin density waves seen by neutron diffraction 

Experimental evidence for static stripe-like order was first 
found in neutron-scattering experiments on LNSCO with dop- 
ing level 0.12 ( fl9t |20T) . which is a superconductor with 
an anomalously low T c of roughly 5K. Those experiments 
detected static spin correlations with an onset temperature 
of about 55 K, which were peaked at wavevectors Q sx — 
2tt(0.5 ± e s ,0.5) and Q sy = 2tt(0.5,0.5 ± e s ), i.e., at four 
spots slightly away from the (tt, it) AF order of the parent an- 



tiferromagnet. At the same time, neutron scattering was used 
to locate the LTO-LTT structural transition at 70 K and the on- 
set of charge order slightly below this temperature, with order- 
ing wavevectors Q cx — 27r(±e c ,0) and Q cy = 2-7r(0, ±e c ). 
Within the experimental accuracy, e s = e c /2 = 0.12, where 
e s = e c /2 is expected on symmetry ground for coexist- 
ing collinear SDW and CDW orders, Sec. III. El This type 
of order was later found in LNSCO also for doping levels 
0.08 < x < 0.20, using both neutron and x-ray scattering 
(I77l;l78l) . The incommensurability roughly follows the doping 
level, e s ss x for x < 0.12, whereas it tends to saturate for 
larger x, with e s « 0.14 at £ = 0.20 (Fig.©. 22 

In fact, incommensurate dynamic correlations were de- 
tected much earlier in LSCO at various doping levels, with 
a similar wavevector dependence on doping (192c 1931) . Sub- 
sequent elastic neutron scattering experiments found quasi- 
static magnetic order, with wavevectors obeying e s w x, in su- 
perconducting LSCO for x < 0.13 (80; 94). This makes clear 
that LSCO is closely located to a QCP associated with incom- 
mensurate SDW order. Indeed, inelastic neutron scattering 
on LSCO-0.14 reported scaling behavior of the magnetic ex- 
citation spectrum, consistent with a nearby QCP (I95I) . Low- 
energy incommensurate fluctuations then occur as precursors 
of static order at these wavevector, see Fig.|6]and Sec. IIII.B72l 

Impurity doping experiments of LSCO, with a few percent 
of Cu atoms replaced by non-magnetic Zn, support the idea of 
nearly ordered stripes in LSCO: Zn induces elastic intensity 
for x — 0.14 d lOOt) . whereas it somewhat broadens the elastic 
peaks at x = 0.12 as compared to the Zn-free sample d87E 



In a picture of charge stripe order, the doping dependence of the incom- 
mensurability, e s (x), can be translated into the doping dependence of both 
stripe distance and stripe filling. This allows an interesting connection to 
the doping dependence of the chemical potential, see Sec. lV.G.H 
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FIG. 5 Summary of experimental data illustrating the doping de- 
pendence of incommensurability e s in the cuprates (the so-called 
"Yamada plot"). Results have been obtained by different groups : 
LNSCO CI; HI) (squares); LSCO {H; S iiiil 
(triangles); La 2 Cu0 4+i Ref. (OS ^> (rhombi); Zn-doped LSCO 
lU m> (circles); YBCO dHHifl HElb (crosses). For LSCO, the 
plot also shows the incommensurability y2e' s of the diagonal mod- 
ulations observed for x < 0.06 (filled triangles). For both YBCO 
and La2Cu04+a, the x values corresponding to the hole doping level 
are estimates only. Note that most data points are slightly below the 
line e s — x (solid). The vertical dashed line is the boundary between 
insulating and superconducting phases in LSCO. 



1881) . These observations support the notion that Zn impurities 
tend to pin stripes, although it should be noted that Zn i s also 
know n to induce magnetic moments in its vicinity d 1 lk Il02t 
Il03l) which certainly contribute to the slowing down of spin 
fluctuations. 

Static incommensurate SDW order, similar to that of 
LNSCO, has also been dete cted in neutron scattering on 
La 2 Cu0 4+(5 <H), LBCO-1/8 (flOl) . and LESCO-0.15 ©. In 
the two latter compounds, the full doping range of static mag- 
netic order has been mapped out by other probes: In LBCO, 
magnetization and /iSR measurements have established mag- 
netic order for 0.095 < x < 0.155 (Slz!), while //SR mea- 
surements give evidence for magnetic order in LESCO over 
the entire doping range up to x = 0.20 d55l) . The La2Cu0 4 +,5 
system, where interstitial oxygens are positioned in every 
fourth La2C<2 layer for 8 w 0.12 (the so-called stage-4 struc- 
ture), displays an incommensurability of e s w 0.12 at a hole 
doping level of roughly 0.15. Remarkable, SDW order and 
bulk superconductivity appear simultaneously at T c w 42 K 
dS6h . However, the behavior of this compound has been at- 
tributed to phase separation into magnetic and superconduct- 
ing domains dl05b . 23 

To fully characterize the SDW order common to the 214 



compounds, more information is required. The first ques- 
tion, triggered by the finding of four magnetic Bragg peaks, 
is whether the order consists of two types of uni-directional 
stripe domains with a single Q vector each, or whether it is 
of checkerboard type with two Q vectors. The experiment did 
not detect magnetic peaks along the diagonal direction, i.e., 
at locations 27r(0.5 ± e',0.5 ± e'); this rules out a checker- 
board with modulation directions along (1,0), (0,1), but could 
be compatible with a checkerboard of diagonal stripes (f88b . 
However, such a structure would lead to diagonal CDW peaks, 
which were not detected. 24 Hence, the most plausible inter- 
pretation is that of two types of large domains (i.e. horizon- 
tal and vertical) in which uni-directional SDW and CDW co- 
exist. In fact, if stripes follow the structural distortion pattern 
of the LTT phase, then the stripe direction can be expected to 
rotate by 90° from plane to plane. 25 

The second question concerns the size and orientation of 
the spin moments in the SDW phase. The size of the or- 
dered moment per Cu site was determined from neutron scat- 
tering in LSCO to vary between 0.04 . . . 0.07 for 0.06 < 
x < 0.1 and to reach a maximum of 0.1 hb for LSCO-0.12 
dH); in LNSCO-0.12 and La 2 Cu0 4+< 5 the neutron scatter- 
ing results yield moment sizes of 0.1 and 0.15 hb ( I20l 1861) . 
In comparison, /iSR experiments suggest a moment size of 
0.3 . . . 0.35 /is, i.e. roughly 60% of the value in the undoped 
parent compound, for LNSCO, LBCO-1/8, LSCO-0.12, and 
La2Cu0 4 +5 Note that the /xSR moment sizes are 

usually inferred from a simulation of theuSR signal taking 
into account the spatial stripe structure d 1 08k 1 1 id) , which con- 
sists of small and large moments, and the quoted numbers de- 
note the maximum (not average) moment size in the stripe 
state. The remaining discrepancy regarding the moment size 
is not fully understood, and disorder in the moment directions 
as well as calibration issues of neutron scattering may play a 
role. 

The moment orientation was studied on LNSCO-0.1 2 us - 
ing polarized neutron scattering by Christensen et al. (II 1 II) . 
The moments are primarily lying in the CUO2 plane, as in 
the undoped parent compounds, consistent with conclusions 
drawn from susceptibility measurements (|72|) . However, the 
neutron result was not entirely conclusive w.r.t. the in-plane 
orientation: The polarization signal was found consistent with 
either a collinear single-Q (i.e. stripe-like) structure or a non- 
collinear two-Q structure, but was inconsistent with a single- 
Q spiral or a collinear checkerboard order. Together with the 
simultaneous existence of charge order, the data are again con- 
sistent with collinear ordering, but a direct proof is still miss- 
ing. 

The third question concerns the correlation length of the 



With an eye towards La2Cu0 4 ^j, Ref. llOd) discussed a scenario of com- 
peting magnetism and superconductivity, where - in the case of phase co- 
existence at low temperatures - both orders may set in at the same temper- 
ature. 



24 See Sec. lII.EI for the Landau theory argument regarding the relative orien- 
tation of spin and charge peaks. 

25 In both LSCO-0.12 and La2Cu0 4+ ^, the quasi-elastic magnetic peaks do 
not exactly lie along the high-symmetry directions of the crystal, but are 
rotated by approximately 3° (8a : 1 1071). T his is consistent with the crystal 
symmetry being LTO instead of LTT <66l) . 
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SDW order, both in-plane and along the c axis. In 214 com- 
pounds, the in-plane order reaches correlation lengths £ a f, of 
200 A and beyond, however, is often not resolution limited 
(whic h points toward disorder effects and/or glassy behavior 
(II 12b ). In contrast, the correlation length £ c in c direction is 
typically only one interplane distance, with the exception of 
the field-induced signal in LSCO-0.10 d 1 1 3h where £ c corre- 
sponds to 6 interplane distances. 

The weak c-axis correlations of stripes can be rational- 
ized considering that, with stripe directions alternating from 
plane to plane, interactions between 2nd-neighbor layers are 
required to align the stripe pattern along the c axis. Such inter- 
actions will be small and have to compete with pinning forces 
from defects - it is no surprise that the latter are dominant. In 
contrast, in the case of field-induced order, "correlated" pin- 
ning by vortex lines can enhance three-dimensional correla- 
tions. 

A fourth question is: How important is the LTT distor- 
tion for thejropejmmce of static stripes? A number of stud- 
ies dl04t ll 14t 1 1 1 5b have been devoted to the compounds 
Lai5/gBai/8_ x Sr x Cu04 with doping level 1/8, where the 
LTT distortion of LBCO-1/8 weakens with Sr doping and 
disappears at x ~ 0.09. These studies suggest that static 
stripe order is tied to the LTT distortion, i.e., it disappears 
in favor of stronger superconductivity around x sw 0.09. 
In contrast, a very recent pressure study on LBCO-1/8 re- 
vealed that stripe domains still occur in the high-pressure 
high-temperature tetragonal (HTT) phase above pressures of 
2GPa dl 16l) . If verified, this would constitute an exciting 
example of simultaneous, electronically driven, spontaneous 
breaking of rotation and translation symmetry. 
Jn the insulating regime of both LSCO dgolS and LNSCO 
1 at very small doping, x < 0.055, a different type of mag- 
netic order has been found. Here, elastic peaks were found 
at Q s = 27r(0.5 ± £^,0.5 ± e^), i.e., in the diagonal direc- 
tion of the Cu square lattice. The incommensurability follows 
e' s w x/V2 for 0.02 < x < 0.055; thus the peak distance 
to (tt, 7t) appears to follow the same linear x dependence as 
for superconducting compounds with 0.055 < x < 0.125. 
This may suggest a common origin of the SDW in both dop- 
ing ranges, and consequently, the low-doping SDW order has 
been interpreted in terms of diagonal stripes. However, no ev- 
idence for charge order has been reported for x < 0.055. An 
alternative scenario consists of spi ral SDW order without ac- 
companying CDW ([OS [HI (HI)- To my knowledge, polar- 
ized neutron scattering - which would be able to distinguish 
the two alternatives - has not been performed to date. 

In other cuprate families, static incommensurate spin or- 
der has not been detected, with the exception of strongly un- 
derdoped YBCO: A conclusive set of data stems from de- 
twinned crystals of YBCO-6.35 and 6.45 where neutron scat- 
tering measurements of Hinkov et al. and Haug et al. 
J96b detected incommensurate quasi-static order at wavevec- 
tors 27r(0.5 ± e s , 0.5). However, the data on YBCO-6.45 sug- 
gest the order is weak and glassy: Quasi-static order in neu- 
tron scattering sets in below 30 K, with a correlation length 
of about 20 A only, and ^SR measurements detect static or- 
der only at 1.5 K ©. The findings of Refs. dH |H) are 
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FIG. 6 Comparison of constant-energy scans at 3 meV through an 
incommensurate magnetic peak (along path shown in inset) for (a) 
Lai.85Sro.i5Cu04 and (b) Lai.48Ndo.4Sro.i2Cu04. Both scans are 
at T = 40 K > T c and illustrate the similarity of low-energy fluctua- 
tion in non-stripe-ordered (a) and stripe-ordered (b) 214 compounds 
(reprinted with permission from Ref. (112), copyright 1999 by the 
American Physical Society). 



broadly consistent with earlier neutron scattering reports us- 
ing twinned crystals @ |H IH). While Mook et al. © 
suggested the presence of stripe order in YBCO-6.35, Stock 
et al. 1981) reported the existence of a "central mode" 
in quasi-elastic scattering, centered at (tt, tt) and correspond- 
ing to glassy short-range order. Note that the effect of twin- 
ning on the in-plane geometry has been nicely demonstrated 
in Ref. (91), i.e., twinning smears an incommensurate signal 
such that it appears as a single peak at (tt, tt). In a stripe pic- 
ture, the incommensurate SDW order in underdoped YBCO 
would correspond to stripes along the b axis. In YBCO-6.45, 
e s w 0.045, while the nominal hole doping is x 0.085. 
These values are inconsistent with the relation e s w x, estab- 
lished for stripe order in LSCO with x < 1/8. 

Finally, we briefly mention the neutron-scattering search 
for commensurate (tt, tt) antiferromagnetic order in super- 
conducting cuprates, which was primarily motivated by the 
proposal of staggered loop-current order (dubbed rf-density- 
wave) in the pseudogap regime There were several 

observations of AF order in the underdoped YBCO samples 
(EH [H2 [HI), with ordered moments of about 0.02. . . 0.05 
fis- However, the polarization neutron analysis revealed that 
the order is dominated by moments aligned in the CUO2 plane, 
as in the insulating parent compound. Considering, on the 
one hand, that other high-quality single crystals with a sim- 
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ilar doping level do not show a similar order (1 123b and, on 
the other hand, that impurity substitution can induce anti fer- 
romagnetic order at 300 K even at optimal doping dl24l) . it 
is likely that the observed AF order is not a generic property 
of the underdoped state. Subsequently, Mook et al. d 125b re- 
ported polarized neutron experiments suggesting the existence 
of a weak AF quasi-2d order with moments perpendicular to 
the CuC>2 plane. However, the estimated magnitude of the or- 
dered moment is 0.0025 /is, which is close to the experimen- 
tal threshold of detection ( 125). No subsequent experiments 
have reported conclusive evidence for magnetic order at or 
near (ir, ir) setting in at the pseudogap scale. 



2. Charge density waves seen by x-ray scattering 

Although the picture of stripe order, i.e., coexisting SDW 
and CDW, was already proposed in the context of the first 
magnetic neutron scattering results in 1995, the first unam- 
biguous observation of charge order was only made by Abba- 
monte et al. in 2005 dl32l) . The primary reason is that neu- 
tron and non-resonant x-ray scattering can only detect charge 
order indirectly by the associated lattice distortion, as these 
techniques are mainly sensitive to the nuclear scattering and 
the core electron scattering, respectively. 

Nevertheless, results from neutron and non-resonant x-ray 
scattering have provided valuable information on charge or- 
dering, because lattice distortions can be expected to follow 
charge order regarding both the amplitude and the tempera- 
ture dependence. In particular, high-energy x-rays allow the 
use of a synchrotron source to obtain high intensity and high 
momentum resolution dl30h . Using these techniques, charge- 
order superlattice peaks were found in LNSCO for doping lev- 
els 0.08 < x < 0.20 (US HE EU and recently in 
LBCO-1/8 (l73ETT27h. 

In this context, studies using extended x-ray absorption fine 
structure (EXAFS) are worth mentioning, which provide in- 
formation on local lattice displacememnts. In LNSCO-0.12, 
lattice fluctuations have been found t o stro ngly increase be- 
low the charge-ordering tempera ture dl28l) . A similar effect 
was found before in LSCO-0. 15 dl29h . supporting the idea of 
dynamic stripes in this material. 

The most direct information on the charge modulation in 
cuprates can be obtained by resonant soft x-ray scattering us- 
ing photon energies at the O K and the Cu L edge d 1 3 lh . This 
technique was applied to stripe-ordered LBCO-1/8 by Abba- 
monte et al. (|132b and to LESCO at various doping levels by 
Fink et al. d69b. 

In LBCO-1/8, charge superlattice peaks were detected at 
in-plane wavevector 27r(0.25 ± 0.02, 0) with an in-plane cor- 
relation length £ a t, of about 480 A (=125 lattice spacings). 
The dependence on vertical momentum was consistent with 
a period-2 order in c-axis direction, with a small correlation 
length of £ c < 2 interplane distances (which is consistent with 
the x-ray study of Ref. dl27b ). The onset temperature of the 
signal was T c h ~ 60 K which almost coincides with the tem- 
perature of the LTO-LTT transition in this sample. The form 
factor of the scattering signal was used to extract an estimate 



of the actual modulation amplitude in the charge sector. As- 
suming the orbital charge-order pattern obtained in a slave- 
boson analysis of stripes in the three-band Hubbard model 
(EH, the authors obtained a large modulation of oxygen hole 
densities varying between 0.03 and 0. 12, i.e., by a factor of 4, 
within a unit cell. However, it has to be emphasized that this 
analysis is strongly model-dependent, and the assumed mod- 
ulation pattern, although bond-centered, appears inconsistent 
with that obtained in STM d37h . 

Recently, resonant soft x-ray scattering was also performed 
on LESCO (69), with some results shown in Fig. [7] For dop- 
ing x = 1/8, the order appeared at T c h = 80 ± 10 K at a 
wavevector Q c — 27r(0.228, 0), whereas for x = 0.15 T c h 
was 70±10KandQ c = 2tt(0.254, 0). The correlation length 
£ a b was about 80. . . 100 lattice spacings, while £ c was of order 
unity. In contrast to LBCO, the LTO-LTT transition in LESCO 
occurs at significantly higher temperature (Tltt ~ 125 K), 
while spin ordering has been detected at 25 K (45 K) by /iSR 
(neutrons) u55c 1701) . Thus, LESCO displays a sequence of 
well-separated phase transitions with T sp < T c h < Tltt- 

Taken together, the neutron and x-ray experiments estab- 
lish a number of important characteristics of the order in 214 
cuprates: (i) CDW order sets in at higher temperatures than 
SDW order, and both coexist at low T. Thus, CDW order 
is unlikely to be only a subleading consequence of collinear 
SDW order in the sense of Eq. ( fTQb . (Of course, this does 
not exclude that stripe formation is driven by antiferromag- 
netic exchange.) (ii) The ordering wavevectors are related, 
e s = e c /2 within error bars, (iii) The wavevector dependence 
on doping e c (x) is inconsistent with a simple nesting scenario 
of CDW formation, because e c increases with x whereas the 
distance between the antinodal Fermi surface decreases (see 
e.g. Ref. dUi). 

It should be noted that no signatures of charge order have 
been reported for LSCO, even in the doping regime x < 0.13 
where the SDW order appears static. Also, resonant x-ra y 
scattering failed to detect charge order in CCOC-1/8 (1134b . 
where STM has established the existence of a period-4 charge 
ordering pattern d37b 1 1 35D . Here, two explanations are possi- 
ble: Either the STM pattern exists at the surface only, or the 
correlation length of the charge order is too small to be de- 
tectable by x-rays at present. Indeed, the intensity of the x-ray 
signal scales with £^ b , thus assuming £ a b of 10 lattice spacings 
(137b renders the signal 100 times smaller than for LBCO-1/8 
(for comparable CDW amplitudes) which is below the sensi- 
tivity limit of the experiment. 

3. Magnetic order at Q — 

Over the last few years, indications for a distinct type of or- 
der have been found in several experiments, which were mo- 
tivated by the proposal d3~Tl |32|) for the pseudogap phase in 
terms of a spontaneous loop-current order within the unit cell 
of CUO2 planes, with ordering wavevector Q = 0. 

The most direct indication for circulating-current order 
comes from recen t elastic polarized neutron scattering exper- 
iments d~L37t 1 1 38b Il39l) . Such measurements are difficult, as 
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FIG. 7 Resonant soft x-ray scattering results on LESCO (reprinted with permission from Ref. ( l69r) . copyright 2009 by the American Physical 
Society). Left: Temperature dependence of h scans along (h, 0, 1) showing superstructure reflections of LESCO-1/8 and LESCO-0.15 using 
O K (I = 0.75) and Cu L3 (I = 1.6) photon energies. The curves are vertically shifted for clarity; the solid lines are fits to the data. Right: 
Temperature dependence of the intensities of the superstructure reflections, normalized to the intensity at T — 6 K. The dotted line is a 
s/Tch — T fit. The estimated spin-ordering temperatures are marked by arrows J55H70I) . 



the magnetic signal is located on top of a large nuclear Bragg 
peak. Therefore, the magnetic Bragg peak 26 has to be ob- 
tained as a difference between spin-flip and non-spin-flip neu- 
tron scattering signals. As the spin-flip ratio is an unknown 
constant, the subtraction is done by suitable re-scaling of one 
of the signals such that the difference at high temperatures is 
zero. While the initial experiments on YBCO at doping lev- 
els 8 = 0.5 . . . 0.75 were somewhat controversial, the exper- 
iment was repeated on a larger sample of YBCO-6.6 (1 1 3 8h . 
with consistent results: The data suggest the magnetic order 
at wavevector (0, 0) at a doping-dependent temperature which 
varies from 300 K for YBCO-6.5 to 170 K for YBCO-6.75; 
these values appear to match the accepted pseudogap temper- 
atures T* for these samples. (Note that the subtraction pro- 
cedure renders the determination of a sharp onset tempera- 
ture difficult.) Recently, a similar neutron scattering signature 
of magnetic order in the pseudogap state was obtained in the 
tetragonal single-layer HgBa2Cu04+,5 compound (Il39l) . 

In all cases, the moment amplitude is of order 0.1/zs, and 
the moments are oriented roughly in a 45° angle w.r.t. the 



The actual measurement is done at finite wavevectors (2n, 0) or (An, 0). 



planes. As susceptibility measurements appear to exclude 
ferromagnetic order of this magnitude, the ordered moments 
within a unit cell apparently compensate each other. One pos- 
sibility is spin moments on oxygen atoms with opposite direc- 
tions, another one is given by the loop-current order of Varma 
dUl l32t) . In the latter case, however, the moment directions 
suggest that the current loops involve out-of -plane oxygen or- 
bitals. 

A number of other experiments are worth mentioning: 
Early ARPES experiments on BSCCO using circularly po- 
larized photons reported a dichroic signal indicating time- 
reversal symmetry breaking in the pseudo-gap state (Il40b 
(however, at the time, the result was questioned by others). 
A search for static fields in the pseudogap regime using /iSR 
of LSCO gave a null result (11421) . Static screening of the muon 
charge, leading to a local change in the doping level, has been 
invoked to explain the absence of a signal ( 1 43b . However, 
a null result was also reported from NMR measurements in 
the pseudogap regime of Y2Ba4CuyOi5_5 (144). Very re- 
cently, extremely sensitive polar Kerr effect measurements 
detected signatures of broken time reversal in the pseudogap 
regime in a series of underdoped YBCO crystals (114ll) . The 
straightforward interpretation of the Kerr signal is in terms 
of /erromagnetic order, however, with a tiny magnitude of 
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1CP 5 hb per Cu atom. A puzzling aspect of the Kerr mea- 
surements is that memory effects have been found to survive 
up to room temperature far above T*. At present, the relation 
of the Kerr effect to the neutron observations is not settled. 
Finally, careful measurements of the magnetic susceptibility 
in underdoped YBCO revealed a small kink in samples with 

0. 4 < 5 < 0.8 (fl4l . The kink may signify a thermody- 
namic phase transition and occurs at a doping-dependent tem- 
perature which appears to match the onset temperature of the 
described neutron signal. 

Although these results are exciting, further experiments are 
needed to check whether loop-current order is indeed a com- 
mon feature of high-temperature superconducting cuprates, 
and also to establish the relation between different possible 
experimental signatures. A brief theoretical discussion of 
loop-current order is in Sec. lVJI 

B. Inelastic neutron scattering 

High-resolution inelastic neutron scattering, used to probe 
the spectrum of magnetic fluctuations as well as of phonons, 
has been performed extensively on 214 cuprates, on YBCO, 
and, to a lesser extend, on BSCCO. For other materials, the 
lack of sufficiently large single crystals limits the available 
data. 

While the most direct oberservation of translation sym- 
metry breaking is via Bragg peaks in elastic scattering, the 
spectrum of finite-energy excitations can provide information 
complementary to that of elastic scattering probes, such as 
the energy range and character of the fluctuations. Moreover, 
the presence of excitations at very low energies and specific 
wavevectors is usually a precursor of an ordered state, i.e., in- 
commensurate low-energy spin excitations will occur close to 
stripe ordering. 

Inelastic scattering is also a suitable probe for rotation sym- 
metry breaking: An anisotropic fluctuation spectrum must 
arise from an anisotropic state, as the local excitation created 
by the external perturbation cannot change the symmetry. As 
explained above, such a situation can only be expected in the 
presence of a (structural) anisotropy field, as otherwise both 
types of domains will be present with equal weight. The can- 
didate material is YBCO, where indeed signatures of rotation 
symmetry breaking have been identified, see Sec. IIII.B.4I 

1 . Magnetic excitations of ordered stripe phases 

For static spin stripe order, the existence of low-energy 
magnetic excitations (spin waves) follows from the Gold- 
stone theorem. 27 Low-energy incommensurate excitations 
have been probed in a number of stripe-ordered 214 cuprates, 
but extensive data over a wide range of energies are available 
only for LBCO-1/8. 



Magnetic anisotropics will induce a small gap for lattice-pinned commen- 
surate order. 



Such data, taken by Tranquada et al. are shown in 

Fig. [8] Remarkably, rather well-defined collective excitations 
are visible at all energies up to 200 me V. At low energies, four 
spots near (tt ± tt/4, tt) and (tt, tt ± 7r/4), i.e. the ordering 
wavevector, are observed. The spots are found to disperse 
towards (tt, tt) with increasing energy, where they meet at 
around 50 meV. This is essentially consistent with spin-wave 
theory, although the expected spin-wave cones (i.e. intense 
ellipses in constant-energy scans) are never observed, possi- 
bly due to a combination of broadening and matrix element 
effects. Interestingly, the spectrum at elevated energies does 
not appear to follow a simple spin-wave dispersion. Instead, 
the excitation branch above 50 me V is very similar to that of a 
two-l eg spin ladder with an exchange constant J w 100 meV 
(fl46l) . Constant-energy cuts show four intense spots in di- 
agonal direction from (tt, tt), i.e., the scattering pattern has 
rotated by 45° from low to high energies. At an energy of 
200 me V, the experimental intensity distribution has reached 
the boundary of the magnetic Brillouin zone, but at the same 
time becomes damped rather strongly. The location of the in- 
tensity maxima in constant-energy cuts, shown in Fig. [8] trace 
out an "hour-glass" dispersion, Fig. [9] 

As will be discussed in detail in Sec. IV.C.ll the low- 
temperature excitation spectrum of LBCO-1/8 can be nicely 
described by simple models of coupled spin ladders dl47l : l 148c 
Il49l Il50l) . provided that the response of horizontal and ver- 
tical stripes is summed up, Fig. [8] In the absence of perfect 
charge order, a model of fluctuating (or disordered) stripes ap- 
pears more appropriate dl5ll) . which can account for the data 
as well. 



2. Incommensurate spin excitations and hour-glass spectrum 

Investigations of the spin excitation spectrum have played a 
prominent role in cuprate research, mainly because spin fluc- 
tuations are a candidate for the glue that binds the cooper 
pairs. 

Historically, the observation d 152b of the so-called "reso- 
nance mode" in optimally doped YBCO below T c , located at 
41 meV and Q = (tt,tt), triggered enormous activities. The 
doping dep endence of the resonance energy was mapped out 
(Il53l Il54l). and a s imilar resonance was also found in other 
cuprates (Il55l Il56l) . On the other hand, the 214 family of 
cuprates displayed low-energy excitations at incommensurate 
wavevectors (I^l80t [8lh. As a result, it was believed that the 
underlying magnetism is very different, and that the behavior 
of 214 compounds is rather special. 

This view was challenge d by a series of de t ailed neu- 
tron scattering experiments (fl46l [l57t (TM EU [l60l fl6lt 
Il62t 1 1 63b . which mapped out the spin excitations of vari- 
ous, mainly underdoped, cuprates over a wide energy range. 
For YBCO, dispersive excitations were found which em- 
anate from the resonan ce peak and disperse both upwards and 
downwards in energy (fl57l [l58t [l6lt [l62t [l63T) . For stripe- 
ordered LBCO, the incommensurate low-energy excitations 
were found to merge at (ir, tt) at 50meV, and an upper exci- 
tation branch emerges which is well described by the spec- 
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FIG. 8 Left: Neutron-scattering results from ordered stripes in LBCO-1/8 (from Ref. dl46l) . reprinted by permission from Macmillan Pub- 
lishers Ltd: Nature 429, 534, copyright (2004)). Shown are constant-energy cuts through the magnetic excitation spectrum x"(<7j w )- Energy 
has been integrated over the ranges indicated by the error bars. Panels a-c were measured with an incident neutron energy of 80 meV, panels 
d-g with 240 meV, and panel h with 500 meV. Right: Theoretical result from a model of coupled spin ladders (reprinted with permission from 
Ref. ( 1471) . copyright 1999 by the Ame rican Physical Society), as described in Sec. IV.CI The signals from horizontal and vertical stripes have 
been added. (Results obtained in Refs. ( 148; 149; 150) are very similar.) Note that only the reduced Brillouin zone of (n, n) magnetic order is 
shown. 



trum of a spin ladder (fl46h . as discussed above. For LSCO 
a similar excitation structure was found, with incommensu- 
rate low-energy excitations dispersing to wards (tt , tt) a nd dis- 
persing outwards again above 50 me V 

HElSEll)- Very 
recently, signatures of an hour-glass spectrum were reported 
as well for optimally doped BSCCO, although the features in 
constant-energy cuts are relatively broad dl65l) . 

Taken together, these experiments provide evidence for 
an excitation spectrum of hour-glass form being common to 
many cuprates. Indeed, plotting the scattering intensity max- 
ima as function of momentum and re-scaled energy, i.e., rel- 
ative to the exchange constant of parent compound, results 
in Fig. [9] which suggests universality of the hour-glass spec- 
trum. In addition, it shows that the relevant energy scale for 
magnetic fluctuations is J. 

This prompts the question for a common microscopic origin 
of the excitation spectrum. While a weak-coupling descrip- 
tion based on RPA captures some features of the YBCO data, 
a picture of stripes is clearly appropriate for 214 compounds 
where low-energy incommensurate excitations are clear pre- 
cursors of ordered stripes, Fig. [6] It has been proposed early 
on that dynamic stripes are responsible for incommensurate 



excitations below the resonance energy in YBCO as wel l 
(I167t 1 1681) . but opposite views were also put forward (11571) . 
I shall return to this discussion in Sec. lV.C.2l 

Despite the similarities there are, however, a number of im- 
portant differences between the cuprate families which should 
be mentioned: (i) While YBCO displays a sizeable spin gap 
at low temperatures except for small dopings below i5 < 0.5, 
the spin gap in 214 materials is small, (ii) In YBCO, constant- 
energy cuts at the resonance energy show a very sharp peak at 
(tt, tt) (EH EH EH EH fl69l) . In contrast, in LSCO the 
peak is bro ader and well fit by the sum of two incommensu- 
rate peaks (EHEH- (iii) In YBCO, again with the exception 
of small dopings below S < 0.5, the spectra above and be- 
low T c are significantly different. At optimal doping, a clear 
spin gap opens below T c together with the resonance appear- 
ing, while most of the structured response disappears above 
T c . In both YBCO-6.6 and YBCO-6.95, the lower branch of 
the hour gl ass i s only present below T c , but strongly smeared 
above T c dl69t \VKX> . In contrast, in LSCO and LBCO the 
spectral changes at T c are weaker dl46t Il64t 1 1 7 it EH), (iv) 
To my knowledge, an hour-glass dispersion has not been mea- 
sured in other cuprates, with the exception of optimally doped 
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FIG. 9 Universal hour-glass disperison of magnetic excitations in 
cuprates (Fig. 6.3 of Ref. dl661) . reproduced with kind permission 
of Springer Science and Business Media). Comparison of measured 
dispersions along Q = 27r(0.5 + h, 0.5) in LSCO-0.10 (up trian- 
gles) and LSCO-0.16 (down triangles) from Christensen etal. il59h . 
in LBCO-1/8 (filled circles) from (HI), and in YBCO-6.5 (squares) 
fro m Sto ck et al. ( 163) and YBCO-6.6 (diamonds) from Hayden et 
al. dl58l) . The energy has been scaled by the superexchange energy 
J for the appropriate parent insulator. For YBCO-6.6, the data at 
higher energies were fit along the [1,1] direction; the doubled sym- 
bols with bars indicate two different ways of interpolating the results 
for the [1,0] direction. The upwardly-dispersing dashed curve cor- 
responds to the result for a two-leg spin ladder, with an effective 
superexchange of ~ | J; the downward curve is a guide to the eye. 



BSCCO dl65l) . mainly due to the lack of large single crystals 
required for neutron scattering. 

It should be pointed out that the approximate universality 
of the hour-glass spectrum, when plotted as in Fig. [9] only ap- 
plies to doping levels between 10% and optimal doping. For 
smaller dopings, both the incommensurability e s and the en- 
ergy E CIOSS , where upward and downward dispersing branches 
meet, decrease. In LSCO for x < 1/8, both vary roughly 
linearly with doping, and £? cross /400 meV w x has been re- 
ported. Interestingly, this also includes the insulating small- 
doping phase of LSCO, which is characterized by static order 
with diagonal incommensurate peaks and displays an hour- 
glass-like spectrum as well d 1 73b . 

For LSCO at large doping, incommensurate low-energy 
spin excitations have been found to survive throughout the 
entire superconducting phase, disappearing only at x = 0.30 
where the sample becomes non-superconducting (11741) . In this 
overdoped regime, the wavevector-integrated magnetic inten- 
sity below 60meV drops with doping, and almost vanishes 
at x = 0.30 as well dl75l) . These findings point to an in- 
triguing relation between magnetism and superconductivity. 
Moreover, they strongly argue against simple Fermi surface 
nesting as source of magnetic fluctuations even in overdoped 
cuprates, because the magnetic intensity drops within a small 
doping window despite the Fermi surface becoming more well 
defined with increasing x (Il75l) . 



3. Phonon anomalies 

Charge order naturally couples to the lattice: Static charge 
order will lead to periodic atomic displacements and cause 
backfolding of the phonon branches. One can expect that 
slowly fluctuating charge order induces corresponding precur- 
sors, which should be manifest in anomalies in the phonon 
dispersion and linewidth near the charge ordering wavevec- 
tor. Hence, strong stripe-related effects should be visible in 
a momentum scan of the phonon dispersion perpendicular 
to the stripe direction, with a strength dictated by the over- 
all electron-phonon coupling and by matrix elements which 
depend on the specifics of the charge order and the phonon 
mode. 

While various reports on phonon anomalies are in the liter- 
ature, it has not been conclusively established which of them 
are related to stripe order. Essentially all studies focused on 
optical phonon branches, which exist in the energy range of 
60 . . . 100 meV and can be easily distinguished from other ex- 
citations seen in neutron scattering. 

For LSCO at x = 0.15, McQueeney et al. (Il79l) reported 
a discontinuity in the longitudinal optical (LO) dispersion 
branch at 10K, which was interpreted as evide nce for unit 
cell doubling. However, a subsequent study (180) of the same 
compound arrived at a different conclusion: The data speak 
in favor of a continuous dispersion of this bond-stretching 
mode. In addition, the strongest broadening was observed 
near wavevectors (0.57T . . . 0.67T, 0), which can be taken as a 
precursor to stripe ordering with a periodicity of 4 ... 3 lattice 
spacings. The softening of the LO phonon at (tt/2, 0) is in fact 
found in LSCO over a wide doping range, including strongly 
overdoped samples at x — 0.29 d 1 8 lb . suggesting that it is 
unrelated to translational symmetry breaking tendencies. For 
YBCO, conflictin g results on phonon anomalies have been re- 
ported (1 1671 : Il82l) . their relation to stripe order being unclear 
at present. 

A comprehensive neutron-scattering analysis of bond- 
stretching phonons i n str i pe-or dered cuprates was performed 
by Reznik et al. OH fl77h . Both LBCO and LNSCO 
at doping x = 1/8, known to have static stripe order, dis- 
played strong anomalies near wavevector (tt/2, 0). The broad 
lineshapes were interpreted in terms of two (instead of one) 
phonon branches. A subsequent x-ray scattering study with 
higher resolution dl78l) essentially confirmed the broad line- 
shapes, but showed that the two-peak interpretation is not jus- 
tified. The observed phonon anomaly is strongest at the lowest 
T of 10 K, and the phonon lines sharpen with increasing tem- 
perature. The authors extended the measurements to LSCO at 
doping levels x = 0.07, 0. 15, and 0.3, and found broad lines at 
(•7r/2, 0) for the x = 0.07 and 0.15 samples, while narrow lines 
where found in the non-superconducting x = 0.30 compound. 

It appears plausible to associate the bond-stretching phonon 
anomalies with the tendency toward stripe order. A theoreti- 
cal study of the phonon dynamics in the presence of static 
stripes d 1 83b seem to support this assertion, although the in- 
terplay of fluctuating stripes and phonons has not been in- 
vestigated theoretically. There are , however, serious caveats 
with this interpretation: Static stripes should cause multiple 
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phonon branches due to backfolding, which are not clearly ob- 
served experimentally (this could be related to matrix-element 
effects). More importantly, the stripe ordering wave vector in 
the 214 compounds is known to vary with doping, both in the 
spin and charge sector (see Fig. [5]), but the phonon anomalies 
do not show a similarly large shift in momentum space . Al so, 
a somewhat similar LO phonon anomaly was found dl84l) in 
optimally doped YBCO-6.95, where there is otherwise little 
evidence for stripe behavior. 

An interesting alternative interpretation dl84l Il85l) of the 
phonon anomalies invokes one-dimensional physics: It has 
been suggested that the observed anomalies are due to 2kp 
effects of essentially one-dimensional metallic stripes, and 
hence occur in a direction parallel to the stripes. This scenario 
would explain the weak doping dependence of the anomaly 
wavevector, as, in such a one-dimensional picture, kp within 
a stripe is expected to vary weakly with x for x < 1/8. It 
remains open why no backfolding effects are observed. 

In the context of phonons coupling to stripes, thermal con- 
ductivity measurements may provide additional information. 
The phonon thermal conductivity of LSCO was found to be 
strongly suppressed at low temperature, with the suppression 
being correlated with superconductivity dl 86b . Remarkably, 
the suppression was absent in non-superconducting stripe- 
ordered LNSCO and LESCO. While these findings were orig- 
inally interpreted d 1861) as evidence for phonon scattering off 
fluctuating stripes, a careful re-analysis d 1 871) showed that soft 
phonons, caused by the structural instability of the LTO phase, 
provide a scattering mechanism which can account for the ob- 
served suppression of heat transport. Conceptually, a separa- 
tion of soft-phonon-induced and soft-stripe-induced scattering 
is difficult. 

As there can be no doubt about a strong coupling of stripes 
to the lattice, as evidenced e.g. by the large isotope effect 
J68h in stripe-ordered LESCO, more experimental and theo- 
retical work is required to elucidate the interplay of stripes 
and phonons. 



4. Anisotropic magnetic spectra and nematic order 

A state with broken rotational symmetry causes a neu- 
tron scattering intensity XsGfj^) w i tn directional q space 
anisotropy, and vice versa. However, such a spectral asym- 
metry is wiped out if domains of different orientation co-exist 
and are probed by the neutron beam. This complication can 
be avoided by a small symmetry-breaking field. Among the 
cuprates, de-twinned crystals of YBCO have precisely this 
property: For dopings 6 > 0.4 the presence of the CuO chains 
induces a structural orthorhombic distortion, which should be 
able to align nematic domains. 

In a remarkable experimental effort, the spin-fluctuation 
spectrum of de-twinned superconducting YB CO crystal s has 
been studied in detail by Hinkov et al. (HH 1 1 571: Il69h . In 
YBCO, one expects bilayer splitting of all magnetic modes; 
all the following information apply to magnetic excitations 
which are odd under bilayer exchange unless otherwise noted. 

In moderately underdoped YBCO-6.85 and YBCO-6.6 the 



spin excitations are gapped. Below T c , the spectrum is con- 
sistent with the hour-glass shape described above. While its 
high-energy part above the resonance energy approximately 
obeys the square-lattice symmetry, the low- ener g y pa rt is 
significantly anisotropic in both compounds d 1 571 : Il69t) : In 
YBCO-6.6, the intensity at 33meV displays strong incom- 
mensurate peaks along the a axis, whe reas the corresponding 
peaks along the b axis are weak dl69h - this anisotropy may 
be consistent with nematic behavior. However, as function of 
temperature, the spectra change smoothly except at the super- 
conducting T c , i.e., no signature of a nematic order ing transi- 
tion has been detected. The results of Hinkov et al. 
are consistent with earlier reports of anisotropic spin fluctua- 
tions at 24meV on a partially de-twinned YBCO-6.6 crystal 
dl67h . 

More intriguing is the behavior in strongly underdoped 
YBCO-6.45 © with T c = 35 K: At this doping concen- 
tration, the neutron spectrum is essentially gapless at 40 K 
and below. At low energies of e.g. 3 and 7 meV, the inten- 
sity distribution in x"(Qi w ) takes the form of an ellipse in 
momentum space around (it, w), i.e., the intensity is broadly 
distributed along the a axis while it is less broad along the 
b axis, Fig. [10] At energies of 3 meV and below, the in- 
tensity along the a axis is well be fitted by the sum of two 
Lorentzians at incommensurate wavevectors 27r(0.5±e s , 0.5). 
Plotting the incommensurability e s at 3 meV as function of 
temperature shows an order-parameter-like behavior, i.e., it 
decreases with increasing temperature and vanishes at around 
150 K, i.e., far above the superconducting T c , Fig. [10] Al- 
though the orthorhombic structural distortion is expected to 
smear out a sharp nematic phase transition, a reasonable inter- 
pretation of the data is in terms of a spontaneous onset of both 
magnetic incommensurability and magnetic anisotropy of the 
electronic system at around 150 K. This is suggestive of a ne- 
matic phase transition. The strong increase of the in-plane 
transport anisotropy below 200 K (188) of a YBCO sample 
of similar doping level appears consistent with this interpre- 
tation, see Sec. IIII.D.ll (Whether the NMR/NQR work on 
YBCO of Ref. ( 189), showing two inequivalent planar O sites 
but only a single type of Cu site, is related to nematic order is 
not known.) 

The quasi-static neutron signal shows a significant upturn 
below 30 K, but, according to /iSR, static magnetic order sets 
in only below 2 K, with an ordered moment of about 0.05 hb- 
This is consistent with very slow, perhaps glassy, spin dynam- 
ics at low temperatures (I91I) . 

Interpreting the incommensurate magnetism in YBCO-6.45 
in terms of stripes, then the stripes run along the b direction, 
which appears consistent with the resistivity anisotropy, p a > 
Pb. However, as noted above, reports on static charge order in 
underdoped YBCO are controversial. 

The anisotropics of the spin fluctuation spectra in YBCO 
suggest the existence of a nematic QCP around S c = 0.5 J9ll) 
(although one cannot exclude the possibility that nematic tran- 
sitions also exist at larger S, but with a strong smearing due to 
the large orthorhombicity). A magnetic QCP may exist at the 
same or a slightly smaller doping, but static magnetic order 
for S < S c and zero field (284) appears to be restricted to ex- 
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FIG. 10 Signatures of nematic order in YBCO-6.45 (from Ref. (91), Science 319, 597 (2008), reprinted with permission from AAAS). Left: 
Energy and temperature evolution of the a-b anisotropy of the spin correlations. Full squares and empty circles represent data points measured 
at fixed K along a and at fixed H along b, respectively. Scans are normalized and background-subtracted. Solid lines represent the results 
of fits with one or two Gaussians. Right: Temperature and energy evolution of parameters characterizing the spin excitation spectrum. The 
parameters are the results of fits to the raw data. Shown are incommensurability S = e s (red symbols), half-width-at-half-maximum of the 
incommensurate peaks along a , black symbols) and along b i^ 1 , open blue symbols) in reciprocal lattice units, measured at 3 meV. 



tremely low temperatures, while nematic order extends well 
into the pseudogap regime. 



C. Scanning tunneling microscopy 

The techniques of scanning tunneling microscopy and spec- 
troscopy (STM/STS) have contributed enormously to the ex- 
citing progress in the field of cuprates over the past decade 
J27l) . With the caveat of being only sensitive to the physics 
of the sample surface, STM has provided real-space images 
of apparently intrinsic inhomogeneities in BSCCO, allowed 
for a detailed analysis of local impurity physics, and unrav- 
eled the tendency towards ordering phenomena accompanied 
by lattice symmetry breaking - the latter shall be summarized 
here. As STM is essentially a static measurement, "fluctuat- 
ing" stripes can only be detected if pinned by impurities. 

High-resolution STM experiments require high-quality 
sample surfaces, which currently restricts the application to 
BSCCO and CCOC. Although STM data have been reported 
for YBCO and 214 compounds as well, atomic resolution is 
often not achieved, and the quality of the surface layer can be 
problematic. 

STM experiments measure the current I between tip and 
sample as function of voltage V and position r. Assuming 
an energy-independent electronic density of states in the tip, 
the measured dl/dV is equivalent to the spatially resolved lo- 
cal density of states (LDOS), p(f, E), up to an r-dependent 
tunnel m atrix element which depends on the set-point condi- 
tions d37tll90l) . For cuprates, it is commonly assumed that the 
measured LDOS at low energies reflects the properties of the 
CUO2 layers, as the additional layers between crystal surface 
and topmost CUO2 plane are insulating. However, the tunnel- 
ing path through those layers may non-trivial, as discussed in 
Sec. EH 



1 . Quasiparticle interference vs. charge order 

In order to extract possibly periodic signals from LDOS 
maps p(r, E), a Fourier transformation to momentum space 
is routinely used. The resulting quantity p(k 1 E), dubbed 
FT-LDOS, can show well-defined structures for various rea- 
sons. 28 (i) If a modulation in the charge sector with wavevec- 
tor Q is present, then p(k, E) will show "Bragg" peaks at 
k = Q at all energies E. (ii) Friedel oscillations caused by 
impurities will contribute to momentum-space structures in 
p(k, E) as well. Importantly, these structures will be energy- 
dependent due to the energy-momentum dispersion of the 
single-particle excitations dl9ll) . 

The first observation of such generalized Friedel oscilla- 
tions, or "quasiparticle interference" (QPI), made by Hoff- 
man et al. dl92l) in the superconducting state of optimally 
dope d BSCCO, was subsequently confirmed and extended 
d 1 93b Il94t 1 1 95b . In interpreting the experimental data, a 
simple recipe, dubbed "octet model", was used to extract in- 
formation about the single-particle dispersion from the peak 
locations in the FT-LDOS. The extracted Fermi surface and 
superconducting gap show good agreement with results from 
ARPES dl95l) . which is remarkable, as STM is usually not a 
probe with momentum-space resolution. There are, however, 
a number of problems with the octet-model interpretation of 
QPI, and we will give a theoretical discussion in Sec. lV.EI 

The presence of QPI phenomena poses a serious problem 
in the search for charge order: As doped cuprates are intrinsi- 
cally dirty, the QPI signals are not weak and therefore are not 
easy to disentangle from modulations due to collective charge 



The analysis of the FT-LDOS is usually restricted to its power spectrum 

p(k, E)\ 2 . 
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ordering, in particular if the charge order is only short-ranged 
due to strong disorder pinning. Therefore, the experimen- 
tal search for charge order concentrated on identifying non- 
dispersive peaks in the FT-LDOS p(k,E), in particular near 
the wavevectors (tt/2, 0), (0, 7r/2) corresponding to period-4 
charge order as known from 214 cuprates. 



2. LDOS modulations 

Spatial modulations in the dl/dV signal, suggestive of 
charge order, were first detected near the vorte x cor es in 
slightly overdoped BSCCO in a applied field of 7 T (H%1) . The 
spatial pattern resembled a checkerboard, for more details see 
Sec. |III.E.3| below. 

Subsequently, different measurements of the FT-LDOS 
in zero magnetic field in optimally doped and underdoped 
BSCCO le ad to some contr overs y: While Howald et al. 
dl97l 1 198b and Fang et al. dl99l) interpreted their data as 
evidence for an underlying charge-density mod ulation (co- 
existing with QPI features), Hoffman et al. ( 192) and McEl- 
roy et al. (193) asserted that their data are consistent with 
QPI, not showing signs of a CDW. In this debate, Howald et 
al. have pointed out that the pea ks in the F T-LDOS which 
were associated with QPI in Refs. dl92Hl93l) display a much 
weaker dispersion below 15meV compared to what is ex- 
pected from the octet model (Fig. ITTb . Whether this is sim- 
ply a failure of the octet model or indeed evidence for charge 
ordering is difficult to decide, see Sec. IV.EI Moreover, Fang 
et al. ( 199) show that the height of the superconducting co- 
herence peaks display a modulation very similar to that of the 
low-energy LDOS, with a period of 4.5 a®. 




V(mV) 



k x (2;r/a ) 

FIG. 1 1 FT-LDOS data, p(k, E), of optimally BSCCO, taken at 8 K 
(reprinted with permission from Ref. (197J), copyright 2003 by the 
American Physical Society). Shown are line scans as a function of k x 
along the (1,0) direction, and as a function of energy (color scale). 
Top: LDOS (dl/dV). Bottom: LDOS integrated up to the given 
energy (I/V). The peak at 27r(0.25, 0) displays little dispersion. 



Unambiguous evidence for charge oder (albeit with a short 
correlation length) came from LDOS measurements in the 
pseudogap regime of BSCCO by Vershinin et al. d200h . An 
underdoped BSCCO crystal with a T c of 80 K, measured at 
100 K, displayed non-dispersive peaks in the FT-LDOS be- 
low 40meV, corresponding to a modulation period of 4.7 ao, 
Fig. Q~2] The observe d FT-LDO S signal was shown to be in- 
compatible with QPI (200; 201), essentially because the octet 
wavevectors qi and 175 did not disperse at all. (The measured 
data below T c displayed dispersing peaks, compatible with 
Ref. dl92l) .) A plausible interpretation is that QPI signatures 
become progressively weaker with increasing temperature, as 
quasiparticles are ill-defined in the pseudogap regime over a 
large fraction of the Brillouin zone. As a result, it becomes 
easier to identify signatures of CDW formation. Whether, in 
addition, charge order is enhanced in the pseudogap regime is 
unclear. 

More recent STM measurements of various underdoped 
BSCCO samples have established that non-dispersive FT- 
LDOS peaks also occur in the superconducting state, namely 
at elevated energies above 30 ... 40 meV. While those peaks 
were originally associated with pseudogap patches of the in- 
homogeneous sample dl941) . they have been recently argued 
to be generic features of the approach to the Mott insulator 
(|202|) (see also next subsection). A connection between strong 
gap inhomogeneities and the presence of non -dispersive FT- 
LDOS peaks has been pointed out in Refs. d203c 12041) : In 
both the superconducting and the pseudogap regimes, non- 
dispersive peaks occurred predominantly in strongly inhomo- 
geneous samples or sample regions. This appears plausible 
under the assumption that pinning due to disorder plays a 
dominant role in inducing the static modulation. 

A charge-order signal much stronger than that in BSCCO 
was identified in the low-temperature LDOS of unde rdoped 
samples of CCOC with dopings x = 0.08 ... 0.12 (fl35l) . Here, 
the LDOS at all energies below the pseudogap energy of 
100 meV displayed a clearly visible checkerboard modulation 
with spatial period 4. 29 ARPES experiments on similar sam- 
ples of CCOC detected nearly nested, but incoherent, antin- 
odal Fermi surface pieces with almost doping-independent 
nesting wavevectors close to (7r/2,0), (0, 7t/2). The au- 
thors proposed that these antinodal regions are responsible 
for c harge ordering in the spirit of a weak-coupling scenario 
d206h . 

Very recently, non-dispersive FT-LDOS peaks at energies 
below 30meV were reported in single-layer BSCCO-2201 
(12071) . The modulations displayed a doping-dependent spa- 
tial period between 4.5 ao (underdoped) and 6.2 ao (optimally 
doped). This tendency is opposite from the one of the CDW 
wavevector in 214 compounds (where the modulation period 
becomes smaller with increasing doping), and consequently 
a weak-coupling scenario has been suggested as the origin of 



Early STM data taken on CCOC at a bias voltage of 200 meV have been 
interpreted in terms of nanoscale inhomogeneities <205f) - The more recent 
high-quality data of Ref. ( 37) show that those inhomogeneities correspond 
to uni-directional bond-centered modulation patterns. 
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FIG. 12 Fourier analysis of DOS modulations measured in slightly underdoped BSCCO (T c = 80 K) in the pseudogap regime at 100 K 
(from Ref. (200), Science 303, 1995 (2004), reprinted with permission from AAAS). (A) Fast Fourier transform (FFT) of an unprocessed 
conductance map acquired over a 380 A by 380 A field of view at 15 mV. (B) The FFT has peaks corresponding to atomic sites (colored black 
and labeled A), primary (at 27r/6.8ao) and secondary peaks corresponding to the b-axis supermodulation (colored cyan and labeled S), and 
peaks at « 2-7r/4.7ao along the (H, 0) and (0, H) directions (colored red and labeled Q). (C) The energy evolution of the peaks in (B), scaled 
by their respective magnitudes at 41 mV. (D) Two-pixel-averaged FFT profiles taken along the dashed line in (B) for the DOS measurement at 
15 mV shown in (A) and measurements acquired simultaneously at mV and —15 mV. The positions of key peaks are shown by dashed lines 
and labeled according to their location in (B). 



the modulations (12071) . 

Field-induced LDOS signals near vo rtex cores have been 
investigated in more detail in BSCCO J201 12091: l2loh and 
also observed in YBCO d21 ll) . and will be described in 
Sec.lIILEJbelow. 

In the experiments described so far, clear-cut evidence for 
(locally) broken rotation symmetry was lacking (although the 
data of both Hoffman et al. d 196b and Howald et al. dl97l) in- 
dicate a weak local breaking of C4 symmetry). Superficially, 
this might be more consistent with checkerboard than stripe 
order. However, as discussed in Sec. III.G1 pinning of stripes 
by impurities can result in checkerboard (instead of stripe) 
patterns, in particular if the clean system is on the disordered 
side of a stripe-ordering transition (65; 66). 



3. Tunneling asymmetry 

A critique which has been voiced against the analysis of 
modulations in dl/dV, as described in the last subsection, 
is related to the possibly position-dependent tunneling matrix 
element in the STM experiment. The standard measurement 
protocol is to adjust the tip height at each position such that a 



constant current flows at some (high) set-point voltage. In the 
presence of significant charge inhomogeneities, the tip height 
will then be modulated as a function of r as well, rendering 
dl/dV inequivalent to the LDOS. Moreover, this effect will 
depend on the set-point conditions. 

To separate physical modulations from set-point effects, it 
was proposed to study the LDOS ratios 

Z<r,E) = -02-, R{r,E)= [»' L ! PM . (11) 



p{r,-E) : 



J-e dup{r 7 u) 



The physical content of these ratios, which measure spectral 
particle-hole asymmetry, is non-trivial. For one-band models 
of weakly doped Mott insulators, both Z a nd R hav e been 
argued to be proportional to the hole density d212tl213l) . 

Maps of R(r, E), taken in underdoped BSCCO and CCOC, 
have been analyzed by Kohsaka et al. (37). The tunneling 
asymmetry at 150mV was found to be strongly inhomoge- 
neous, with a spatial pattern interpreted as "electronic cluster 
glass". A number of properties of the observed R(f, E) are in- 
deed remarkable: (i) The modulations are centered on the Cu- 
O-Cu bonds and are strongest on the O (instead of Cu) sites, 
with a contrast of up to a factor of two. (ii) The modulation 
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FIG. 13 Valence-bond glass as seen in STM (from Ref. {37]), Science 315, 1380 (2007), reprinted with permission from AAAS). A) and D): 
R maps of Na-CCOC and Dy-Bi2212, respectively (taken at 150 mV). The blue boxes in A), D) indicate areas shown in panels B), C) and 
E), F). B and E: Higher-resolution R map within equivalent domains from Na-CCOC and Dy- Bi2212, respectively. The locations of the Cu 
atoms are shown as black crosses. C) and F) Constant-current topographic images simultaneously taken with panels B) and E), respectively. 
The markers show atomic locations, used also in B) and E). 



pattern locally breaks to C4 symmetry down to C2, i.e., uni- 
directional domains, suggestive of stripe segments, are clearly 
visible, (iii) The modulations in the R map of BSCCO and 
CCOC are essentially indistinguishable (although differences 
exist in the energy-resolved LDOS spectra). In both cases, the 
spatial correlation length of the modulation pattern is of or- 
der 10 lattice spacings. (iv) In both cases, the order co-exists 
with well-established superconductivity, with a T c of 45 K and 
21 K for BSCCO and CCOC, respectively. 

Subsequently, maps of Z(r, E) w ere used to analyze QPI 
in slightly underdoped CCOC (190) and strongly underdoped 
BSCCO J2021) . In CCOC, the spatial period-4 modulations 
are so strong that QPI features are hard to detect in the LDOS 
(Il35l) . In contrast, QPI peaks were observed in the LDOS 
ratio Z(k, E) below 20meV dl90l The fact that Z is more 
sensitive to QPI features than the LDOS itself was ratio- 
nalized arguing that charge order causes modulations in the 
LDOS which are approximately in-phase between positive 
and negative energies, whereas QPI modulations occur ap- 
proximately anti-phase. In BSCCO, dispersing QPI peaks in 
Z(k, E) were observed at low energies as well d202l) . How- 
ever, above a doping-dependent crossover energy the Z(k, E) 
spectra change: Some QPI peaks disappear, while others 
cease to disperse. Interpreting the dispersive peaks within the 
octet model suggests that the QPI signal exists for quasipar- 
ticles inside the antiferromagnetic Brillouin zone only (12021) . 
The non-dispersive features at elevated energies correspond to 
a modulation period of approximately 4 ao, with little doping 
dependence. Interestingly, the modulations in Z(r, E) appear 
strongest (as function of energy) at the local pseudogap energy 
Ai which is of order lOOmeV. 



Although a comprehensive understanding of the asymme- 
try maps Z{r, E) and i?(r, E) is difficult, as it necessarily 
involves Mott physics, the described results strongly suggest 
that a tendency toward bond-centered stripe order is present 
at the surface of both underdoped BSCCO and CCOC. This 
stripe-like order, primarily visible at elevated energies, co- 
exists with well-defined low-energy quasiparticles and bulk 
superconductivity. 30 In both BSCCO and CCOC, the charge 
order exists without long-range magnetic order, although spin- 
glass-like magnetism has been reported in CCOC ( 12141) . The 
result in Fig. [13] showing that the main modulation is on oxy- 
gen, suggests that a quantitative description requires a three- 
band instead of a one-band model for the G1O2 planes. 



D. Other probes 

Besides neutron scattering and STM, a broad variety of 
other experimental techniques have been employed to detect 
and investigate ordering phenomena accompanied by lattice 
symmetry breaking. Here we shall give a quick overview, 
without pretense of completeness. 



It should be noted that a precise relation between the observed modulations 
in the tunneling asymmetry and those in the LDOS has not been established 
experimentally. 
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1. Transport 

Transport measurements have been used to search for ten- 
dencies towards order. Two possible signatures are obvious: 
(i) If the ordered phase is less conducting than the disordered 
phase, then the resistivity will show an upturn below the or- 
dering temperature, (ii) If the ordered phase globally breaks 
rotation symmetry, then the resistivity will develop an in-plane 
anisotropy. 

With the exception of the small-doping insulating regime, 
i.e., x < 0.055 in LSCO, stripe phases appear generally con- 
ducting, i.e., are not of Wigner crystal type. Temperature- 
dependent resistivity measurements often show a small upturn 
feature upon cooling below the stripe ordering temperature 
T c h. Although dp/dT < at low T in stripe-ordered cuprates 
(above T c ), an exponential rise of the resistivity signaling in- 
sulating behavior is never seen. 

Optical conductivity measurements in LBCO-1/8 d215l) find 
a residual Drude peak even in the charge-ordered regime at 
low T, which has been interpreted as a "nodal metal" state, 
i.e. ungapped nodal quasiparticles coexist with stripe order. 
This finding is consistent with the absence of a gap in ear- 
lier optical-conductivity data taken on LNSCO-1/8 d216h . The 
presence of nodal quasiparticles in stripe phases is also com- 
patib l e wit h resu lts fromphotoemission and STM experiments 
d 1 3€t 1 1901: 12171). Note that, on the the ory s i de, a rguments 
against d218t |219|) and in favor 4H HH El iHl of nodal 
quasiparticles in stripe phases have been put forward. For 
small stripe amplitude, the survival of nodal quasiparticles im- 
mediately follows from the fact that the ordering wavevector 
Q does not connect the nodal points d47t 12221 : |223|) . while for 
larger stripe amplitude a substantial <i-wave com ponen t of the 
charge order can protect nodal quasiparticles (136c 12241) . 

The optical conductivity data in LBCO-1/8 (12151) show a 
rapid loss of spectral weight below 40 meV which occurs be- 
low about 60Jf ~ T c h, consistent with the development of 
an anisotropic gap. However, the existence of strong far- 
infrared peaks (between 20 and 100 cm -1 ) at low tempera- 
tures, present in earlier data on LSCO and LNSCO (12 161 : 12251) 
and interpreted in terms of stripe pinning d226l) . could not be 
confirmed. This disagreement has been tentatively attributed 
to surface problems (215). 

A detailed resistivity study of LSCO at finely spaced dop- 
ing s J227I) identified specific doping levels, defined by peaks 
of the resistivity at fixed temperatures of 50 and 100 K as func- 
tion of doping. These peaks, albeit being weak, have been in- 
terpreted as magic doping fractions 1/8, 1/16, 3/16, 3/32, 5/32, 
and suggested to be a signature of a hierarchy of checker- 
board charge-ordered states. While this idea is interesting, 
other measurements in 214 materials are more consistent with 
stripe instead of checkerboard order. The presence of special 
commensurate doping levels with enhanced ordering is also 
expected for lattice pinning of stripes, although the set of spe- 
cial doping fractions may be different. 

Let me now turn to more direct ordering signatures. States 
with broken in-plane 90-degree rotation symmetry will dis- 
play anisotropic transport properties, e.g., a d.c. resistivity 
tensor with p xx ^ p yy . This applies both to nematic and to 



stripe phases, but, as discussed in Sec. III. Gl requires a single- 
domain sample. Hence, de-twinned YBCO is again the prime 
candidate. The magnitude of the anisotropy will depend on 
details of the order and the electronic scattering mechanisms: 
A metallic nematic state has a full Fermi surface (the same 
applies to most small-amplitude stripe states), hence velocity 
and scattering anisotropics become important. 

In-plane resistivity anisotropics in de-twinned YBCO with 
dopings S — 0.35 ... 1 have been investigated by Ando et al. 
(11881) . While an anisotropy is present for all temperatures 
and dopings due to the presence of the CuO chains in this 
material, the anisotropy decreases with cooling below 200 K 
for 5 > 0.6, but increases for S < 0.6. This effect is par- 
ticularly significant for 6 — 0.35 and 0.45, where the ratio 
Pal Pb increases from 1.3 at room temperature to 2 ... 2.5 at 
low T. This result is consistent with the onset of nematic or- 
der inferred from neutron scattering in YBCO-6.45 (91), see 
Sec. IIII.B~4l Note that the presence of a sharp onset tempera- 
ture cannot be expected, because rotation symmetry is broken 
by the orthorhombic distortion of the crystal from the outset. 

In stripe-ordered 214 compounds, global resistivity 
anisotropics cannot be observed, due to the plane-to-plane al- 
ternation in the stripe direction in the LTT phase. Noda et 
al. d228l) instead measured magnetotransport in stripe-ordered 
LNSCO at various doping levels. For doping x < 1/8, the 
data show a distinct drop of the Hall coefficient at a tempera- 
ture of order 80 K where stripe charge order is believed to set 
in. This has been interpreted as evidence for one-dimensional 
charge transport deep in the stripe-ordered phase of LNSCO 
for x < 1/8, with a crossover to more two-dimensional trans- 
port for a; > 1/8. However, subsequent theoretical work indi- 
cated that this picture is too simple: A nearly vanishing Hall 
coefficient can be obtained in a quarter-filled stripe phase in- 
dependent of the total doping level, not as a result of one- 
dimensional transport, but of an appro ximate pa rticle-hole 
symmetry at this particular stripe filling (229; 230). 31 A rapid 
drop in the Hall coefficient upon cooling has also been ob- 
served in LBCO-0. 1 1 d233l) . with the Hall coefficient becom- 
ing negative at lowest temperatures. A simplistic interpreta- 
tion would be in terms of electron pockets arising from band 
backfolding (234), but correlation effects may again change 
the picture d230l) . 

Measurements of the Nernst effect in both LNSCO and 
LESCO showed a distinct low-temperature enhancement of 
the Nernst signal for dopings between 0.12 and 0.20, with the 
signal onset occurring at a t emperature T„ which tracks the 
pseudogap temperature d235l) . (Interestingly, the T v of LSCO- 
0.12 and LESCO-0.12 are essentially equal, around 140 K.) 
As T v roughly follows the onset of stripe order, the Nernst 
signal was interpreted as evidence for a stripe-induced Fermi- 



31 An approximate particle-hole symmetry was also inferred from ther- 
mopower measurements in the antiferromagnetic low-doping state of 
YBCO 1 231). Its relation to stripe physics is unclear. A distinctly doping- 
dependent behavior of the thermopower was measured in stripe-ordered 
LNSCO, which was interpreted as evidence for a Fermi-surface change 
near optimal doping d232l) . 
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surface reconstruction. However, the present data appear as 
well consistent with a pseudogap-induced Nernst signal of 
non-stripe origin. Measurements at lower dopings might help 
to disentangle these possibilities. 



2. Photoemission 

Angle-resolved photoemission (ARPES) is the method of 
choice to obtain momentum- and energy-resolved informa- 
tion on the single-electron spectral function. 32 Compared to 
neutron scattering, ARPES is only sensitive to a surface layer 
of the sample with a thickness of order 5 ... 20 A, depend- 
ing on the photon energy. Quality and resolution of ARPES 
data have improved tremendously over the past decades, al- 
though the current ARPES energy resolution, being typically 
2. . . lOmeV, is not as good as that of state-of-the-art neutron 
scattering. An experimental problem is that ARPES requires 
samples with high-quality surfaces; those are routinely avail- 
able for BSCCO and CCOC. Extensive ARPES studies have 
also been performed on a number of 214 cuprates, with by 
now comparable data quality. In contrast, YBCO suffers from 
a charge imbalance at the surface, rendering the interpretation 
of corresponding ARPES difficult. 

I start with a description of a few general features of 
ARPES spectra on underdoped cuprates, primar ily obtained of 
BSCCO crystals. A common observation dH ^ [2061: l23~~t 
|237|) is that low-energy electronic states along the diagonal 
(nodal) direction in momentum space appear to be rather well 
defined, i.e. produce sharp peaks in energy-distribution curves 
(EDC). In contrast, states near (±7r, 0) and (0, ±7r) (antinodal 
points) are broad and appear gapped even above T c . The dis- 
parate behavior of the two regions in momentum space, also 
dubbed "nodal-antinodal dichotomy", is also reflected in the 
temperature dependence of the linewidths: Upon lowering T, 
nodals become sharper, while antinodals tend to be broader. 

Below T c , experiments find a gap consistent with d-wave 
symmetry, sometimes mixed with higher harmonics such that 
the gap near the nodes is smaller compared to a pure <i-wave 
shape. Moving into the pseudogap regime above T c , antin- 
odals remain gapped, whereas near-nodal states are essentially 
gapless, leading to apparent segments of Fermi surface. These 
observation lead to the concept of "Fermi arcs", with a tem- 
perature and doping dependent arc length d2~"~h . From theo- 
retical studies, it has been proposed that the arcs may in fact 
be Fermi pockets centered along the nodal direction, with the 
outer part of the pockets being near ly invisible to ARPES due 
to matrix-element effects (1238 1239). The distinct temperature 
and doping dependence of the antinodal and near-nodal gaps 
has also prompted proposals of a two-gap scenario, with the 
gaps in the two regions in momentum space being caused by 
different underlying physics J240h . 

The absence of well-defined antinodal quasiparticle has 



sometimes been interpreted as (indirect) evidence for order- 
ing tendencies (e.g. of CDW type) primarily carried by antin- 
odals d206), which could be consistent with an RPA picture of 
CDW formation. However, opposite experimental views have 
also been put forward d24ll) . 

I now turn to the question how ARPES spectra are directly 
affected by symmetry-breaking order: Broken rotation sym- 
metry should be visible directly in momentum-resolved data 
of properly aligned crystals, and broken translation symmetry 
should lead to Bragg gaps and multiple bands via band back- 
folding. 

To my knowledge, unambiguous evidence for sponta- 
neously broken rotation symmetry has not been detected in 
cuprate ARPES to date. As discussed above, the simultane- 
ous presence of domains with different preferred directions - 
either in a single plane or in adjacent planes - makes rotation 
symmetry breaking difficult to observe (unless the experiment 
would be sensitive to the topmost CUO2 layer only, which in 
addition had to be single-domain). The exception would again 
be YBCO, with complications arising from the presence of 
CuO chains; however, high-quality photoemission data for un- 
derdoped YBCO, showing the full temperature dependence of 
the spectrum, is not available to my knowledge. 

With broken translation symmetry (or other stripe signa- 
tures) in ARPES, the situation is more inv olve d. An early 
experiment on stripe-ordered LNSCO-0.12 d242l) at 20 K has 
found distinct signatures near the Fermi energy, namely a 
cross-shaped intensity distribution in the Brillouin zone and 
no evidence for nodal quasiparticles. These features were at- 
tributed to stripes with strong charge modulation; in particular, 
the straight horizontal and vertical Fermi-surface pieces near 
the antinodal points were interpreted as evidence for promi- 
nent one-dimensional behavior. 33 This view was subsequently 
supported by model calculations which started from weakly 
coupled chains or ladders and assumed a strong ld-like mod- 
ulation (Sec. IV.D~ l. Howeve r, the results of ( 12421) remained 
controversial ( 244 ; 12451 ; |246|) . and later experiments with im- 
proved samples and energy resolution could neither verify the 
strong cross-shaped intensity pattern nor the absence of nodal 
quasiparticles d245[); instead, Fermi arcs were found at 15 K 
above T c = 7K J246I) . However, clear-cut stripe signatures 
were not detected. 

For stripe-ordered LBCO-1/8, an ARPES study (fl36l) at 
16K, much below the charge-ordering temperature, also 
failed to detect signatures of stripe-like superstructures, but 
instead reported an anisotropic gap along the Fermi surface, 
with angle dependence being consistent with a rf-wa ve form 
factor. A recent higher-resolution experiment d247l) mainly 
confirmed this data, but also indicated that the gap slope 
changes rather abruptly along the Fermi surface, suggestive 
of a two-component gap. While the large antinodal gap may 
be interpreted in terms of pseudogap physics, the origin of 



32 Complications arising from non-trivial final states and from a dispersion 
perpendicular to the C11O2 planes shall be ignored here. 



The conclusions in Ref. (242) were mainly drawn from ARPES data which 
had been integrated over an energy window of 0.5 eV around the Fermi 
level. A theoretical consideration on energy-integration effects in ARPES 
in the presence of fluctuating charge order is in Ref. 12431) . 
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the smaller near-nodal gap is unclear at present, with fluctuat- 
ing d-wave pairing (1431) or d-wave stripe order ( 36) being two 
candidates. 

Other recent ARPES measurements are worth mentioning: 
Data obtained on stripe-ordered LESCO-1/8 (12171) displayed 
weak signatures of both rotation symmetry breaking and 
band backfolding, which may be consistent with the expected 
period-4 charge order. However, not all features are consis- 
tent with model calculations (Sec. IV.DI ). For stripe-ordered 
LNSCO-1/8, Ref. d246h verified the existence of Fermi arcs 
and, in addition, observed a second Fermi surface crossing 
near the Brillouin zone diagonal, suggestive of band back- 
folding leading to Fermi surface pocket. The location of this 
second branch appears consistent with (it, ir) order, i.e., a unit 
cell doubling. The branch is visible up to 1 10 K, i.e., far above 
the charge ordering and LTT transitions, and a similar branch 
is seen in LSCO-0.12 (but not at other LSCO doping levels). 
We note that hints for the existence Fermi surface pockets 
have also been detected in ARPES in the pseudogap regime 
of underdoped BSCCO (12481) . however, with the pocket not 
being centered around (tt/2,tt/2). At present, the explana- 
tion of these observations is open. 

In summary, a satisfactory experimental verification of the 
theoretical expectation for electronic spectra of stripe phases 
is still missing. Clearly, more experiments are called for, 
with candidate materials being the stripe-ordered 214 materi- 
als, studied systematically as function of doping, temperature, 
and field, and underdoped YBCO, where signatures of rota- 
tion symmetry breaking have already been detected in neutron 
scattering. 



3. fxSR, NMR, and NQR 

Techniques sensitive to local properties have been very use- 
ful in characterizing ordering tendencies in cuprates. /iSR has 
been used to detect magnetic order via static (or slowly fluc- 
tuating) hyperfine fields; magnetic ordering as determined by 
/iSR typically gives lower ordering temperatures T sp as com- 
pared to elastic neutron scattering (e.g. Fig. @), due to the 
different frequency windows of the two techniques. 34 /iSR 
also allows to determine the magnitude of the magnetic mo- 
ment relative to some known reference (e.g. the undoped 
parent compound) as well as the volume fraction of the or- 
der. For striped 214 compounds, the moment size is typically 
found to be half of that in the und oped compound, but with re- 
duced volume fraction d 1 08k 1 1 09b . /iSR has been used to map 
out the magnetic phase diagram of LESCO, where spin stripe 
ord er rep laces superconductivity over a large doping rang e 
dUdla, as well as of YBCO, LSCO, and LBCO £]3t.\2A%. 

NMR and NQR provide related information and have been 
used extensively in the context of stripes. The broadening of 
NMR lines at low temperature in LSCO d250t; l25ll) provides 



evidence for local spatial inhomogeneities. NQR measure- 
ments have detected a so-called "wipeout effect", i.e. the 
gradual loss of the Cu NQR signal below some tempera- 
ture. While this was origi nally interpreted as direct mea- 
sure of charge stripe order (1252b , subsequent work has es- 
tablished that the wipeout is due to a slowing down of spin 
fluc tuations which accompany the tendency toward stripe or- 
der I2l^l2l^l255l) . 

Subsequent NMR/NQR studies p rovided details about the 
spatial distribution of doped holes (fl89t [256t l257h in LSCO 
and YBCO, giving evidence for spatial inhomogeneitities at 
low doping. In underdoped YBCO, with oxygen content be- 
low 6.4 and located close to the transition to the insulator, 
/iSR has detected static magnetic order o f glas sy character 
which co-exists with superconductivity (19 It 12581) . For YBCO 
at larger doping, it has been argued from NQR measurements 
that charge inhomogeneities arise only from the presence of 
local oxygen vacancies (HI); this would be consistent an 
SDW and/or CDW critical point being located at very small 
doping. A particular remarkable finding is that of two inequiv- 
alent planar O sites, with only one type of Cu site, in YBCO 
over essentially the entire doping range d 1 891) . This interest- 
ing result may be related to the tendency toward electronic 
nematic order. Unfortunately, a detailed experimental study 
of the temperature dependence of oxygen NMR/NQR has not 
been performed to date. Together with an in-depth theoreti- 
cal analysis, this would be highly revealing about nematic or- 
dering. NMR measurements have also explo ited the field de- 
pendence of magnetism in YBCO d260H26l1) . see Sec. IIII.E31 
below. 

NMR has been as well applied to LESCO dHHHl) which 
is known to display static stripes, with charge order setting in 
around 75K at x = 1/8. The results show a distinct change 
in the O NMR spectra below 80 K, and have been interpreted 
(12631) in terms of a correlation between local hole doping and 
domain walls in the spin modulation, as expected from stripe 
phases. Also, the static stripe order and/or the LTT distortion 
characteristic of LESCO apparently act to suppress the very 
slow spin fluctuation s, which a re present in LSCO, in favor of 
ordered magnetism d262tl263l) . 



E. Relation to superconductivity 

The experiments described so far have established the ex- 
istence of stripe-like order in certain underdoped compounds, 
and the data suggest that the phenomenon is in fact common to 
a variety of cuprate families. Then, a central question is about 
the relation between stripes and superconductivity. A num- 
ber of relevant experiments have been performed that will be 
summarized in the following. A detailed theoretical discus- 
sion will be given in Sec. IV.FI 



1 . T c suppression and 1/8 anomaly 



Further differences in T sp may originate from the different behavior of 
polycrystals and single crystals, compare e.g. Refs. 1 75) and f76l) . 



Do stripes promote or inhibit superconductivity? First in- 
formation can be obtained by looking at the transition temper- 
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atures of both ordering phenomena. The bulk superconduct- 
ing T c as function of doping follows the well-known parabolic 
dome shape - this applies in particular to compounds with 
little stripe signatures, e.g. multilayer BSCCO. In contrast, 
in the single-layer 214 materials which tend to stripe forma- 
tion, the situation is different: The superconducting materials 
LSCO and LBCO display a pronounced T c suppression at ap- 
proximately 1/8 doping, first observed in LBCO (71) (which 
was later found to display static stripe order). This phe- 
nomenon has been termed "1/8 anomaly". If the LTT distor- 
tion is stabilized by co-doping as in LNSCO and LESCO, then 
superconductivity is replaced by non-superconducting stripe 
order over a large part of the doping range. Notably, also bi- 
layer YBCO displays features which may be of similar origin 
as the 1/8 anomaly in LSCO: The doping dependence of T c 
in YBCO shows a shoulder near oxygen content 6.4, a nd th e 
penetration depth is found to have a distinct maximum (264). 

As stripe order is known to be particularly strong near 1/8 
doping, the 1/8 anomaly is commonly attributed to stripes 
competing with superconductivity. In line with this inter- 
pretation, most cuprate materials display either quasi-long- 
range stripe order (like LESCO and LNSCO for x < 1/8) 
or well-developed bulk superconductivity. Exceptions, with 
static stripes apparently co-existing with superconductivity at 
lowest T, are LBCO near x — 1/8 as well as LNSCO and 
LESCO at larger hole doping, x > 1/8. 35 In LSCO with 
x < 1/8 superconductivity co-exists with quasi-static incom- 
mensurate spin order, but here no direct evidence of charge 
order has been reported. In superconducting LSCO, spin or- 
der can be induced for larger x by applying a magnetic field, 
see Sec. IIII.E.3I In underdoped BSCCO and CCOC, the 
signatures of charge order seen by STM co-exist with well- 
developed superconductivity. 

For LBCO with 0.11 < x < 0.15, the stripe ordering tem- 
perature is above the superconducting T c . The most extensive 
data set is available at doping x = 1/8, displa ying a remark- 
able hierarchy of temperature scales (F73i : |266|) which demon- 
strates that physics is far more complicated than a simple com- 
petition of two ordering phenomena would suggest. The struc- 
tural transition to the LTT phase is located at 55 K, immedi- 
ately followed by the charge ordering transition around 54 K. 
Spin order sets in around 42 K, accompanied by a large drop in 
the in-plane resistivity p a i,. However, p a b only vanishes below 
16 K, whereas p c vanishes below 10K. Finally, a bulk Meiss- 
ner effect is seen below 4 K. This set of highly non-trivial data 
is puzzling: In addition to competing orders, it likely involves 
dimensional crossover and disorder effects. The temperature 
regime above 10 K has been interpreted in terms of fluctuat- 
ing 2d superconductivity, and we will come back to this in 
Sec. IV.F.2I A detailed experimental discussion is in Ref. ©. 

The competition between superconductivity and stripes, 
which are stabilized by the LTT distortion of some of the 214 
compounds, is visible in a number of other effects: In LESCO, 



where stripes exist for x < 0.2, the oxygen isotope effect on 
the superconducting T c was studied in in Ref. (68). For a non- 
stripe-order x — 0.24 sample essentially no isotope effect on 
T c was seen, whereas T c decreased significantly by 1 .5 K in an 
x = 0.16 sample upon substituting 16 by 18 0. A plausible in- 
terpretation is that the isotope substitution primarily strength- 
ens stripe order, thereby suppressing superconductivity. This 
would then also be a nice confirmation of stripe-lattice cou- 
pling, although a detailed modeling is not available. 

Related are t he pressu r e studies of 214 compounds in 
Refs. dH El |267t [26^ Hydrostatic pressure 

applied to stripe-ordered LNSCO- 0.12 l eads to a large in- 
crease of the superconducting T c (268; l270h . simila r dat a 

"" 2671) and LESCO (HH). 



have been obtained for LBCO (1 11 61: 



Uni-axial in-plane pressure increases the superconducting T c 
in both non-stripe-ordered LSCO and stripe-ordered LESCO 
and LNSCO. However, while the pressure directions [100] 
and [110] are equivalent for LSCO, there is a significant 
anisotropy for stripe-ordered LESCO-0.16 d58l) . with the T c 
increase occurring for pressure along the [110] direction. 
From these experiments, it is believed that pressure reduces 
the tilt angle of the CuOg octahedra, i.e. tends to stabilize 
the orthorhombic symmetry of the lattice, thereby suppress- 
ing stripe order and hence enhancing superconductivity. The 
maxi mum T c is typically reached at high pressures in the HTT 
phase (270) when the CUO2 planes are flat. Interestingly, the 
pressure-induced increase in T c is r ather slow for x — 1/8 as 
compared to x = 0.12 or 0.13 d267l) . and signatures of stripe 
ord er hav e been found to survive in the HTT phase of LBCO- 
1/8 (Gil). 



2. Temperature dependence of stripe signatures 

Insight into the interplay of ordering phenomena can be 
gained by studying their temperature dependence. Distinct 
behavior is found depending on whether stripes or supercon- 
ductivity dominate. 

I begin with the materials showing well-developed stripe 
order according to scattering probes. For LBCO with 0.11 < 
x < 0.15, LNSCO, and LESCO charge order sets in below the 
LTT transition followed by quasi-static spin order at lower T, 
and superconductivity with a very small (or zero) T c . While 
T c h < Tltt in LBCO and LNSCO, the two transition are 
well separated with T" c h < Tltt in LESCO. All compounds 
have T S p < T c h, i.e., charge order exists without spin order. 36 
No appreciable changes of charge and spin order have been 
reported upon crossing T c . 

The materials LSCO-0.12 (HI) and La 2 Cu0 4 .i 2 dsl), dis- 
playing SDW order (without reported signatures of charge 
order), share an interesting feature: Within the experimen- 
tal accuracy, the superconducting T c equals T sp . This may 
be related to a property of LBCO- 1/8 noted above: The spin- 



For LNSCO and LESCO, it has been debated whether the superconductiv- 
ity co-existing with stripes is of true bulk nature I; 263) . 



36 Note that T sp depends on the employed probe, suggesting a wide regime 
of slow fluctuations. 
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ordering transition at T sp is accompanied by a large drop in the 
in-plane resistivity (but bulk superconductivity is not reached 
here). Whether these findings are coincidental or rather sig- 
nify a positive correlation between magnetic order and pair 
formation is unclear at present. Note, howeve r tha t a mag- 
netic field separates T sp from T c in La2Cu04.i2 d27ll) and T sp 
from the resistivity drop in LBCO-1/8 (266). Phase separa- 
tion has been under discussion for La2Cu04+5 dl05l) : inter- 
estingly, a theoretical scenario based on phase coexistence of 
competing SDW and superconducting orders h as be en shown 
to allow for T c = T sp for a range of parameters (106). 

I now turn to the temperature dependence of incommen- 
surate spin excitations. For stripe-ordered LBCO-1/8, the 
evolution of the low-T hour-glass spectrum was studied in 
Ref. d272l) . The low-energy incommensurate signal is well 
visible at 65 K above T c h w 55 K, but deteriorates at 300 K. 
However, the upper excitation branch above 50 meV changes 
little from 5 K to 300 K. This was interpreted as evidence for 
a stripe liquid in Ref. d272l) . In LSCO, the temperature depen- 
dence of the excitation spectrum is similar: While excitations 
at high energies display only a weak temperature dependence, 
the low-energy excitations become much more structured at 
low T. While in optimally doped LSCO this crossover is 
tied to T c dl72l) . it occurs at a higher temperature in under- 
doped LSCO-0.085, which has been attributed to the pseu- 
dogap ( 164). Remarkably, incommensurate excitations below 
T c are even seen in overdoped LSCO dl74l) . but disappear for 
doping x = 0.30 outside the superconducting dome. 

In optimally and slightly underdoped YBCO, a clear hour- 
glass dispersion with a sharp resonance at (tt, it) is observed 
at low T. In YBCO-6.95, the lower branch disappears com- 
pletely above T c , with damping being rather strong (I170I) . In 
YBCO-6.6 with T c = 61 K, the low-T hour glass is replaced 
by a "Y"-shaped dispersio n at 70 K, and the resonance in- 
tensity is strongly reduced (fl69l) . Turning to strongly under- 
doped YBCO, a "Y"-shaped dispersion is seen in YBCO-6.45 
(T c = 35 K) both above and below T c , with the low-energy 
intensity being strongly enhanced at low T, however, this in- 
crease does not appear tied to T c d9ll) . 

It is worth noting an interesting aspect of the tempera- 
ture dependence of spin fluctuations: While in optimally and 
overdoped compounds, cooling and in particular the onset of 
superconductivity usually reduces the magnetic low-energy 
spectral weight (i.e. often a clear spin gap opens below T c ), 
this is not the case in strongly underdoped materials, where 
instead magnetism appears enhanced. 

Finally, I briefly discuss the CDW signatures as detected 
with STM. Here, only a few studies a t elevated temperatures 
are available and restricted to BSCCO d200l: l204h . The data in- 
dicate robust period-4 CDW features in the pseudogap regime, 
likely pinned by sample inhomogeneities. The comparison 
with data obtained at low T is complicated by the signatures 
of quasiparticle interference at low T which are weak or ab- 
sent at elevated T: At low T, it is difficult to assign peaks in 
the FT-LDOS uniquely to either quasiparticle interference or 
charge order (while this ambiguity does not exist above T c ), 
see Sec. IIII. CI Whether this implies that CDW tendencies are 
enhanced in the pseudogap regime as compared to the low-T 



superconducting phase is not known. 

While a one-line summary is difficult, a common feature 
seems that in underdoped compounds the ordering tendencies 
develop below the pseudogap temperature (although long- 
range order may set in only at very low T). Moving towards 
optimal doping, ordering tendencies weaken, and the influ- 
ence of the superconducting T c becomes more significant, 
with strongly damped behavior above T c . 



3. Magnetic field tuning of order 

The application of a c-axis magnetic field to superconduct- 
ing cuprates has been found to enhance the tendency towards 
spin and charge order in a variety of cuprates. These experi- 
ments provide strong evidence for a competition of stripes and 
superconductivity, as will be discussed below. The concrete 
modeling of this competition is described in Sec. lV.F.ll 

A series of experiments is concerned with enhanced mag- 
netism in an applied field. Particularly striking are low- 
temperature neutron scattering experiments on 214 com- 
pounds. Where well-developed superconductivity and quasi- 
static incommensurate spin order coexist, the latter is found 
to be significantly enhanced by a field of order 1 0T. T his ap- 
plies to LSCO with x < 1/8 (El [ml Hi! and to 
La 2 Cu0 4 +5 dJzHEzl. 

For LSCO-0.10, the ordered moment 
in creas es from 0.03 /is at zero field to 0.15 /ib per Cu at 14.5 
T (1273b . Inelastic measurements show that the spectral-weight 
transfer is restricted to energies below 1 .5 meV, albeit with a 
non-trivial temperature dependence dl7ll) . In superconduct- 
ing compounds without quasi-static order, but located close to 
a magnetic instability, an applied field can reduce the gap to 
spin excitations and eventually drive the system into a mag- 
netically ordered sta te. This is again nicely seen in LSCO 

dnl l274tj276l 12771 1273) . see Fig. Q3]for data from LSCO- 
0.145 (274). For this composition, with T c = 36 K, inelastic 
neutron scattering show that the low-temperature spin gap de- 
creases continuously as function of the applied field and closes 
at about 7T, consistent with a field-induced quantum critical 
point. In fiSR experiments on 214 comp ounds field-enhanced 
magnetism is observed as well d76l : l274[) . 

In contrast to these observations of field-enhanced or- 
der, in compounds with well established stripes the field ef- 
fect in neutron scattering is either weak, as in LBC O-1/8 
d279b . or absent, as in LBCO-0.095 Q, LNSCO-0.12 d274h 
and LNSCO-0.15 d280l) . High-resolution x-ray scattering in 
LBCO-1/8 found a slight field enhancement of the correlation 
length of the c harge order, but no measurable increase of the 
intensity (l28ll) . Note, however, that transport experiments on 
LBCO and LNSCO ha ve been in terpreted in terms of field- 
enhanced charge order (282; 2831). 

Oxygen NMR experiments in fields up to 42 T, mainly 
performed on YBCO, found two components with distinct 
temperatur e dependencies in the spin-lattice relaxation rate 
(12601 l26ll) . The results were interpreted as a enhanced anti- 
ferromagnetic fluctuations in the vortex cores, and have been 
linked to the enhanced magnetism in neutron scattering, de- 
scribed above. 
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FIG. 14 Left: Field tuning of magnetic order in LSCO (reprinted with permission from Ref. J274[) . copyright 2008 by the American Physical 
Society). Shown are the field dependences of the elastic neutron scattering intensity at Q s for LSCO with x = 0.105, x = 0.12, x — 0.145, 
and LNSCO-0.12. The solid lines are fits to the theory of Demler et al. d44lh . see Eq. j I St in Sec. lV.F.ll Moment sizes are plotted relative 
to LNSCO-0.12; the moments at zero field were estimated from /iSR and the gray colors indicate the corresponding error. The LSCO-0.145 
data are presented in arbitrary units. Right: Low-energy LDOS, integrated between 1 and 12meV, of slightly overdoped BSCCO in a field of 
5 T, as measured by STM (from Ref. dl96l) . Science 295, 466 (2002), reprinted with permission from AAAS). A checkerboard-like modulation 
around each vortex is clearly visible. 



A clear-cut neutron scattering observation of field- 
enhanced static order in YBCO is that of Haug et al. d284l) 
in de-twinned YBCO-6.45. At 2K, the quasieleastic inten- 
sity at incommensurate wavevectors was found to increase by 
a factor of 2 from zero field to 15T, with the intensity in- 
crease 1(B) being approximately linear. In turn, the inelastic 
response between 2.5 and 4meV decreased with field. From 
/iSR, the zero-field ordered moment in YBCO-6.45 is rather 
small, about 0.05 fi B at 2K, indicating that YBCO-6.45 is 
closely located to a magnetic QCP. 

A search for field-induced ordering with the additional ad- 
vantage of spatial resolution has been performed using scan- 
ning tunneling microscopy. Hoffman et al. (Il96l) observed 
a checkerboard modulation in the local density of states in 
slightly overdoped BSCCO in an applied field of 5 T, Fig. [14] 
The LDOS, integrated over an energy window between 1 and 
12 meV, showed a modulation with period 4.3 an around each 
vortex, with a decay length of approximately 8 an, signifi- 
cantly larger than the vortex core size. 

As with the zero-field modulation, it is interesting to check 
whether the spatial periodicity of the field-induced modula- 
tion is energy-dependent or not. For BSCCO, Levy et al. 
(12081) reported an energy-independent modulation with p eriod 
(4.3 ± 0.3) ao and fourfold symmetry. Matsuba et al. d209l) 
detected a small breaking of the C4 symmetry towards C2, and 
furthermore emphasized that the modulations at positive and 
negative bias are anti-phase. In optimally doped YBCO, STM 
experiments in a field of 5 T also found signatures of charge 
order with a period of (4.25 ± 0.25) a d21 lb . While these 
experiments suggest field-induced CDW order near the vor- 
tices, a recent detailed analysis of quasiparticle interference in 



a field in CCOC (HH) was interpreted differently: The field 
enhancement of the QPI vectors q% and 05 may also induce a 
checkerboard-like structure in the LDOS. 

At this point, a crucial difference between the spin and 
charge channels in an applied field has to be emphasized: 
While static order in the spin channel is necessarily accompa- 
nied by spontaneous symmetry breaking, charge fluctuations 
are pinned by the vortices which break the lattice translation 
symmetry and couple linearly to the CDW order parameter, 
similar to impurities. Hence, incommensurate fluctuations (in 
either the spin or the charge channel), which would be gapped 
in the absence of vortices or impurities, cause weak inhomo- 
geneous modulations in the charge channel (i.e. of both the 
density and the LDOS) due to pinning. In such a case, the 
LDOS modulations are not expected to be energy-independent 
(HHHH), seeSec.EH 

Summarizing, in superconducting compounds where spin 
or charge order in zero field is weak or absent (but "almost" 
ordered), an applied magnetic field operates to stabilize the 
order, whereas fully developed stripe order displays little field 
dependence. In principle, the physical origin of these phe- 
nomena could be either a direct coupling of the field to the 
spin/charge order parameter, or an indirect coupling where the 
field-induced vortices in the superconducting order parameter 
shift the balance between superconductivity and stripe order 
towards stripes. The strength of the field effect and the nearly 
linear increase or the order parameter with field point towards 
the second, indirect, possibility, because effects of a direct 
coupling between field and CDW/SDW would be quadratic in 
the field. Indeed, a phenomenological theory based on com- 
peting order parameters d44lt 1442 ) account for a large part of 
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the data, see Sec. lV.F.ll and Fig. [14] 

Finally, I note th at, for electr on-do ped cuprates, reports 
both pro (287; 288) and contra d289h field-enhanced mag- 
netism are in the literature. However, in all cases the spin cor- 
relations were found at the commensurate wave vector (w, 7r). 



F. Glassy behavior 

Glassy behavior, usually arising from a manifold of low- 
lying many-body states which are separated by large energy 
barriers, has been seen in a variety of underdoped cuprates at 
low temperatures. While canonical spin-glass physics arises in 
spin systems under the influence of both disorder and geomet- 
ric frustration, the situation in cuprates is more involved due 
to the presence of charge carriers. 

A spin-glass phase has been identified early on in LSCO 
in the insulating doping range 0.02 < x < 0.055 d290t l29lt 
|292|) . Irreversibility and glassy freezing, as evidenced e.g. in 
magnetization and ac susceptibility, are seen below tempera- 
tures of 5 ... 10 K. Remarkably, similar glassy behavior is also 
seen for larger doping, wher e spin-gla s s ord er appears to co- 
exist with superconductivity (12491 12501 : 12931) . Taking together 
NMR and NQR results with those of neutron scattering, one 
concludes that this phase is not a "simple" spin glass with 
little spatial correlations, but instead displays well-defined in- 
commensurate antiferromagneti c cor relations. This has been 
dubbed "cluster spin glass" d52tl290h . 

In fact, the incommensurate magnetic order in the 214 com- 
pounds as well as in underdoped YBCO is characterized by a 
wide regime of very slow fluctuations. The best indication is 
that the apparent ordering temperatures extracted from differ- 
ent probes are significantly different: In elastic neutron scat- 
tering, the onset of quasi-static order is often seen at temper- 
atures by a factor of two (or more) higher than the onset of 
order in uSR. More over, the ordered state is often character- 
ized d98t ll 121 : Il63l) by only medium-range magnetic corre- 
lations and sometimes even by almost isotropic spin orienta- 
tions (dubbed "central mode" behavior by Stock et al. dl63l) 
for YBCO). All these properties point towards the glassy na- 
ture of the incommensurate order. The slowing down of spin 
fluctuations in the low-temperature lim it has been arg ued to 
persist essentially up to optimal doping (293; 294; 295). Spin- 
glass-like behavior was also reported in CCOC in the doping 
range . 5 < x < . 1 2 ; for x — . 1 2 it co-exists with supercon- 
ductivity d2l4 . 

A number of observations have also been attributed to 
glassy behavior in the charge sector. In the insulating regime 
of LSCO, resistivity noise measurements suggest the exis- 
tence of a charge glass, but only below 0.3 K, i.e., far be- 
low the spin freezing temperature d296l) . This is sup porte d 
by detailed impedance spectroscopy in LSCO-0.03 (297). 
For stripe-ordered LESCO, NMR experiments d255l) revealed 
a distinct slowing down of the spin dynamics below 30 K 
(whereas glassy spin freezing only appears at 5 K) - this was 
attributed to a glassy "stripe liquid" forming at 30 K, which 
in turn slows down the spins. In addition, resistivity measure- 
ments in LSCO in both the insulating and superconducting 



regimes have been interpreted in terms of glassy behavior in 
the charge sector, see Ref. (1295b for a discussion. 

Finally, the STM results of Kohsaka et al. (f37b suggest 
the existence of a period-4 valence-bond glass in BSCCO and 
CCOC, based on the fact that the observed order is static, but 
short-ranged. However, dynamical aspects of glassy behavior 
were not studied there. 

On the theory side, relatively little work has been done to 
investigate these collective glass states. Glassy behavior in 
the charge sector is not unexpected, given the presence of 
strong collective effects in the spin-singlet channel, i.e., va- 
lence bonds and stripes, which couple linearly to disorder of 
potential-scattering type. The vast amount of literature on 
the random-field Ising model is relevant in this context, see 
Ref. d62l) for an overview. However, the dynamics in the 
quantum regime is not well understood. The interplay of sin- 
glet formation and disorder has recently been investigated in 
a large-A^ framework and shown to lead to a distinct valence- 
bond glass phase d298l) . 

The origin of glassy behavior in the spin sector is clearly re- 
lated to carrier doping, as the underlying square lattice struc- 
ture of cuprates is unfrustrated. In the insulating regime at low 
doping, a picture dl 18l) of holes that are localized at random 
position s and induce local spiral distortions of the spin back- 
ground (117; 299) appears appropriate. At higher doping, the 
situation is less clear. However, it is conceivable that ordering 
tendencies in the charge sector, which are short-range due to 
quenched disorder, plus competing interactions cause frustra- 
tion in the spin sector, driving it glassy. 

I note that behavior reminiscent of randomness-induced 
glassiness can even arise in the absence of quenched disorder, 
i.e., purely from interactions. This remarkable phenomenon 
has been theoretically demonstrated in Ref. d64l) in a sim- 
ple stripe model with competing (i.e. frustrated) interactions, 
but i s poo rly understood in general. The NMR results of 
Ref. d293l) for LSCO-0.12, showing spin-glass features which 
depend only weakly on the level of chemical disorder, may 
point in this direction as well. 

In summary, both the spin and charge order observed in 
underdoped cuprates display glassy features, which, however, 
strongly mix with collective effects. More theoretical work is 
needed, e.g., in order to understand how the glassy behavior 
influences the low-temperature magnetic excitation spectrum 
measured in neutron scattering. 



IV. MICROSCOPIC MECHANISMS 

After the survey of experimental evidence for lattice sym- 
metry breaking in cuprates, we now turn to the theory side. 
The most important questions appear to be: (i) Which mi- 
croscopic ingredients (on the level of model Hamiltonians) 
are required to obtain symmetry-broken phases as observed 
in experiment? (ii) Which general principles can be identified 
as driving forces of such symmetry breaking? (iii) Beyond 
symmetries, concrete theoretical results for observables (e.g. 
neutron or STM spectra) have to be contrasted with experi- 
mental data. In this section, I shall address (i) and (ii), cover- 
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ing mainly stripes (but briefly also spiral and nematic states), 
while (iii) is subject of Sec.[Vl 

Remarkably, early mean-field studies of the Hubbard model 
predicted the formation of inhomogeneous states at small dop- 
ing, with stripes being one variant. These theory papers, pub- 
lished in 1989-1990, preceded the 1995 experimental reports 
on evidence for stripes, with an enormous subsequent growth 
of theory activities. Naturally, this section can only give a par- 
tial coverage, and I refer the reader to previous review articles 
for alternative expositions {THE!; EH)- 



A. Microscopic models: Hubbard and t-J 

Based on band structure calculations, the low-energy elec- 
tronic properties of the cuprate superconductors are com- 
monly assumed to be dominated by the Cu 3d x 2_ y 2 and O 
%Px,y orbitals of the Q1O2 planes, Fig. Q] Supplementing 
the one-particle terms by local Coulomb repulsion results in 
a three-band Hubbard or "Emery" model d300l) . At half fill- 
ing, the system is a charge-transfer insulator, with one hole 
per Cu orbital and filled O orbitals. 

It is commonly assumed 37 that under certain conditions, 
namely the hopping matrix elements being smaller than both 
the charge transfer energy and the on-site Coulomb repulsion, 
the low-energy physics of such a model can be mapped onto a 
simpler one-band Hubbard model d302l) . 

h = - *y( c L c ^ + c ]o c ™) + U Y1 ' < 12 ) 

defined on the square lattice formed by the copper sites. The 
Hubbard repulsion U is roughly given by the charge-transfer 
energy of the original CuC>2 system; for the cuprates U = 8t 
is a common choice, where t is the nearest-neighbor hopping. 

For large Coulomb repulsion, U ^> t, the one-band 
Hubbard m odel can in turn be mapped onto a t-J model 
(13031 l304l) . with J i3 = 4t?-/Z7. In this mapping, doubly oc- 
cupied sites are excluded from the Hilbert space, and virtual 
hopping processes between neighboring sites transform into 
an antiferromagnetic Heisenberg exchange, leading to 

h = - £ t« (aU«r + h .c.) + Ju 0i ■ Sj - ^) 

(13) 

where the electron operators c* i(J exclude double occupan- 
cies, c| = c' ia (l — n^-u). At half filling, all sites are 
singly occupied, hopping processes are suppressed, and the 
t-J model reduces to the well-known square lattice spin-^ 
Heisenberg model. The mapping from the Hubbard model 
to the t-J model can be understood as expansion in t/U; 



Ref. 1 301) has shown that an exact mapping from the three-band model to a 
one-band Hubbard or t-J model is not possible even in suitable parameter 
limits. This result, however, does not imply that the low-energy physics of 
three-band and one-band models is necessarily different. 



we note that higher orders in this expansion generate addi- 
tional terms in Eq. (TTJt like ring exchange and spin-dependent 
three- site hoppin g which are in general not negligible, see e.g. 
Refs. (13051 13061) . Very accurate perturbative mappings up to 
high order can be obtained using the method of continuous 
unitary transformations; this has been done for the Hubbard 
<— > t—J mapping in Ref. (l306h . At half filling, a direct map- 
ping from the three-band Hubbard model to the Heisenberg 
model has also been discussed d307t) . It is worth emphasizing 
that the discussed mappings imply transformations for all op- 
erators, i.e., care has to be taken when calculating observables 
with one of the effective models. 

In general, models like (fT2b and ([T3l have to be supple- 
mented by longer-range Coulomb interactions, as metallic 
screening can be expected to be poor in particular close to 
the Mott insulating phase. However, such terms are often ne- 
glected for simplicity. 

At present, an open question is how much of the cuprate 
physics is captured by one-band models of Hubbard or t-J 
type. This concerns, e.g., the fundamental issues of supercon- 
ductivity and of the strong asymmetry between hole-doped 
and electron-doped materials. Numerical methods applied to 
the one-band models have reproduced a number of salient fea- 
tures of the cuprate phase diagram (see e.g. Ref. d308T0 . in- 
cluding superconductivity, but the results also indicate that 
other properties are not properly reproduced, see Sec. IIV.CI 
below. However, other workers have raised doubts regard- 
ing the existence of supercond uctiv ity in the one-band mod- 
els, based on numerical results d309l) not showing sufficiently 
strong pairing. 38 While this debate is open, physics of the 
three-band model will eventually be important to fully under- 
stand cuprate properties. 39 

B. Stripes: Concepts and mean-field theories 

Nearly in parallel, Zaanen and Gunnarsson dl5h . Poilblanc 
and Rice (Qil, Schulz 03), and Machida (fl8l;l312l) employed 
mean-field approximations to the one-band and three-band 
Hubbard models to study the physics slightly off half-filling. 
Their real-space Hartree-Fock approaches allowed for inho- 
mogeneous charge and spin distributions, with mean fields 
representing on-site densities. The results showed a clear 
tendency of holes preferring to agglomerate in the presence 
of an antiferromagnetic background. The most favorable 
configurations were such that the holes formed parallel one- 
dimensional lines, later on dubbed "stripes". The Hartree- 
Fock stripes have a number of interesting properties: (i) The 



Based on variational Monte-Carlo studies showing superconductivity in 
the Hubbard model, Ref. ( 310) suggested that the absence of superconduc- 
tivity in the simulations in Ref. i309l) is in fact related to the neglect of 
second-neighbor hopping t' and the large doping values of £ > 0.18 stud- 
ied in Ref. 1303) : In Ref. (31 CD , superconductivity disappears for x > 0.18 
for t' = 0, whereas it survives up to x = 0.25 for t' /t = —0.3 relevant for 
cuprates. 

The importance of three-band physics has also been emphasized in the 
context of more phenomenological approaches, see e.g. Refs. t3ll : l31 ll) . 
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stripes form 7r (or"antiphase") domain walls of the antiferro- 
magnetic order (Fig. |5J. (ii) The number of holes per unit 
length of stripe, pi, is unity, such that a stripe distance of n 
lattice spacings gives a total doping of 1 jn. Those stripes are 
often called "filled stripes" (or "empty" in an electron picture), 
and result in insulating behavior. 

The formation of such stripes reflects a tendency of the 
antiferromagnetic background to expel holes. Indeed, in the 
t-J model at low doping, phase separation into hole-rich 
and hole-poor regions was found d313l) . and hole droplets 
also appeared in numerical studies of the three-band Hubbard 
model d314l) . For the mean-field Hubbard stripes, it was ar- 
gued that quantum fluctuations lead to a bunching of stripes 
(1521) . In all situations, the presence of additional long-range 
Coulomb interactions is expected to disfavor phase separation, 
and stripe-ordered states may result. Such frustrated phase 
separation (HH) has become an important concept in correla- 
tion physics. Subsequent phenomenological theory work sup- 
ported this notion: The presence of interactions with two dis- 
tinct l ength sca les was shown to give rise to stripe-like struc- 
tures iSES). 

The experiments of Tranquada et al, (fl9t I20I) . establish- 
ing stripe order in LBCO and LNSCO, uncovered deficien- 
cies of the early Hartree-Fock stripes. Experimentally, stripes 
were found to be conducting instead of insulating and, consis- 
tent with this, the ordering wavevectors were more compatible 
with half-filled instead of filled stripes. At lowest tempera- 
tures, stripes were even found to co-exist with bulk supercon- 
ductivity. 

As a result, analytical and numerical activities were di- 
rected to find and study metallic and superconducting stripe 
phases. While the numerical works will be summarized in the 
following subsections, I shall briefly discuss alternative mean- 
field approaches here, most of which are based on slave-boson 
formulations of the Hubbard or t-J models. 

Metallic stripes with a filling close to 1/2 were found in 
a Gutz wille r variational treatment of the one-band Hubbard 
model Metallic vertical stripes compete with insulating 

diagonal stripes, and long-range Coulomb interaction is re- 
quired to stabilize the former. Later, a systematic slave-boson 
analysis of the three-band Hubbard model dl33l) found nearly 
half-filled metallic stripes to be stable without long-range 
Coulomb interactions. These stripes were bond-centered, with 
a doping evolution of the stripe distance reminiscent of what is 
found experimentally, i.e., the Yamada plot of Fig. [5] Similar 
metallic stripes were also obtained in mean-field treatments 
of on e-band Hub bard models including second-neighbor hop- 
ping d32lt [322). Lattice anisotropics, as exist in the LTT 
structure of 214 materials, were taken into account in the real- 
space Hartree-Fock study of Ref. (323). As expected, both 
hopping and exchange anisotropies are quite effective in sta- 
bilizing stripe structures. In all cases, the stripe order co-exists 
with well-developed antiferromagnetism in the hole-poor re- 
gions. 

While the mean-field theories sketched so far used decou- 
pling fields on the lattice sites, conceptually different mean- 
field descriptions employ bond var iable s. Those appear in 
the resonating-valence-bond (RVB) d324h and large- mean- 



field theories (HI [H3) of the t-J model. The Sp(2N) 
large- A limit of the t-J model with mean fields be- 

ing an omalous bond v ariables, was used to describe stripes in 
Refs. (HH HH HH). In the Sp(2N) approach, fermionic 
pairing is naturally implemented, while magnetic interactions 
lead to singlet dimerization, i.e., spin-Peierls-like bond or- 
der (111), instead of antiferromagnetic long-range order. In this 
mean-field theory, states with finite doping are superconduct- 
ing at T = 0. The dimerization of the undoped paramagnetic 
Mott insulator survives for a finite doping range, while at large 
doping a homogeneous rf-wave superconductors emerges. At 
small doping, the system is unstable towards phase separa- 
tion. Upon including long-range Coulomb interactions, super- 
conducting stripe states are found, with bond order inherited 
from the undoped system. Thus, the width q and periodicity 
p of the stripes are even, with p and q depending on doping 
and Coulomb repulsion. For t/J values relevant to cuprates 
and q — 2, the stripe filling is close to 1/2, and the doping 
evolution of l/p similar to the Yamada plot. On the basis 
of these results, a global phase diagram for doped Mott in- 
sulators on the square lattice was proposed in Ref. d318h . see 
Fig. [15] Superconducting checkerboard states were investi- 
gated as well d319T) . which turned out to closely compete with 
stripes. Noteworthy, indications for paramagnetic stripes be- 
ing favored over antiferromagnetic stripes were found before 
in numerical calculations for coupled t-J ladders d325l) . 

The evolution from a bond-ordered Mott insulator to a d- 
wave superconductor was investigated in more detail using 
a bond-operator approach d326l) . Ot her m odels for striped 
superconductors have been proposed d327l) . and correspond- 
ing mean-field theories have been used extensively to de- 
scribe concrete experiments, see Sec. [V] More recently, 
renormalized mean-field and variational Monte-Carlo calcu- 
lations have been used extensively to study in detail the real- 
space structure of p aramagnetic superconducting stripe states 
d4l: [3li:[329t[33a) . 

A first-principles view on stripes was obtained from an 
LDA+U calculation of a CUO2 plane, focused on a period- 
4 stripe state (in the charge sector) at 1/8 doping (133 lb . The 
self-consistent stable stripe solution was found to be bond- 
centered and strongly antiferromagnetic with period 8 in the 
spin sector (note that singlet formation cannot be described by 
LDA+U). The calculated photoemission intensity, integrated 
over a large energy window, was found to be roughly consis- 
tent with the early ARPES result of Ref. J242h . 

A common deficiency of most mean-field-like approaches 
to stripes is that the charge modulation within a unit cell is 
very strong, i.e., the hole-poor regions are essentially un- 
doped. Quantum fluctuations can be expected to reduce the 
modulation, see Fig. [16] for a result beyond mean-field. Fur- 
ther, it should be emphasized that the filling of stripes, pi, can- 
not be sharply locked to a particular doping-independent value 
(e.g. pi = 1/2): If this were the case, then it could only be 
the consequence of the existence of an incompressible state at 
this filling, which, however, then would be insulating instead 
of metallic. 

From a one-dimensional perspective, metallic stripe phases 
may display further instabilities: Apart from superconductiv- 
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FIG. 15 Left: Glob al gro und-state phase diagram for doped Mott insulators on the square lattice, proposed in Ref. d3 1 8h (reprinted with 
permission from Ref. ( 318), copyright 1999 by the American Physical Society), as function of doping and the amount of quantum fluctuations 
in the spin sector. The latter is tuned by TV in an Sp(2N) generalization of the spin symmetry or by the strength J'/J of a competing 
frustrating interaction. C (A4) denote the translation (magnetic) symmetries which can be spontaneously broken. All states at finite doping 
are superconducting; insulating Wigner crystal states will occur at very small doping (not shown). The insets illustrate the type of lattice 
symmetry breaking, with circle sizes denoting the amount of hole doping per site, and line widths denoting the strengths of the square-lattice 
bonds. Superconducting stripe states with periodicity p, with p even, occur naturally as a result of hole doping into a paramagnetic dimerized 
Mott insulator. At large doping, the state becomes a d-wave superconductor without additional symmetry breaking. Right: Phase diagram 
emerging from the mean-field treatment of an extended Sp(2N) t-J model (reprinted with permission from Ref. i319h . copyright 2002 by 
the American Physical Society), as function of t/J and the strength of the long-range Coulomb interaction V/t. at a fixed doping of 0.20. 
All states with the exception of the Wigner crystal are superconducting. "Full stripes" refers to states where the charge modulation in the 
large-TV limit is maximal, i.e., the hole density is zero in the hole-poor regions, whereas the "partial stripe" states have a finite hole density 
there. Plaquette (or checkerboard) states generically compete with stripes. 



ity, density waves with 2kp or Akp along the s tripe direction 
may occur. Those have been seen in DMRG d332l) . but not 
experimentally to date. 

Some theory papers have also discussed Wigner-crystal 
states, e.g. Wigner crystals of Cooper pairs (333), in partic- 
ular in connection with checkerboard structures observed in 
STM. However, considering the metallic character of stripes 
in 214 cuprates, this type of approach may be too far on the 
strong-coupling side. 



C. Numerical results 

Numerically exact solutions of microscopic models for the 
CuC<2 planes play an essential role in the cuprate high-T c re- 
search. Naturally, stripe (and other symmetry-broken) states 
have been searched for, with limited success. 

Of the available methods, exact diagonalization is limited to 
system sizes below 40 sites, being too small to detect clear-cut 
stripe signatures. Quantum Monte Carlo (QMC) and density 
matrix renormalization group (DMRG) techniques can treat 
larger systems, with the former usually being unable to access 
low temperatures and the latter being restricted to the ground 
state and to quasi- I d systems. Dynamical mean-field theory 
(DMFT) d34 it. |342|) . extended to a small number of inequiv- 
alent sites, may be used to search for stripes. However, such 
a method does not capture inter-site correlations beyond the 
single-particle level. This problem can be overcome using 
modern cluster extensions (340) of DMFT, provided that the 



clusters are sufficiently large to host stripes. However, these 
computationally demanding method share the QMC problem 
of being restricted to elevated temperatures. In this section, 
we shall discuss results, obtained with the listed techniques, 
relevant for lattice symmetry breaking in cuprate supercon- 
ductors - the focus will be on the existence of and conditions 
for symmetry-broken phases. 

Due to the limited system or cluster sizes, currently all nu- 
merical techniques have severe difficulties in providing de- 
tailed momentum- and energy-resolved spectral information, 
which could be compared to experimental data. Therefore, 
simplified effective models are commonly employed for this 
purpose, with an overview given in Sec.lVlbelow. 



1. Stripes in DMRG 

In an effort to search for stripe states, White and Scalapino 
have applied D MRG to the t-J model on clusters of sizes up 
to 19 x 8 sites (HHHMHH). DMRG, being most suitable 
for the investigation of Id systems, was applied with periodic 
boundary conditions in y direction, but open boundaries in x 
direction. 

The initial calculations d334l) showed stripe states for J/t = 
0.35 and doping 1/8, with a periodicities in the charge and 
spin sector of 4 and 8, respectively, i.e., very similar to the ex- 
perimental data on 214 compounds. Depending on details of 
the boundary conditions, the stripes were either site-centered 
or bond-centered. In all cases, the amplitude of the charge 
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density modulation was about ±40 . . . 50%. Such half-filled 
stripes, i.e. with pi ps 1/2 hole per unit length of stripe, were 
found to be present for all dopings x < 1/8. In contrast, for 
dopings 0.17 < x < 0.3 stripes with pi « 1 were found, and 
the region with 1/8 < x < 0.17 displayed phase separation 
between pi ps 1/2 and pi ~ 1 stripes. 

In Ref. d332h . the calculations were extended to include 
a next-neighbor hopping t 1 and to study pairing along with 
stripe formation. 40 In general, t' /t > was found to en- 
hance d- wave-like pairing while t'/t < had the opposite ef- 
fect, t' with both signs suppressed the tendency toward stripe 
formation, such that stable stripes only occurred for \t'\ < 0.2. 
A plausible interpretation is that stripes compete with super- 
conductivity, but the two may co-exist, with superconductiv- 
ity here b eing two-dimensional rather than one-dimensional 
d332tl338h . Also, the dependence on boundary conditions was 
examined in more detail. For stripes along the x direction, the 
open boundaries at the stripe ends induce a2kp density wave 
in stripe direction. 

Motivated b y the observation of checkerboard-like struc- 
tures in STM (Il35l) . the authors also investigated the pos- 
sible occurrence of checkerboard modulations, i.e., coexist- 
ing vertical and horizontal CDW. Somewhat surprisingly, true 
checkerboards (i.e. with C4 symmetry intact) were not found 
as low-energy states. However, approximate checkerboard 
patterns can arise from stripes with additional CDW forma- 
tion along the stripes d337l) . 

The DMRG results have been criticized for their de- 
pendence of boundary conditions. On the basis of e xact- 
diagonalization results, Hellberg and Manousakis con- 
cluded that the ground state of the t—J model is not striped. 
Instead, stripe states may appear as excited states, which then 
in turn could be favored by suitable boundary conditions or 
lattice anisotropies; the effect of the latter was explicitly stud- 
ied in DMRG in Ref. (US). 

While this debate d332t 13361) has, to my knowledge, not 
been completely settled, the results indicate that tendencies 
toward stripe formation are part of the t—J model physics. 
Note that long-range Coulomb interaction was not included in 
the cited calculations. Hence, the scenario of frustrated phase 
separation (see Sec. HV.Et does not literally apply. 

A disturbing feature of the DMRG results is that no robust 
pairing is obtained for the physical parameter regime oft'/t = 
—0.3. This may cast doubts on whether the one-band t—J 
model contains the essential ingredients required to describe 
cuprate superconductivity. 

Subsequently, DMRG has been applied to study stripes in 
the Hubbard model. A large-scale study of six-leg Hubbard 
ladders d339h investigated stripe states, modulated along the 
leg direction, for a doping of x = 9.5% at U/t = 3 and 
12. Finite-size scaling indicated that stripes are stable in the 
infinite-length limit for U/t = 12, whereas those at U/t = 3 



The experimentally measured Fermi-surface shapes of hole-doped 
cuprates are reproduced if one assumes t'/t < 0, with t'/t = —0.3 being 
a typical value. 



are artifacts of DMRG boundary conditions. For U/t = 12, 
this result suggests a small value of the Luttinger parameter 
K p and hence only weak superconducting correlations. 



2. DMFT and quantum cluster methods 

While single-site DMFT is designed to treat homoge- 
neous systems, it can be easily generalized to inhomoge- 
neous or modulated situations, by assuming a local, but site- 
dependent electronic self-energy. Then, the lattice problem 
maps onto a set of single-impurity problems supplemented by 
self-consistency conditions. The inhomogeneous DMFT ap- 
proach treats correlation effects only locally, and modulations 
are captured in mean-field fashion. This implies that valence- 
bond physics is not part of this method, and tf-wave supercon- 
ductivity cannot be described beyond mean-field (in contrast 
to cluster DMFT to be discussed below). 

Fleck et al. (1343b have applied an inhomogeneous DMFT 
method, using supercells up to 36 x 8 sites, to the single-band 
Hubbard model in the absence of superconductivity. They 
found metallic stripes for doping levels 0.03 < x < 0.2, 
with a number of remarkable properties: The stripe orien- 
tation changed from diagonal to horizontal/vertical around 
x = 0.05, and the stripe separation followed closely the Ya- 
mada plot. For x < 0.17, the stripes were site-centered - this 
is not surprising given the absence of bond singlet physics 
in the approach. Also, nodal quasiparticles were gapped by 
the stripe order; while this was consistent with some early 
ARPES results d242l). it i s inconsistent with more recent data 
HEM El |245t|24$. 

The problem of missing non-local correlations is overcome 
in powerful cluster extensions (13401) of DMFT, with the dy- 
namic cluster approximation (DCA) and the cluster DMFT 
(CDMFT) being the most prominent variants. DCA is based 
on coarse-graining the momentum dependence of electronic 
self-energy, and the lattice problem is now mapped onto a 
correlated cluster embedded in a dynamic bath in the single- 
particle sector. In CDMFT, the cluster is formed in real 
space instead of momentum space. Solving the cluster impu- 
rity problems, typically using QMC or exact diagonalization 
methods, is numerically much more demanding than treating 
a single impurity, which limits the applicability of the method. 
Most practical calculations use N c = 4 cluster sites - the 
minimum required for <i-wave superconductivity - but calcu- 
lations up to N c = 64 (albeit restricted to high temperatures) 
have also been reported. Although the results should become 
increasingly reliable with growing cluster size, the mean-field 
character of the method has to be kept in mind, i.e., the self- 
consistency loop will converge to dynamic saddle points of 
the problem. 

DCA and CDMFT have been quite successfully applied 
in the study of cuprate superconductivity (308). Starting 
from the standard 2d Hubbard model, supplemented by a 
nearest-neighbor hopping t', these methods have reproduced 
key features of the cuprate phase diagram as function of dop- 
ing and temperature: Antiferromagnetism and superconduc- 
tivity occur in the roughly correct doping and temperature 
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FIG. 16 Bond-centered stripe state in the t—J model as obtained from DMRG (reprinted with permission from Ref. d335l) . copyright 1999 by 
the American Physical Society). Left: Hole density and spin moments for the central 8x8 region of a 16 x 8 t—J system. The diameter of the 
gray circles is proportional to the hole density rih, and the length of the arrows is proportional to (S z ), according to the scales shown. Right: 
Hole density (solid circles) and spin structure function Sir (open squares), averaged over each column on the 16 x 8 system. 



ranges (I344l>. a pseudogap occurs in the underdoped regime 
(1341 l346t l347h . and the normal state around optimal dop- 
ing is character ized by i ll-defined quasiparticles and a large 
scattering rate 13491) . The dichotomy between nodal and 
antinodal quasiparticles in the underdoped regime, i.e., well- 
defined nodals and broad incoherent antinodals, has been dis- 
cussed in detail (350). 

So far, states with broken lattice symmetries have not been 
unambiguously identified in DCA or CDMFT. Clearly, the 
smallest cluster size of 2 x 2 is insufficient for stripes, as 
non-magnetic stripes require at least a 4 x 4 cluster (and 
low temperatures). What has been found are strong tenden- 
cies towards bond order in the Hubbard model: The parent 
Mott insulator with additional magnetic frustration displays 
a dimerized phase in CDMFT on a 2 x 2 cluster, but also 
the neighboring phases are characterized by soft singlet-to- 
singlet fluctuations (1352b . Using DCA, the doped Hubbard 
model has been shown to display a bond-order susceptibil- 
ity which diverge s as T — > in the underdoped regime of 
doping x < 0.22 d347l) . This has been taken as evidence for 
a quantum critical point at optimal doping between a Fermi 
liquid and a bond-ordered phase. While the precise interpre- 
tation of these numerical results may be problematic due to 
finite-size and finite-temperature effects, they support the no- 
tion that valence-bond physics is a relevant player in under- 
doped cuprates. 

The self-energy-functional theory of Potthoff allows 
to construct more general dynamical cluster theories. Here, 
correlated clusters coupled to uncorrelated bath sites, all 
with variationally determined parameters, are solved exactly. 
DMFT and CDMFT arise in the limit of an infinite number of 
bath sites. In contrast, without bath sites one obtains the so- 
called variational cluster approach ( VCA) (354). Compared to 
CDMFT, larger clusters can be used and access to low temper- 
atures is simplified. This m ethod has been e mployed recently 
to the 2d Hubbard model d354t 13551 13561) . and the general 



structure of the phase diagram has been obtained in overall 
agreement with DCA and CDMFT. So far, stripe states have 
not been systematically searched for. A recent extension al- 
lows for the calculation of two-particle correlations beyond 
the bubble approximations d357l) . and results will be discussed 
below. 

Interestingly, the present results from quantum cluster ap- 
proaches indicate that the difference between electron and 
hole doping in the cuprates is not fully captured by the single- 
band t-t' Hubbard model: Although antiferromagnetism is 
found to be more stable for electron doping, in agreement with 
experiment, the low-temperature amplitude of the supercon- 
ducting d-wave order para meter is roughly equal for electron 
and hole doping (308l: l356h . in disagreement with experi men t. 
An initial DCA study of the three-band Hubbard model (l35ll) . 
with microscopic parameters extracted from LDA bandstruc- 
ture calculations and limited to small cluster of four Cu sites, 
has indeed shown that the superconducting T c depends sensi- 
tively on microscopic parameters. 



D. Weak-coupling approaches 

For small band filling, or equivalently large doping, it is 
natural to expect weakly correlated Fermi liquid behavior in 
the Hubbard and t—J models. Assuming that this persists up to 
band fillings of ~ 0.7, weak-coupling many-body techniques 
may be used to access the physics of high-T c cuprates from 
the overdoped side. 



1 . Random-phase approximation 

The simplest way to capture collective instabilities of a 
Fermi liquid is the random-phase approximation (RPA). This 
has been primarily employed to describe the spin excitations 
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in hole-doped cuprates and their interplay with superconduc- 
tivity. RPA is an infinite resummation of bubble diagrams for 
the susceptibility, leading to 



1 - 9(q)Xo(g,v) 



(14) 



where xo is the Lindhard susceptibility of the non-interacting 
systems and g is the fermionic four-point vertex in the singlet 
(triplet) channel for the charge (spin) susceptibility. 

RPA can predict instabilities towards ordered phases: 
Xrpa (q, ui = 0) diverges when g(q)xo (<f, w = 0) reaches unity 
at some wavevector q — Q. In addition, RPA is able to de- 
scribe collective modes: Sharp modes are reflected in poles of 
Xrpa at some finite frequency, which then have to lie outside 
the particle-hole continuum of Xo> whereas damped modes 
can exist inside the continuum. In fact, RPA has been initially 
used to model the emergence of the so-called spin resonance 
mode in the superconducting state of cuprates: This mode at 
(tt, 7r) is overdamped in the normal state, but pulled below 
the particle-hole continuum by the onset of superconductiv- 
ity, see Sec. IV.CI 

For an interaction vertex g with little structure in mo- 
mentum space, RPA predicts the instability wavevector Q to 
be given by the maximum of xoiQi 1 ^ = 0), i.e., Q is de- 
termined by Fermi-surface properties. From the measured 
Fermi surfaces of underdoped cuprates, there are two candi- 
date wavevectors Q where xo displays peaks: (i) a wavevec- 
tor close to (tt, 7r), connecting regions near the nodal points 
- this has frequently been related to incommensurate SDW 
order, and (ii) a short wavevector in (1,0) or (0, 1) direction 
which connects the often nearly straight antinodal segments 
of the Fermi surface and which may be related to CDW or- 
der. With increasing hole doping, the magnitude of both vec- 
tors decreases. For the SDW nodal nesting wavevector this 
means that the difference to (tt, tt) increases, i.e., the mod- 
ulation period w.r.t. commensurate antiferromagnetism de- 
creases, which is in qualitative agreement with spin fluctua- 
tions observed in various cuprates. In fact, RPA has been used 
quite successfully to describe features of the cuprate spin fluc- 
tuation spectrum, with a brief overview given in Sec. IV.CI In 
contrast, the possible CDW period deduced from the antin- 
odal nesting wavevector increases with hole doping, in sharp 
contrast to ordered stripes in 214 cuprates. As charge order 
sets in at a higher temperature than spin order (see Sec. IIII.Al i, 
this indicates the failure of RPA in the underdoped regime, at 
least for the 214 compounds. (Note, however, that a CDW 
wavevector decreasing with doping has been extracted from 
STM measurements in BSCCO-2201 J207h .) A failure of 
RPA has also been pointed out f or th e dynamic spin excita- 
tions in optimally doped YBCO ( fl70h and BSCCO (fl65l) as 
well as in electron-doped Pro.8sLaCeo.i2Cu04 based 
on strong mismatch between RPA and experiment regarding 
the spin fluctuation intensities and the structure of the signal 
above T c , A more detailed critical discussion of the applica- 
bility of RPA is in Sec. IVHl 

In applications to the Hubbard model a common choice 
is g = U, whereas for the square-lattice t—J model g = 



J(cosk x + cosfcj,). In practical applications, the amplitude 
of g is often treated as a fit parameter to adjust the position of 
the instability, diminishing the predictive power of RPA. 

Although it is frequently assumed that RPA is a controlled 
approximation for weakly interacting systems, this belief is 
incorrect: In situations with competing instabilities, RPA is 
unreliable in principle even at weak coupling. This can be eas- 
ily seen for the half-filled Hubbard model with perfect nesting, 
where the U dependence of the antiferromagnetic transition 
temperature is predicted incorrectly by RPA. As one of the 
main features of cuprates is the presence of multiple instabili- 
ties, approaches beyond RPA are required. 



2. Functional renormalization group 

A systematic and consistent treatment of competing weak- 
coupling instabilities is provided by the functional renormal- 
ization group (fRG) treatment of the underlying microscopic 
model. The practical applications of the fRG method to the 
Hubbard model have been summarized in detail by Hon- 
erkamp (HH). 

The starting point is an exact formulation of the renormal- 
ization group in terms of flow equations for A^-particle irre- 
ducible vertex functions, which are truncated such that only 
the flow of the one-particle vertex functions is kept. To solve 
the flow equations numerically, the frequency dependence of 
the vertex functions is neglected, and the momentum depen- 
dence is discretized via suitable patches in momentum space. 
Different flow schemes have been developed, with either the 
momentum-space cutoff, the temperature, or the interaction 
as flow parameter (13581) . In the temperature-flow scheme, one 
starts at high temperatures T and calculates the temperature 
dependence of the momentum-resolved interaction vertices 
in the particle-particle and particle-hole channels. Some of 
these vertices will grow upon reduction of T, with a diver- 
gence indicating an instability towards an ordered state. Usu- 
ally, this stops the applicability of the fRG scheme. 41 

Importantly, the fRG method is able to deal with compet- 
ing ordering phenomena in an unbiased manner, in contrast to 
RPA. The perturbative character of the fRG limits its applica- 
bility to small and moderate values of U, practically values of 
U/t up to 4 have been used. It furthermore implies that the 
physics of the Mott insulator cannot be described, although 
precursors of Mott physics driven by strong Umklapp scatter- 
ing have been discussed (13581) . 

The pub l ished fRG calculations for the 2d Hubbard model 
d360t l36lt l362t l363t |364|) consistently show antiferromag- 
netism (at small doping) and <i-wave pairing (at larger doping) 
for both signs of the second-neighbor hopping t', correspond- 
ing to electron and hole doping, respectively. The hole-doped 
regime with t 1 ft < is addition characterized by a so-called 
saddle-point regime where many instabilities compete and re- 



41 A proposal to circumvent this problem is the application of a sma ll 
symmetry-breaking field which is kept as additional flow parameter )359l) . 
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enforce each other, in addition to antiferrmagnetism and d- 
wave pairing these are a d-wave Pomeranchuk (nemati c) in - 
stability and, to a weaker extent, d-density wave order d363h . 
In the overdoped regime for 0.15 < x < 0.30, a recent fRG 
calculation (13651) found strong angle-dependent scattering at 
U/t = 4. This leads to an apparent violation of Fermi-liquid 
behavior, consistent with transport experiments on overdoped 
T^I^CuOg+z (366), but does not exclude that the underly- 
ing low-temperature state is a Fermi liquid. 

In present fRG studies, instabilities toward charge order are 
found to subleading. Interestingly, a recent study d367l) of a 
Hubbard-like model supplemented by electron-phonon inter- 
actions found that phonons strengthen the tendency towards a 
CDW with a d-wave form factor and a wavevector given by 
the nested antinodal Fermi-surface pieces. However, also this 
instability appears to be magnetism-driven. In summary, gen- 
uine CDW instabilities do not seem to occur in the Hubbard 
model at weak coupling. 



E. Origin of stripe formation 

After having reviewed results from microscopic calcula- 
tions, I will try to summarize the current phenomenological 
ideas on the origin of stripes. This is particularly important, 
as collected experimental and theoretical results described in 
Sees. [Ill] and [IV] leave us with what looks like a conflict: On 
the one hand, stripes have been unambiguously identified in 
experiment, but, on the other hand, stripes do not readily 
appear in theoretical treatments (beyond mean-field) of the 
popular microscopic models. An incomplete list of possible 
explanations is: (i) Without the LTT distortion of some 214 
compounds, stripes are never the lowest-energy state of CuC<2 
planes. While this is certainly a possibility, tendencies toward 
stripes have also been identified in cuprates with tetragonal 
CuC»2 planes, and hence should be detectable in numerics, (ii) 
The present limitations inherent to the numerical methods, i.e. 
cluster size, temperature etc., do not yet allow to detect stripes, 
(iii) The investigated single-band models do not contain all in- 
gredients required for stripe formation - this appears to be a 
likely explanation. 

About the origin of stripes: Phenomenologically, popular 
lines of thought - all from a strong-coupling perspective and 
not being mutually exclusive - have been (A) frustrated phase 
separation (or micro phase separation), (B) spin-charge "topo- 
logical" properties, and (C) valence-bond solid formation, all 
of which have been invoked to rationalize the formation of 
conducting stripes in cuprates. While (A) is general, (B) and 
(C) refer to more microscopic aspects. 

Frustrated phase separation 1521 13151 : 13161) requires little 
explanation: A system of particles on a lattice, moving un- 
der the influence of short-range attractions and long-range re- 
pulsions, minimizes its energy by forming linear domains of 
enhanced and reduced density. 42 This argument invokes long- 



Ref. (368) has ai'gued, based on a Ginzburg-Landau analysis, that long- 



range Coulomb repulsion as crucial, and leaves open the issue 
of how the short-range physics exactly looks like. Being clas- 
sical, it generically results in conducting instead of insulating 
stripes. 

Microscopically, frustrated phase separation in cuprates is 
related to the fact that an undoped Mott insulator (partially) 
expels holes. 4344 A resulting stripe state then not only bal- 
ances magnetic (i.e. short-range attractive) and Coulomb (i.e. 
long-range repulsive) energies, but also the hole kinetic (i.e. 
quantum) energy: Holes can gain kinetic energy by moving 
along the hole-rich stripes, but at the same time leave the 
hole-poor domains essentially unaffected. Due to the pres- 
ence of quantum effects, long-range Coulomb interaction may 
not be absolutely necessary for stripe formation (although it is 
in some mean-field approaches which underestimate quantum 
effects). 

Inspired by the physics of one-dimensional systems, spin- 
charge "topological" properties have been discussed as driv- 
ing force of stripe formation. Hole-rich stripes form charged 
domain walls of the background antiferromagnetic order, 
which may be inter preted as a two-dimensional generaliza- 
tion of holons d37ll) . Within this picture, stripes form elas- 
tic s trings on the lattice, which may be pinned by disorder 
(HH l372t 13731) . and the destruction of a stripe states is driven 
by transverse stripe fluctuations or topological defects. The 
term "topological" is somewhat misleading here, as domain 
walls are topological defects of Ising antiferromagnets, but 
not of ones with Heisenberg symmetry. Hence, there is no 
topological protection for the strings. Microscopic considera- 
tions of the motion of holes along AF domain walls within the 
t—J model show that motion both along the domain wall and 
into the AF domains (which leads to some spin-polaron-like 
dressing) is required to compensate for the cost of the domain 
wall d374l) . A general problem with this set of ideas is that it 
is intimately tied to magnetic order. Experimentally, charge 
stripes exist without long-range magnetic order (Sec. [TTTb. al- 
though some short-range order is certainly present. 

In contrast, valence-bond solid formation is a concept re- 
lying on physics in the singlet sector. It starts from the ob- 
servation that the destruction of antiferromagnetic order in an 
undoped Heisenberg magnet on the square lattice is accompa- 
nied by lattice symmetry breaking, mos t likely in the form of 
a columnar valence-bond solid (3; 376). It has been argued 
that this ty pe of symmetry breaking survives with hole doping 
(13181: 13751) . leading to a VBS metal or superconductor. More- 
over, the resulting system, effectively consisting of coupled 
two-leg ladders, may display a stripe (or frustrated phase sep- 
aration) instability, because weakly doped ladders have strong 



range physics alone is insufficient to produce the experimentally observed 
antiphase magnetic stripe order. A more formal reasoning on frustrated 
phase separation can be found in Ref. i369l) . 

This is most easily seen by the fact that hole motion in an antiferromag- 
netic background is frustrated due to the creation of spin defects. 

In the phenomenological SO(5) theory, a tendency towards phase separa- 
tion is found as well 1370I) . here with superconducting and antiferromag- 
netic domains. 
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rungs, whereas kinetic energy is better gained by motion along 
the legs. The STM data of Ref. ©, showing bond-centered 
stripes with strong bond modulations, appear to support this 
concept. It should be pointed out that valence-bond stripes, 
although not relying on magnetic long-range order, are pri- 
marily driven by magnetic exchange interactions, as those set 
the scale for VBS formation. 

In both microscopic pictures, the doping dependence of 
the stripe ordering wavevector is explained by the hole-rich 
stripes having a filling close to half a hole per unit length 
of stripe which weakly depends on the overall doping x for 
small x. As a result, Q c w (4nx, 0). For larger x and short 
stripe distances, stripes strongly repel each other, resulting in 
a plateau in Q c (x), before stripes eventually disappear. 

While stripes may in principle arise from weak-coupling 
Fermi-liquid instabilities, many cuprate experiments point to- 
ward a strong-coupling picture being more appropriate, as de- 
tailed in Sec. IV.HI However, it is well conceivable that nearly 
nested antinodal Fermi surface pieces have a share in driving 
charge order. This speculation is particularly appealing as an 
explanation of increased stability of stripes at x = 1/8. This is 
commonly attributed to lattice commensuration effects, may 
be further enforced by Fermi surface nesting. In this scenario, 
x = 1/8 is the doping concentration where the antinodal nest- 
ing wavevector and the strong-coupling stripe wavevector co- 
incide. 



F. Spiral magnetism 

Holes doped into a quantum antiferromagnet ha ve been ar - 
gued to induce an instability towards a spiral state (ll 17b l299h . 
It has been proposed that the physics of spin spirals domi- 
nates the insulating regime of cuprates at low doping (111 8b . 
In this scenario, holes are pinned at random positions by the 
electrostatic potential of the dopant atoms. The local spi- 
ral distortions of the classical spin background lead to elas- 
tic peaks in neutron scattering around wavevectors Q s = 
27r(0.5 ± e' s , 0.5 ± e^), with a doping evolution of e' s as ob- 
served experimentally. Although a doped spiral state is un- 
stable w.r.t. phase separation, this may be circumvented here 
by strong pinning. A distinct but related scenario was pro- 
posed in Ref. (11 19b : Magnetic anisotropies were argued to 
stabilize a canted Neel state at doping x < 2%, whereas they 
lead to a helicoidal magnetic phase at larger doping (but still 
in the insulating phase), with a small out-of-plane magnetiza- 
tion component. This proposal appears more consistent with 
susceptibility measurements in weakly doped LSCO (13771) . 

In the superconducting state at higher doping, mobile holes 
and quantum effects become more important and drastically 
change the underlying physics. Therefore, the physics of spin 
spirals is commonly assumed to be irrelevant to cuprate su- 
perconductivity. However, the smoking-gun polarized neu- 
tron scattering experiment has not been performed to date, see 
SecHirA] 



G. Nematic order 

As with stripes, both weak-coupling and strong-coupling 
approaches have been employed to search for nematic phases 
in correlated electron models. A common feature of the stud- 
ies described below, which make no reference to fluctuat- 
ing stripes, is that a nematic instability (equivalent to a spin- 
symmetric e?-wave, I = 2, Pomeranchuk instability) is driven 
by strong forward scattering. 

Strong-coupling mean-field calculations in the standard 
RVB slave-boson formalism find a nematic phase in the 2d 
t-J model, whi ch compe tes with homogeneous d-wave su- 
perconductivity (378; l379l) . For the 2d Hubb ard model, weak- 
coupling functional RG has been applied (!360| l36ll 13831) . 
with qualitatively similar conclusions: In particular near the 
van-Hove filling, there is a tendency towards a <i-wave Fermi 
surface deformation. Subsequently, antiferromagnetism or d- 
wave superconductivity may develop at low temperatures in- 
side the nematic state. If, however, antiferromagnetism or su- 
perconductivity set in before the nematic order, then the for- 
ward scattering interactions stop to grow, and a nematic insta- 
bility does not develop. The overall picture of nematic order 
competing with d-wave pairing is supported by exact diago- 
nalization d380l) and variational Monte-Carlo d38ll) studies of 
the t-J model. A nematic phase has also been shown to occur 
in a sp ecial strong coupling limit of the three-band Hubbard 
model (13821) . 

In 3d lattice models, the Pomeranchuk instability for even I 
is generically of first order due to the presence of cubic terms 
in the Ginzburg-Landau theory. It has been argued that even 
in d = 2 s trong fluctuations may drive the transition first order 
(l386t 1387b . Full quantum critical behavior is restored in these 
2d models if a sufficiently strong repulsive term is added to 
the forward-scattering interaction (388). 

Mean-field treatments of nematic ordering often employ ef- 
fective models with a quadrupole-quadrupole inte raction ex- 
plicitly d esigned to produce an I = 2 instability d385t l386t 
I387H388I) . Recently, it has been shown that generic central in- 
teractions in 2d can produce Pomeranchuk instabilities as well 
d389h . 

Notably, nematic order in an isotropic Fermi liquid (i.e. 
without lattice) leads to non-Fermi liquid behavior due to the 
overdamped Goldstone fluctuations of the director order pa- 
rameter (13851) . 

A natural expectation is that the formation of stripes should 
be favored in a <i-wave nematic phase due to the uni-axial 
anisotropy (21). However, to my knowledge, a comprehen- 
sive study of a microscopic model (using mean-field or more 
elaborate techniques) with a sequence of disordered, nematic, 
and stripe phases is not available. 



V. EFFECTIVE MODELS: LINKING THEORY AND 
EXPERIMENT 

In this section, we describe theoretical activities aimed on a 
detailed modeling of experimental data which have been ob- 
tained in (or close to) phases with broken lattice symmetries. 
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Given the difficulties with the numerically exact treatment of 
strongly correlated Hubbard or t-J models, most efforts are 
constrained to simplified effective models representing spin or 
renormalized single-particle degrees of freedom. Most of the 
section will again be devoted to stripe order, but nematic and 
loop-current order will also be covered. 



A. Weak modulations vs. coupled chains and ladders 

Often, effective models assume the presence of static stripe- 
like modulations in both the charge and spin sectors. Two 
seemingly distinct viewpoints can then be employed for model 
building: Either (A) one starts from a two-dimensional sys- 
tem and adds (weak) modulations, or (B) one starts from one- 
dimensional chains or ladders and adds a (weak) transverse 
coupling. The states which result from these two approaches 
are usually equivalent regarding symmetries, but for a given 
experimental situation one of the two may be more appro- 
priate. Intuitively, (A) leads to a weakly modulated state, 
whereas (B) corresponds to strong modulations. Experimen- 
tally, the amplitude of the modulations in the charge sector 
is not precisely known, and extracting modulation amplitudes 
from actual data is model-dependent. From x-ray scattering, 
the variation in the oxygen hole concentration in LBCO-1/8 
has been estimated to be about a factor 4 (11321) : the STM tun- 
neling asymmetry maps show a contrast of a factor 1 .5-2 l37h . 
This suggests that the modulations are not weak, but also not 
close to the maximum limit. 

However, this is not the full story: Id building blocks of- 
ten feature fractionalized excitations described by a Luttinger 
liquid, whereas excitations in 2d are conventional. Then, the 
elementary excitations of the starting points of (A) and (B) are 
different. Nevertheless, upon including the perturbations, the 
low-energy physics is the same in both cases, as inter-chain 
coupling is a relevant perturbation of the Luttinger-liquid fixed 
point (except for certain frustrated cases), and even weakly 
coupled metallic chains form a 2d Fermi liquid in the low- 
energy limit. Consequently, the answer as to whether ap- 
proach (A) or (B) leads to a better description of experimental 
data may depend on the observable and the energy range to 
be considered. In principle, Id spinons and holons could be 
a good description of the excitations at elevated energies pro- 
vided that the modulations are strong - this has been proposed 
(1390b in interpreting certain photoemission data, but is contro- 
versial. 

For some observables, the situation is simpler: For instance, 
the spin excitation spectrum of bond-ordered stripes can be 
described both by coupled spin ladders and by a combined 
theory of CDW and SDW order parameters in a 2d system, 
with essentially identical results, see Sec. IV.C.ll One reason 
for this equivalence is that the spin excitations of a two-leg 
ladder are conventional spin-1 triplons. 



B. Static, fluctuating, and disorder-pinned stripes 

The picture of long-ranged fluctuationless stripe order in 
the charge sector is certainly idealized: Even at lowest tem- 
peratures (i.e. without thermal fluctuations), both quantum 
fluctuations and quenched disorder (arising e.g. from dopant 
atoms) will induce deviations from ideal order. Deep in the 
charge-ordered phase, such deviations are often negligible 
(however, e.g. photoemission spectra are strongly influenced 
even by small amounts of stripe disorder). 

The most interesting regime is close to the charge ordering 
transition, where - in the absence of quenched disorder - col- 
lective degrees of freedom associated with the CDW order are 
slowly fluctuating in space and time. Here, the CDW physics 
on short scales is quantum critical, and a comprehensive the- 
oretical treatment becomes difficult. 

In building approximations, the question arises which type 
of stripe fluctuations dominates. In the language of the CDW 
order parameter <f> c , one can distinguish amplitude and phase 
fluctuations. The latter imply that stripes remain well-defined 
objects, but fluctuate in the transve rse d irection - here the 
picture of elastic strings on a lattice 0372c 13731) has been fre- 
quently used. However, this is not the full story: The phase 
field can display topological defects, corresponding to stripe 
end points. Moreover, in a tetragonal environment, there will 
be generalized phase fluctuations between horizontal and ver- 
tical stripes. In high dimensions, it is known that transitions 
are driven by amplitude fluctuations, in 2d the situation is less 
clear, and the importance of certain fluctuations depends on 
whether the system is more in the weak-coupling or strong- 
coupling limit. 

Including the effects of quenched disorder complicates mat- 
ters, but may also imply some simplification (Sec. III. Gi l. Most 
importantly, disorder is of random-field type in the charge 
sector and tends to pin fluctuations. Then, on the one hand, 
short-range ordered CDW configurations are rendered static 
(and hence visible e.g. in STM). Pinning will act on both 
horizontal and vertical stripes and thus tends to smear the dis- 
tinction between stripe and checkerboard order (165c 1661) . On 
the other hand, the finite-frequency CDW dynamics becomes 
glassy, and the charge-ordering transition is smeared. 

In the literature, various routes have been followed to tackle 
the situation of fluctuating stripes: (i) Ignoring subtleties of 
critical physics and disorder, one may treat the CDW modes 
perturbatively, e.g., by calculating the lowest-order diagrams 
descr i bing the interaction of electrons with CDW modes 
(l39ll [392). This approach is suitable if amplitude fluctua- 
tions of the CDW order are dominant, (ii) Theories of quan- 
tum strings on a lattice have been invoked to discuss the phase 
diagram and the dynamics of fluctuating stri pes in the pre s - 
ence of both disorder and lattice pinning (I2lc l3^[373l:[393l) . 
In particular, disorder pinning may lead to a Bragg-glass-type 
state (12 ll) . Based on a calculation of the ratio between spin 
and charge correlation lengths in scenarios with and with- 
out topological string defects, it was argued that stripes in 
214 cuprates are dominated by non-topological fluctuations 
(13931) . (iii) Numerical simulations of order parameter theories 
can in principle account for all possible fluctuations. As the 
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full dynamical treatment of coupled modes is rather compli- 
cated, a Born-Oppenheimer approac h to the coupled SDW- 
CDW problem was proposed in Ref. dl5ll) : It is assumed that 
the collective charge modes are slow compared to the spin 
modes. The CDW modes are approximated as static on the 
time scale of the spin fluctuations, and the spin sector is solved 
for fixed charge configurations. This approach neglects charge 
dynamics, i.e., inelastic processes involving CDW modes, but 
is otherwise non-perturbative and al lows to ac cess the full spin 
dynamics, as discussed in Refs. dl5lt |224[) . The assump- 
tion of slow charge dynamics appears particularly good in the 
presence of strong pinning by quenched disorder. Then, the 
Born-Oppenheimer approximation may even be used down 
to zero energy, e.g. to calculate photoemission spectra, (iv) 
A strongly simplified version of the Born-Opperheimer ap- 
proach consists of treating static parallel stripes with random 
spacings, which only accounts for a very special type of CDW 
phase fluctuations. 

While these approaches provide some insight into different 
limits of the fluctuating-stripe problem, a concise theoretical 
treatment of the random-field glassy dynamics in the quantum 
regime is not available to date. 



C. Spin excitations 

Magnetic excitations in doped Mott insulators such as 
cuprates can either be calculated numerically from one of the 
relevant microscopic models (with the limitations discussed 
in Sec. HV.Ct or can be obtained from simplified effective de- 
scriptions: Here, magnetism can be modeled in an "itinerant" 
or in a "localized" concept. In an itinerant description, one 
starts from band electrons and captures interaction-generated 
collective dynamics via RPA or more sophisticated methods. 
Alternatively, in a localized picture, the modeling is centered 
around collective modes, typically obtained from Heisenberg 
models of localized spins, and the physics of mobile carriers is 
added perturbatively or ignored entirely. While the latter op- 
tion implies a quite drastic simplification, it may be sufficient 
to describe collective modes (seen as peaks in x" ) which re- 
side on top of a broad continuum of particle-hole excitations 
- this continuum is then simply not part of the model. 

The itinerant description is obviously more appropriate at 
small band filling or large doping, where correlation effects 
are expected to be weak, whereas local-moment collective 
magnetism should prevail in the small-doping regime. In the 
absence of symmetry breaking, the two situations may or may 
not be adiabatically connected - theoretical scenarios for both 
cases are known: In a one-band Hubbard model at half-filling 
and close to nesting, the small-C/ itinerant antiferromagnet 
is continuously connected to the large-C/ local-moment an- 
tiferromagnet d394l) . On the other hand, in two-band (Kondo- 
lattice) models of interacting electrons, metallic paramagnetic 
phases with and without local moments can exist, which are 
not continuously connected: While the latter is a conventional 
Fermi liquid, the for mer i s an exotic state dubbed "fraction- 
alized Fermi liquid" (395). For the cuprates, it is not known 
whether the itinerant and localized descriptions of magnetism 



are adiabatically connected. As detailed below, various ex- 
periments in doped compounds can in principle be described 
within both concepts, however, the standard RPA (or "fermi- 
ology") models fail in some important cases. 



1 . Spin excitations of ordered stripes 

Static stripe order in both the spin and charge sectors, which 
has to be accompanied by linearly dispersing low-energy mag- 
netic excitations due to Goldstone's theorem, 45 lends itself to 
relatively simple theoretical descriptions. 

In particular, in a picture of spin moments localized in the 
hole-poor regions of the square lattice, one can write down a 
Heisenberg model for coupled spin ladders. 

(ij) (ik) 

with the Y^Uj) running over intra-ladder bonds, while the sum 
connects neighboring sites from different ladders. The 
model may be supplemented by longer-range couplings or 
ring-exchange terms. Microscopically, J is the usual antifer- 
romagnetic exchange between neighboring Cu spins, whereas 
J' is an effective exchange which is thought to be mediated via 
the hole-rich stripes, Fig. [17] Note that the sign of J' deter- 
mines whether ordered antiferromagnetism will be in-phase or 
anti-phase. The hole-rich stripes themselves are not part of the 
model; they are assumed to only contribute a broad continuum 
to the spin excitations, and their damping effect is ignored. 

For site-centered stripes with a period 4 in the charg e sec - 
tor, the canonical choice are three-leg ladders dl50l : 13981) . 
"Antiphase" magnetic order, with an ordering wavevector of 
Q s = (37r/4,7r), is obtained from antiferromagnetic inter- 
ladder coupling, J' > 0. For bo nd-centere d period-4 stripes, 
two-leg ladders appear natural dH3 [Hi, but four-leg lad- 
ders have also been used ( 150); anti-phase magnetism requires 
J' < in both cases. For odd-leg ladders, an infinitesimal 
J' is sufficient to drive the coupled system into a magneti- 
cally ordered ground state, while a finite J' is required in the 
even-leg case in order to close the spin gap of the isolated lad- 
der. The resulting spin model can be treated by semiclassical 
spin-wave theory (assuming a magnetically ordered state) or 
by more elaborate methods, to be described below. 

Most of the published calculations were triggered by the 
inela stic neutron scattering experiments of Tranquada et al. 
(fl46h on La 15 / 8 Ba 1 / 8 Cu04, which showed an upper disper- 
sion branch well described by the spectrum of an isolated 
two-leg spin ladder in a situation with period-4 charge or- 
der. Cons e quen tly, coupled two-leg ladders were used in 
Refs. ( fl47t fl48h . Ref. dl47l) employed a linearized bond- 
operator description for a situation inside the magnetically or- 
dered phase close to the QCR i.e., where weak magnetic order 



Magnetic anisotropics will open a gap in the spin-wave spectrum, but are 
often ignored. 
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FIG. 17 Left: 2D square lattice split into two-leg 5 = 1/2 ladders (full dots) coupled by bond-centered charge stripes (open dots). Jn = 
J± = J is the intra-ladder exchange, J cyc is an additional cyclic exchange. Magnetic properties are calculated from the two-leg ladders 
coupled by J'. Right: Excitation spectrum 0Jh,k of coupled two-leg spin ladders as model for period-4 stripes as function of the momenta 
perpendicular (h) and along (k) the ladders. The inter-ladder coupling is ferromagnetic and chosen such that the excitation gap vanishes. The 
spectrum displays low-energy "spin waves" and high-energy "triplons" which meet at a "resonance peak" at (tt, tt). (Reprinted with permission 
from Ref. (148). copyright 2004 by the American Physical Society.) 



co-exists with strong dimerization, and the magnetism is far 
from the semiclassical limit. Ref. (11481) used a more sophisti- 
cated approach of continuous unitary transformations (which, 
however, cannot be easily applied to the ordered phase) right 
at the QCP. 

The two approaches (fl47tfl48l) yield dispersions of effec- 
tively non-interacting bosonic excitations, with rather similar 
results. They nicely show a dual character of the excitations 
(Fig.[T7li: At low energies, the excitations resemble spinwave- 
like Goldstone modes, whereas at higher energies, the charac- 
ter of the triplon dispersion of the single ladder is reproduced. 
The low-energy and high-energy branches meet in a saddle 
point at wavevector (tt, tt). If | J'| is not too large, the energy 
of this saddle point is essentially given by the spin gap of an 
isolated ladder. In fact, the idea of low-energy incommensu- 
rate excitations and the resonance peak being part of a unified 
picture was proposed earlier by Batista et al. (1397b . 

The excitation spectra of the coupled-ladder model com- 
pare favorably to the experimental data of Ref. dl46l) . after 
taking an average of horizontal and vertical stripes (assuming 
stripes running orthogonal to each other in adjacent planes), 
see Fig. [8] in Sec. IIII.BI Once matrix elements are properly 
taken into account, the low-energy signal essentially consists 
of four spots near the ordering wavevectors, forming a square 
with corners along the (1,0) and (0,1) directions. Upon in- 
creasing the energy, the four spots meet at (tt, tt) at what one 
might call a "resonance peak". Above this energy the response 
move away from (tt, tt) and forms a square which is now ro- 
tated by 45 degrees. Within this description, J is taken to be 
130 . . . 160 meV (1 1481) which is the canonical value for 214 
compounds). J' is a fit parameter chosen to place the system 
near criticality; within the ordered phase, the spectrum is only 
weakly sensitive to the value of J'. 

Subsequently, Ref. ffl50h employed spin-wave theory to 
systems of coupled three-leg and four-leg ladders. For both 
situ ation s, reasonable agreement with the experimental data 
of (11461) could be obtained for suitable values of J and J'. 



(For three -leg ladders, the results are qualitatively similar to 
Ref. d398h .) The remarkable fact that semiclassical spin-wave 
theory and strongly quantum mechanical bond-operator the- 
ory produce rather similar results can be related to the fact 
that the anomalous dimension r/ of the magnetic order param- 
eter, characterizing the magnetic QCP, is small (77 = 0.07 for 
the 2d 0(3) model). 

Coupled-ladder models can be employed as well to model 
the magnetism of stripes at smaller doping. Taking the doping 
dependence of the real-space stripe period as input, CDW pe- 
riods of 2M sites are naturally described using coupled even- 
leg ladders with (2A1 — 2) legs; other periods can be mimicked 
by an alternating arrangement of even-leg ladders of different 
width. Such models reproduce essential features of the exper- 
imental spectra, for instance the relation E rrn p R oc x where 
Across is the saddle-point energy at (tt, tt) (13991) . 

Calculations of stripe magnetism beyond localized spin 
models have also been performed. Mean-field plus RPA cal- 
culations of spin and charge-ordere d stripes (i .e. with broken 
symmetries) can be found in Refs. (400c l40lb . The obtained 
spin excitations reasonably describe the low-ene rgy p art of the 
experimental data. However, as seen in Ref. d40ll) . the up- 
per branch does not emerge from the resonance peak and is 
moreover heavily damped, which is not in agreement with ex- 
periment. The problem of not correctly describing the upper 
branch of the hour-glass excitation spectrum is shared by es- 
sentially all RPA calculations, see next subsection. 

A more elaborate time-dependent Gutzwiller approach to 
the Hubbard model was put forward by Seibold and Loren- 
zana (Il49l) . which places the system somewhere in between 
the itinerant and localized regimes. Good agreement with 
the data of Ref. d 1460 was found for bond-centered period-4 
stripes. 

Finally, the collective magnetism of static stripes can also 
be described in a lattice order-parameter theory: Starting from 
a (f) A theory for 2d commensurate magnetism, stripe physics 
can be implemented via uni-directional modulations in mass 
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and gradient terms (f396h . The resulting spin fluctuation spec- 
trum is esse ntiall y iden tical to that of the coupled-ladder two- 
leg model (EH fl48h . with the advantage that the order- 
parameter th eory can be generalized to fluctuating or disor- 
dered stripes dl5ll) . see next subsection. 

What can be learned from these theory excursions? (i) For 
ordered stripes, local-moment approaches work, whereas the 
simplest RPA fails, (ii) While bond-centered stripes work 
marginally better in comparison to experiment, site-centered 
stripes cannot be ruled out (on the basis of the neutron scat- 
tering data). Note that the STM result of ^) on BSCCO and 
CCOC gives very clear evidence for bond-centered charge or- 
der, (iii) Scenarios of checkerboard (instead of stripe) charge 
order can be ruled out (at least for LBCO): A spin-wave calcu- 
lation for checkerboard order (402) gives results in disagree- 
ment with the data, and this conclusion is consistent with the 
calculations in Ref. (115 ID on fluctuating charge order. 



2. Gapped incommensurate spin excitations: RPA vs. 
fluctuating stripes 

The modeling of gapped spin excitations in superconduct- 
ing cuprates goes back to the so-called "resonance mode", 
seen in neutron scatterin g at wavevector (n, t t) and ener- 
gies 30...50meV (EH Ell Ell EH [Hi. RPA and 
related Fermi-liquid-based methods, like slave-bosons plus 
RPA J4-03) or the so-called fluctuation-exchange approxima- 
tion (|404|) . have been very successful in capturing the essen- 
tial features of the early experimental observations: The res- 
onance mode only appears as a sharp mode in the supercon- 
ducting state, where it is pulled below the particle-hole con- 
tinuum (i.e. exists below 2A), whereas only overdamped re- 
sponse exists in the normal state. The energy of the resonance 
roughly scales with the superconducting T c . RPA also de- 
scribes the development of incommensurate excitat ions below 
the resonance energy, albeit with a small inten sity d403l; 1404b . 

T h e neu tron scattering experiments of Ref s. (Il57lll58tll6ll: 
1 1621 : 1 1 63b changed the overall picture. In addition to the 
resonance, both upward and downward dispersing excitation 
branches (the so-called "hourglass") were detected in under- 
doped YBCO. Moreover, sharp features in x" were a ls° de- 
tected above T c . Subsequently, modifications and refinements 
of RPA were proposed, mainly consisting in adjusting the in- 
teraction function g(q) in Eq. (11411 . including its momentum 
dependence. As a result, weak excitations above the reso- 
nance energy could be described d405l) . see Fig. ED RPA was 
also used to describe the aniso tropic low-energy excitations in 
de-twinned YBCO d406t |407|) . see Sec.|VI] However, a com- 
mon feature of all RPA calculations is that a full excitation 
branch above the resonance energy is never obtained, and the 
intensity in this energy range is always weak. 

While Fermi-liquid-based methods like RPA are not ex- 
pected to be applicable for underdoped cuprates, the failure 
of RPA is already apparent in optimally doped YBCO-6.95, 
as discussed in detail in Ref. (170): Neutron scattering ob- 
serves dispersing collective modes above the resonance en- 
ergy, both below and above T c (the latter fact is different from 




FIG. 18 RPA result for y" of a d-wave superconductor (reprinted 
with permission from Ref. (405), copyright 2005 by the American 
Physical Society), as a function of frequency and momentum along 
q = 77(71", 7r) (note that this is not the direction in which the incom- 
mensurate low-energy peaks are commonly observed). The RPA in- 
teraction function was chosen as g(q) — go[l — 0.1(cos q x +cosq y )], 
with go = 0.573 eV. The arrows show the resonance mode around 
(tv, 7r) and the weak excitation above the resonance. 



the early experimental results, possibly due to better resolu- 
tion and crystal quality). In contrast, there is little structure 
in the RPA spin excitations. A similar inapplicability of RPA 
was recently reported for optimally doped BSCCO d 1651) . 

Consequently, strong-coupling theories are called for - 
those should account both for collective spin and charge 
modes. An explicit proposal for a phenomenological order- 
parameter theory of coupled spin and charge fluctuations was 
made in Ref. d396l) . with two crucial ingredients: (i) Both 
spin and charge sectors are defined on the microscopic lat- 
tice; hence the theory accounts for short-wavelength effects 
and lattice pinning, (ii) The bare spin fluctuations live at 
wavevector (it, ir), and any incommensurabilities in the spin 
sector are driven by modulations in the charge sector, via cou- 
plings of the type Re[exp(iQ c ■ x)<fi c (x)]\<f) sa (x)\ 2 : Micro- 
scopically, this reflects that both the spin density and the mag- 
netic couplings get modulated along with the charge. Impor- 
tantly, the magnetic couplings can even change sign, switch- 
ing the spin ordering wavevector from (ir, it) to one dictated 
by the charge order - this renders the effect of the coupling 
non-perturbative. For perfect charge order, the results of this 
approach are essentially identical to that of a microscopic the- 
ory of coupled spin ladders 03961) . 

Ref. (151) proposed an adiabatic approximation for the cou- 
pled order-parameter theory of Ref. d396l) . being appropriate 
for slowly fluctuating or disorder-pinned stripes. By combin- 
ing lattice Monte Carlo simulations for the charge sector with 
exact diagonalization of the spin sector, it was possible to de- 
termine the spin excitations of fluctuating stripes. The full nu- 
merical treatment of the charge sector ensures that positional 
stripe fluctuations as well as fluctuations between horizontal 
and vertical stripes including proper domain walls are cor- 
rectly described. Also, the crossover from stripe to checker- 
board order is part of the theory. Remarkably, the results 
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(0,0) (2it,0) (0,it) q x (2ii,jc) 



FIG. 19 Dynamic spin susceptibility x"(9*i tJ ) f° r bond-centered fluctuating/disordered stripes (reprinted with permission from Ref. il5lh . 
copyright 2006 by the American Physical Society). Left/Middle: cuts at a constant energy, at roughly 60% and 150% of the resonance energy, 
E Tea , as function of momentum. Right: cuts along (q x ,Tr) as function of q x and energy, showing the universal "hour-glass" spectrum, a) 
Strong repulsion between ip x and ip y , correlation length £ ~ 30. b) Weak repulsion, £ « 20. c) Weak attraction, £ ~ 20. d) As in panel b), but 
in the presence of an in-plane anisotropy preferring horizontal stripes. The ratio of the charge gradients in is 1.005. While the anisotropy 
is significant at low energies, uj ~ 60%i?ros (lower panel), it is much less pronounced at 200%-E rC s (upper panel). 



display an hour-glass excitation spectrum over a wide range 
of parameters, see Fig. [19] For strongly repulsive horizontal 
and vertical stripes, Fig.[T9b, the constant-energy cuts display 
distinct high-intensity spots as in the ordered-stripe calcula- 
tion, whereas a situation with weakly repulsive stripes includ- 
ing checkerboard domain walls, Fig. [T9b , results in a rather 
isotropic intensity distribution in q space, despite the under- 
lying stripe physics. The lower branch of the hour glass be- 
comes smeared with decreasing charge correlation length £ c h; 
for £ c h < 10, only a broad vertical feature in u>-q space is 
left at low energies (i.e. a "Y"-shaped dispersion). Including 
a small in-plane anisotropy in the charge sector has a strong 
effect on the magnetic excitations at low energies, whereas 
those at high energies are changed rather little, Fig. [T9ti. This 
energy dependence of the anisotrop y is very s imilar to what is 
experimentally observed in YBCO dl57Ul69l) . Theoretically, 
this behavior can be rationalized considering that the charge 
sector is close to symmetry breaking, and hence is very sen- 
sitive to anisotropies. However, only the lower branch of the 
spin excitations is determined by the properties of the charge 
sector, whereas the upper branch essentially reflects the ex- 
citations of a commensurate gapped antiferromagnet. Note 
that the adiabatic approximation for the CDW modes neglects 
inelastic processes. As discussed in Sec. IV.BI this might be 
particularly appropriate for stripes pinned by quenched dis- 
order, which then are static, but only short-range or dered , as 
also seen in STM ©. While the results of Ref. (fl5ll) are 
encouraging, linking them to microscopic control parameters 
like doping or temperature is difficult. 



A remark on fermionic damping of collective modes is in 
order: Here the RPA and the purely collective-mode descrip- 
tion constitute two extreme cases. In the former, little signal 
is left when the resonance moves into the continuum, whereas 
damping is absent e.g. in the calculation of Ref. dl5ll) (but 
could in principle be included). The experiments give evi- 
dence for collective modes also at elevated energies and above 
T c , albeit with weaker intensity, implying that the truth is in 
the middle. 

In this context, I mention that a recent extension of the 
variational cluster approach (Sec. UV.Cb allows to calculate 
two-particle quantities of the one-band Hubbard model. This 
strong-coupling approach has been sho wn to reproduce an 
hour-glass shaped excitation spectrum d357l) . However, a 
more detailed comparison with experiments has not yet been 
performed. 

Finally, I note that the failure of RPA appears even more 
drastic in electron-doped cuprates. For Pro.ssLaCeo.^CuO.}, 
the intensity of the resonance peak as calculated by RPA is 
by a factor of 10 too small as compared to experiment. Fur- 
thermore, RPA predicts downward dispersing "wings" of the 
resonance which are not seen in experiment (14081) . 



D. Photoemission spectra of stripe phases 

As with neutron scattering, angle-resolved photoemission 
(ARPES) experiments are performed routinely on cuprates. 
Theoretical calculations of ARPES spectra of stripe phases 
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have appeared in a vast number of papers over the last decade. 
Consequently, I will only mention a few important results, and 
their possible connection to experiments. 

On general grounds, translational symmetry breaking in- 
duces Bragg scattering and band backfolding. As a result, 
the Fermi surface breaks up into an arrangement of pockets 
of both e lectro n and hole type as well as some open pieces, 
see Ref. d449l) . However, such a structure is not observable 
in ARPES due to matrix-element effects, as concrete model 
calculations show. 

Mean-field approaches readily provide access to the single- 
particle spectrum. As sketched in Sec. IIV.Bl charge order may 
be implemented by modulated site or bond variables, which 
are either calculated self-consistently within a Hubbard or t- 
J model or imposed ad-hoc onto a single-particle Hamilto- 
nian. Expli cit results for ARPES spect ra were reported e.g. in 
Refs. (ES SHI EH EH EH EH E3). Common prop- 
erties of the results are: (i) ARPES appears little sensitive to 
details of the spatial structure of stripe order, like site vs. bond 
modulations and site vs. bond centering, (ii) Bragg scattering 
opens gaps on quasiparticle branches separated by wavevec- 
tors Q c or Q s . (iii) Shadow features, e.g. Fermi-surface pieces 
shifted by Q c or Q s , are present, but suppressed by matrix- 
element effects, (iv) Stripes induce almost straight Fermi- 
surface segments near the antinodal points located in Q c direc- 
tion, i.e., near (tt, 0) for vertical stripes, (v) For small to inter- 
mediate stripe amplitudes, the Fermi surface crossings along 
the momentum space diagonals are preserved, which also im- 
plies a state with co-existing stripes and <i-wave-like pairing 
possesses zero-energy nodal quasiparticles. For small stripe 
amplitudes, this follows from the ordering wavevector Q not 
connec ting the n odal points of the pure tf-wave superconduc- 
tor J47t 12221 ; |223|) . In contrast, for large stripe amplitudes, the 
results resemble those expected from quasi- Id systems, and 
additional pairing eliminates low-energy quasiparticles. 

In fact , the early photoemission experiments of Zhou et 
al. d242l) on LNSCO-0.12, displaying a cross-shaped low- 
energy intensity distribution and no evidence for nodal quasi- 
particles, were interpreted in terms of Id behavior. The con- 
cept of effectively weakly coupled chains/ladders was fol- 
lowed in the cluster-perturbation theory treatment of the Hub- 
bard model d417l) : Here, ladders were treated by exact diag- 
onalization, whereas the inter-ladder coupling was included 
on the one-particle level in an RPA-like fashion. For half- 
filled stripes with maximal charge modulation, good agree- 
ment with the data of Ref. d242l) were obtained. However, 
as those ARPES results were not confirmed by later experi- 
ment s, see Sec. |III.D~2l the concept of strong Id behavior of 
Ref. ( El7h appears too drastic. 

Other theoretical approaches to stripe ARPES spectra be- 
yond mean field have also been put forward, namely supercell 
DMFT of the Hubbard model (343), a spin-polaron approach 
to the t-J model d419l) . and exact diagonalizatio n of a t—J 
Hamiltonian in the presence of a stripe potential (42c£ l42lh . 
The results are broadly consistent with those from mean-field 
studies, the latter calculations, however, suffered from a poor 
momentum resolution. 

While many theory works only account for perfect charge 



order, static spatial disorder of the stripe pattern can be treated 
at least in simple mean-field theories. The physical picture 
here is that of impurity-pinned static charge order. A simple 
approach i s using a random spacing of uni-directional stripes 
(EH EH EH E3), alternatively short -range-ordered charge 
configurations can be generated from a full Mont e-Ca rlo sim- 
ulation of a stripe order-parameter theory d224t |416|) . From 
the results, it is apparent that stripe signatures in ARPES are 
quickly smeared by a combination of spatial disorder and su- 
perposition of horizontal and vertical stripes. For instance, 
at stripe correlation lengths smaller than 20 lattice spacings, 
clear-cut stripe signatures become essentially invisible d416l) . 
This implies that short-range-ordered stripes, while nicely vis- 
ible in STM, are difficult to detect using probes without spatial 
resolution. The influence of purely dynamic CDW fluctua- 
tions on ARPES lineshapes has also been discussed in a few 
papers. In particular, in has been invoked as explanation for 
broad lineshapes and dispersion kinks 

(EH- 

Interestingly, none of the ARPES experiments on stripe- 
ordered 214 cuprates (after Ref. d242h ) displayed clear-cut 
stripe signatures, perhaps with the exception of recent data 
from LESCO-1/8 d217l) . Here, various features of the data 
have been argued to be consistent with stripes of moderate 
amplitude. However, the behavior near the antinodal points is 
not understood, and the limited experimental resolution ren- 
ders a detailed comparison to theory difficult. 



E. STM spectra 

Measurements of the local density of states (LDOS) using 
scanning tunneling spectroscopy in a number of cuprates have 
triggered enormous theoretical activities. While calculating 
the LDOS of a homogeneous BCS d-wave superconductor is 
a textbook exercise, the most important complications arise 
from intrinsic inhomogeneities and from Mott physics. In the 
following, we shall only discuss theory work on periodic (or 
quasi-periodic) modulations in the STM signal, possibly re- 
lated to symmetry-broken states. In contrast, spatially irreg- 
ular inhomogeneities, reflected e.g. in a broad distribution of 
local gaps in the superconductin g sta te a nd po ssibly related to 
the influence of oxygen dopants d27 u59t 14251) . are outside the 
scope of this review. 

A standard assumption in modeling STM data is that the 
measured signal at low energies locally reflects properties of 
the topmost CUO2 layer. This is not trivial: Commonly used 
crystals of BSCCO cleave such that insulating BiO layers are 
exposed on the surface. This has the advantage of charge neu- 
trality, but implies that electrons may follow a non-trivial tun- 
neling path between the STM tip and the topmost CUO2 layer, 
with the result of a momentum-dependent tunneling matrix 
element (which could also differ between the various cuprate 
families). Theoretical work on this issue has proposed such a 
"filter effect" , with the matrix element being of d-wave type 
for BSCCO d422t |423|) . However, experimental results are 
most consistent with direct, i.e., momentum-independent tun- 
neling - this conclusion can be drawn from the general energy 
dependence of the LDOS, the spatial shape of impurity reso- 
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nances, and the fact that STM results of BSCCO and CCOC 
are grossly similar, despite the different surface layers. Re- 
cent LDA calculations for impurity states also argued 
against a general filter effect, but instead proposed impurity- 
specific features of the wavefunctions to explain details of the 
impurity resonances. 



1. Quasiparticle interference 

Isolated elastic scatterers in an otherwise homogeneous 
metallic system are known to cause Friedel oscillations in 
the local density. An energy-resolved version of this ef- 
fect, dubbed "generalized Friedel oscillations" or "quasipar- 
ticle interference" (QPI), can possibly be observed in the 
LDOS. Then, energy-dependent spatial oscillations in p(f, E) 
will show up in the Fourier- transformed LDOS (FT- LDOS) 
p(k, E) (or its power spectrum). These features typ icall y 
takes the shape of arcs or ridges in momentum space (1 19 lb . 
The existence, intensity, and dispersion of such structures de- 
pend strongly on the shape of the Fermi surface, the presence 
and shape of a superconducting gap, and the nature of the scat- 
terers. Understanding QPI is important in order to be able to 
distinguish it from genuine charge order which also leads to 
real-space oscillations in the LDOS. 

Remar kable , the initial experiments revealing QPI in 
BSCCO dl92l) showed relatively clear peaks in momentum 
space. Those were interpreted using the "joint quasiparticle 
density of states" of a BCS <i-wave superconductor: Peaks in 
p(k, E) were assumed to occur at wavevectors k separating 
those points on iso-energy contours of energy E which have 
a minimal velocity - for a <i-wave superconductor with a typ- 
ical cuprate Fermi surface, these are eight points at the tips of 
the "banana-shaped" iso-energy contours. This idea leads to 
the so-called "octet model": The eight points define a set of 
seven wavevectors qi . . .qj (Fig. [20T>. whose positions appear 
to match the experimental locations of QPI peaks, and agree- 
ment with Fermi sur faces as measured by ARPES has been 
pointed out OH EH). 

In subsequent theory work, the problem of scattering from 
impuritie s in BCS ri-wav e superconductors was studied in de- 
tail (l426t l427t 14281 1429b . For a single point-like impurity at 
site fo, the formal result for the impurity-induced change in 
the LDOS is 



6p(f,E) 



-Im[G(r-fo,E)T(E)G(r -r,E)] (16) 



where G(f, E) is the single-particle Green's function of the 
clean translational-invariant system, and T(E) is the T matrix 
of the scatterer. In the superconducting state, both quantities 
are matrices in Nambu space. This formula clarifies an im- 
portant point: 5p is not only determined by the quasiparticle 
DOS, ImG, but also by real parts of Green's functions. There- 
fore, the octet model is certainly oversimplified. 

Numerical calculations for a single impurity in a <i-wave 
BCS state found structures in the FT-LDOS reminiscent of 
what was seen experimentally. However, the theoretical FT- 
LDOS landscape contained both arcs and peaks, whereas 




FIG. 20 Contours of constant quasiparticle energy for a d- 
wave superconductor with A = 0.6t at energies co/t = 
0.0, 0.075, 0.225, 0.375, 0.5, 0.57. Also shown are the underlying 
normal state Fermi surface, the AF zone boundary, and the seven 
distinct wavevectors q\ of the octet model which connect the banana 
tips. (From Ref. J431» ) 



the experimental data show mainly peaks, and also the in- 
tensi t y dis tribution and precise peak locations did not agree 
(14261 4271) . Interference phenomena be tween mu ltiple impu- 
rities were taken into account in Ref. (EH 429), without a 
significant improvement. It became clear that the momentum- 
space intensity distribution in p(k, E) depends strongly on 
the momentum-space structure of the scatterers. An exten- 
sive theoretical study (14301) of QPI caused by various possi- 
ble scattering sources revealed that smooth scatterers in both 
the particle-hole and the particle-particle channel, possibly 
arising from out-of -plane defects, are crucial in modeling the 
experimental data. Combined with the effect of a small con- 
centration of point-like unitary scatterers, Ref. (|430j) was able 
to reproduce the gross intensity distribution in the experimen- 
tal p(k,E). Nevertheless, a complete understanding of the 
observed FT-LDOS patterns is lacking, although the experi- 
mental existence of dispersing peaks can be associated with 
QPI phenomena with reasonable certainty. 

A few remarks are in order: (i) The octet model describes 
the data remarkably well, with a few caveats: The energy 
dispersion of the experimentally observed qi peaks is signifi- 
cantly weaker, and both <fi and q$ have a higher intensity, than 
predicted by the octet model when using the ARPES disper- 
sions as input, (ii) Why does the octet model work so well? 
This is unclear at present. One ingredient is the strong ve- 
locity anisotropy in the superconducting state, vf S> %, 
which renders the "banana tips" very sharp, but this is not 
sufficient, (iii) Essentially all existing theory work is based on 
non-interacting Bogoliubov quasiparticles of a BCS supercon- 
ductor. The influence of strong correlation effects on QPI has 
not been investigated. Among other things, the quasiparticle 
weights of partially incoherent quasiparticles should enter the 
strength of the QPI signatures. This also implies a strong tem- 
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perature dependence of QPI phenomena. Experimentally, no 
QPI has been detected in the pseudogap regime d200l) . (iv) The 
experimentally observed partial disappearance of QPI signa- 
tures upon approaching the Mott insulator d202l) is not under- 
stood. It is worth pointing out that the absence of QPI peaks 
does not imply the absence of quasiparticles, but may simply 
be related to a change in the structure of iso-energy contours. 
It has been proposed that magnetic ordering tendencies could 
play a role here (l43ll) . but other phases with pocket-like Fermi 
surfaces can lead to a similar effect. 

An exception to (iii) is a recent consideration of QPI near a 
nematic quantum critical point d432l) : Using a one-loop self- 
energy, it was shown that both the velocity anisotropy and 
the dominance of the "banana tips" for the QPI signal are 
enhanced due to scattering off nematic critical fluctuations. 
However, a more detailed study is needed to see whether this 
proposal solves the problems noted above. 



2. Stripes and checkerboards 

Atomically resolved STM spectra from charge-ordered 
states contain a wealth of information on the nature of the 
charge order. The direct real-space picture reveals the spatial 
symmetries (e.g. site centering vs. bond centering), and the 
quantitative information can be used to determine the nature 
of the modulation (e.g. primarily on sites or bonds). 

Within mean-field treatments for both charge order and su- 
perconductivity, the LDOS spectra of stripe states can be eas- 
ily calculated. A comprehensive analysis for weak charge 
order on top of a BCS d-wave superconducor has been pre- 
sented by Podolsky et al. d286l) . In particular, these authors 
calculated the energy dependence of the Fourier component 
of p(k, E) at the charge-ordering wavevector k = Q c which 
was measured in Refs. dl97t Il98l) . Based on their compari- 
son of the spectra of different modulation types, the authors 
concluded that an on-site modulation alone cannot explain the 
STM data, but a sizeable modulation on the bonds, either of 
kinetic energy or pairing strength, is required. Subsequently, 
the LDOS has been calculated for a variety of mean-field the- 
ories of charge-ordered stat es, both f or stripes d319tl415ll433l) 
and for checkerboards (I35u319u333h 

The clear-cut observation of stripe order on the surface 
of both BSCCO and CCOC Q showed that the static or- 
der is only short-ranged, likely due to strong impurity pin- 
ning, is moreover dominated by bond (instead of site) modu- 
lations, and apparen tly co -exists with well-defined low-energy 
quasiparticles ( 190; 202). The STM signatures of short-range 
valence-bond stripe order, together with the interplay wi th im - 
purity scattering, were theoretically investigated in Ref. (|224|) . 
Here, an electronic mean-field theory was combined with 
an order-parameter description of short-range charge order, 
which was used before to model the spin excitations of dis- 
ordered stripes d 1 5 lb . The results show that stripe order is 
strongly visible in the LDOS and Z map at elevated energies, 
in particular near the superconducting gap energy, whereas 
QPI signatures dominate at low energies, Fig. [2T] in good 
agreement with experiment. As pointed out in Sec. III. Al this 
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FIG. 21 Theoretical results for the LDOS p(r, E) at negative (left) 
and positive (middle) bias, together with Z(r, E) (right), for pinned 
short-range ordered stripes with additional impurities on a 64 2 lat- 
tice, at energies 50, 24, 8meV (from top to bottom). While stripe 
signatures dominate at elevated energies, the low-energy data, in par- 
ticular Z( r, E\ display QPI features. (Reprinted with permission 
from Ref. (224), copyright 2008 by the American Physical Society.) 



dichotomy is related to the fact that valence-bond stripes have 
an approximate d-wave form factor and hence display little 
coupling to nodal quasiparticles. 

I now come to field-induced order as seen in STM. The 
initially observed checkerboard patterns around vortices in 
BSCCO d 1 960 triggered a variety of theory works. One 
class of explanations assumed that static collinear SDW or- 
der is induced in or n ear the vortex cores by the applied field 
d434t 14351: l436t 14371) . Then, a static charge-density modu- 
lation is automatically induced. Another class of theories is 
based on pinning of otherwise flu ctuat i ng o rder by a vortex 
core (or another type of impurity) (438; 439). While pinning 
of a CDW trivially induces a charge modulation, the pinning 
effect on a collinear SDW will also lead to a static charge mod- 
ulation, however, static spin order is not required 0439ft . The 
reason is that the pinning potential at position xq couples as 
C0 c (a?o) or C0s(^o) t° CDW and SDW order parameters, re- 
spectively. In the case of pinning of SDW fluctuations, the 
resulting peaks in the FT-LDOS are not energy-independent: 
The collective modes entering e.g. a one-loop self-energy are 
fully dispersive, which causes a dispersion in the FT-LDOS as 
well. This implies a subtle interplay of induced charge order 
and QPI d285h . 

It should be noted that a distinct explanation of the vortex 
checkerboard in terms of fi eld-enhanc ed quasiparticle inter- 
ference has been proposed d2 1 Oc 14401) . without reference to 
collective effects. Indeed, a calculation assuming the exis- 
tence of QPI peaks at the octet model spots i n the FT-LDOS 
shows that both qi and q$ are field-enhanced d440h . This ef- 
fect is related to the coherence factors of the superconduct- 
ing state and may indeed enhance a checkerboard-like signal 
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in the LDOS. While it is at present not unambiguously es- 
tablished which of the interpretations of the vortex checker- 
board is correct, the collected experimental evidence for field- 
induced ordering tendencies in cuprates is overwhelming, see 
Sec. IIII. E. 31 Therefore, it seems very likely that the vor- 
tex STM signal involves an intrinsic ordering phenomenon as 
well, in particular because stripe-like ordering is known d37t) 
to occur at the surface of BSCCO and CCOC, the materials 
where the vortex checkerboard has been studied. 

Finally, let me point out that essentially all model calcula- 
tions for cuprate STM spectra are done for one-band models 
of the CuC>2 plane. Considering that the stripe-like modula- 
tion observed in Ref. (1371) is strongest on the oxygen orbitals, 
this may not be justified: A separate treatment of Cu and O 
orbitals within a three-band model could be required to fully 
understand the physics of bond order in cuprates. 46 



F. Stripes and superconductivity 

Considering that superconductivity and stripes are two 
prominent ordering phenomena in the cuprates, obvious ques- 
tions are: (i) Do stripes and superconductivity co-exist? (ii) 
Do stripes and superconductivity compete or co-operate? (iii) 
Can external tuning parameters employed to tune this inter- 
play? 

Some relevant experiments are described in Sec. Hill above, 
in particular Sec. IIII. El The answer to (i) is yes, most clearly 
seen in LBCO for 0.1 < x < 0.15, where scattering experi- 
ments have established stripe order and thermodynamic mea- 
surements show bulk superconductivity (albeit with a strongly 
suppressed T c ). As for (ii), the suppression of the supercon- 
ducting T c is usually taken as evidence for a competition of 
static stripes and superconductivity. This is consistent with 
the possibility (iii) of magnetic field tuning of incommensu- 
rate spin order, and the fact that (static) stripe signatures have 
not been detected in the multi-layer cuprates with the highest 
T c . 

I note that a few other observations in 214 compounds 
have prompted speculations about a co-operative interplay of 
stripes and superconductivity. For LSCO, the incommensu- 
rate low-energy spin excitations have been found to survive 
into the overdoped regime and to disappear concomitantly 
with superconductivity (Il74l) . However, whether the incom- 
mensurate fluctuations are the cause or the result of pairing 
cannot be deduced from the data. For LBCO- 1/8, the ex- 
perimental data include an unusual gap in the in-plane opti- 
cal conductivity appearing together with charge order (12151) . 
the appare nt d-wave gap in the charge-ordered state above T c 
(1 1361: 1247b. a nd the re sistivity drop at the spin-ordering tem- 
perature d43tl73l : l266h . The latter finding has been interpreted 
in terms of "antiphase superconductivity", to be described in 
more detail in Sec. IV.F.2l below. but alternative explanations 



While bond order can in principle be described in a one-band model, the 
physics of the three-band model is certainly richer. 



have been proposed as well (|444|) . At present, a concise the- 
oretical picture has not emerged, which certainly is related to 
the challenge in understanding the pseudogap regime. 

1. Competing order parameters 

The competition of antiferromagnetism and superconduc- 
tivity has played a central role in the cuprate phenomenology 
early on. It appeared in various t heore t ical flavors, includ- 
ing the SO(5) theory of Zhang dH HqI EH). Also without 
appealing to a higher underlying symmetry, the interplay of 
stripes and superconductivity can be modeled using a coupled 
order-parameter field theory - such an approach assumes the 
presence of the two ordering phenomena without making ref- 
erence to their microscopic origin. 

A concrete theory for the competition of SDW order and su- 
perconductivity in the presence of an ex terna l mag netic field 
has been worked out by Demler et al. (I44U I442I) . with the 
focus on the SDW transition inside the superconducting state. 
The ingredients are a classical Ginzburg-Landau free energy 
for the superconducting condensate if>(x) in the presence of 
an external field H, and a quantum 4 theory for the SDW 
order parameter cf> sa (x,T), with a density-density coupling 
w|V'| 2 |(?! ) sq| 2 - The primary effect of a small applied field is 
to induce vortices in the superconducting order parameter. 
Vortices are accompanied by a suppression of ip in a region 
around the vortex core, such that the balance between SC and 
SDW is locally changed. Importantly, the periodic "potential" 
for the <fr order parameter, resulting from an Abrikosov vor- 
tex lattice, enhances the magnetic fluctuations not only in the 
vortex cores, but over the entire sample - this effect eventu- 
ally causes a field-induced transition from a SC to a SC+SDW 
state. 

The schematic phase diagram from this consideration is in 
Fig-G21 it has been verified by a full numerical analysis of the 
coupled field theory. Among the important results is the be- 
havior of the phase boundary H(s), where s is the tuning pa- 
rameter of the SDW order, near the zero-field quantum phase 
transition at s c : 

H/H a c2 = 2(s - «„)/[« ln(l/(« - s c ))] . (17) 

Remarkably, the phase boundary cannot be obtained from an 
analytic expansion in H, the reason being the infinite diamag- 
netic susceptibility of the superconductor. A related result is 
the behavior of the staggered moment, experimentally mea- 
sured as Bragg peak intensity, as function of the applied field: 

(|0 SQ |) 2 cx s c -s(H), (18) 
s(H) ee S -(vH/(2H° c2 ))HW° c2 /H) 

with s(H) the renormalized cf> s mass, H® 2 the upper critical 
field of the bare superconductor, and z9 a number of order 
unity. Both results (fTTI i and ( fT8l have been found to be quanti- 
tatively obeyedin the neutron scattering experiments of Lake 
et a l. J273t r276h and of Chang et al (274), see Fig. [T4lin 
Sec. IHI.E31 In particular, the field-induced intensity (|0 S al) 2 
increases almost linearly with field near the zero-field SDW 
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FIG. 22 Zero-temperature phase diagram for the field tuning of 
SDW order in the superconducting state (reprinted with permission 
from Ref. J442I) . copyright 2002 by the American Physical Society). 
H is the external field, and s is the SDW tuning parameter, with zero- 
field order for s < s c . The line A-M is the field-tuned SDW tran- 
sition, whereas superconductivity disappears above the line B-M-D. 
Pi and P2 denote paths taken by e xperiment s which start outside 
fvrk HH; HH; Ull; IHl) or inside ( fill Il7lt l27ll ; [2741 ; |275|) 
the SDW phase, respectively. 



transition (12731) . The results of Haug et al. on YBCO-6.45 
(|284|) also appear consistent w ith E q. ([T"8l . 

The theory of Demler et al. j44lH442h . originally restricted 
to T = and two space dimensions, was subse quent ly re- 
fined and extended to include inter-layer coupling d443l) . For 
low order parameter symmetry (Ising or XY), the 3d situation 
admits an interesting quasi- Id ordered phase, where ordering 
occurs within each vortex line individually. Then, the B — ► 
limit of the order-disorder phase boundary does not match the 
B = critical point d443l) . 

A few remarks on the theoretical analysis and its conse- 
quences are in order: (i) The Zeeman effect of the applied 
field has been neglected. Its effect only appears at order H 2 
and is small for H <C J, with J being a typical magnetic ex- 
change energy, (ii) The theory implies that the field-enhanced 
magnetism does not arise primarily from the vortex cores, 
but - via the superflow surrounding each vortex - from the 
entire sample. In fact, if the magnetism would only come 
from the vortex core region, then the observed neutron scat- 
tering intensity would imply large local moments of size 0.6 
lib, as in undoped cuprates, which appears unlikely. How- 
ever, other theories have proposed vortices with antiferromag- 
netic cores d445l) . (Note that the nature of the vortex cores 
in cuprat es is still not understood, but the continuum theory 
of Refs. d44lt |442|) gives a reliable description of the physics 
outside the vortex cores.) (iii) Although the theory focuses 
on the SDW order parameter competing with superconductiv- 
ity, a similar approach could be applied to other field-induced 
orders as well. 

Indeed, an open issue is whether the field directly enhances 
the SDW component of the stripe order or whether the field 
enhances the CDW component. In the CDW is primarily en- 
hanced, one may expect separate transitions associated with 
the CDW and SDW order, with an intermediate SC+CDW 



phase. However, due to pinning of the CDW fluctuations by 
the vortex cores the CDW transition will be smeared. Never- 
theless, a simultaneous study of field-induced SDW and CDW 
intensities in one sample might be able to distinguish between 
the two scenarios. 

Under the assumption that stripes and superconductivity 
compete, the 1/8 anomaly in T c finds a natural explanation in 
terms of the enhanced stripe stability at doping x = 1 /8 (which 
in turn is commonly attributed to lattice commensuration ef- 
fects, but may be enhanced by Fermi surface nesting). I note 
that, in addition, the 2kp CDW instability of one-dimensional 
stripes may contribute to the T c suppression: For nearly half- 
filled stripes, 2kp corresponds to a CDW period of four lattice 
spacings along the stripes, which coincides with the modula- 
tion period perpendicular to the stripes near x = 1/8. Hence, 
coupling between the planes could be efficient in stabilizing 
the 2kp CDW which would then suppress pairing d334l) . 

Last but not least it is worth mentioning that the observa- 
tion of T c w T sp , as observed e.g. in La2CuC»4+,5, is not nec- 
essarily inconsistent with a concept of competing SDW and 
superconducting orders. Within a phenomenological Landau 
theory, which admits a region of phase coexistence of the two 
competing orders at low T, there is a regime of parameters 
with simultaneous onset of superconductivity and magnetism 




2. Fluctuating pairing and antiphase superconductivity 

A remarkable hierarchy of te mper ature scales was found 
in stripe-ordered LBCO-1/8 d73l |266|) . with a sharp resistiv- 
ity drop around 42 K, but no corresponding signature in p c 
and no bulk Meissner effect. This has triggered an interpre- 
tation in terms of 2d fluctu ating superconductivity without 3d 
phase coherence d43l : |266|) . While the resistivity drop seems 
to coincide with the spin ordering temperature in zero field, 
a magnetic field separates the two. The data raise a number 
of questions: (i) What is the nature of the fluctuating pairing 
state and its relation to magnetism? (ii) Why is bulk super- 
conductivity, i.e., a Meissner signal, only established around 
5K? 

A particular scenario, proposed in Refs. J43I : l53l) . is that of 
"antiphase" superconductivity in a stripe state, also dubbed 
"pair density wave". Here, the superconducting condensate is 
modulated, i.e., it is d-wave-like within each stripe, but under- 
goes a 7T phase shift between neighboring stripes. If adjacent 
planes have alternating stripe directions, then it is easy to show 
that the leading-order Josephson coupling between neighbor- 
ing planes vanishes by symmetry, and bulk superconductivity 
will be strongly suppressed. (Higher-order couplings, how- 
ever, will be finite.) 

In the order-parameter language, antiphase superconductiv- 
ity is very similar to an FFLO state, i.e., it is characterized 
by a finite-momentum condensate, (c^ + q j c _£j) 7^ 0, with 

Q p being the pairing modulation wavevector. In the context 
of stripe phases, such a state was first proposed by Himeda 
et al. (42) and then independently by Berg et al. d43h . A 
condensate at wavevector Q p naturally couples to the charge 
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density at wavevector 2Q p , see Sec. III. El hence the pairing 
modulation in a period-4 charge stripe state is of period 8. 
Spectral propertie s of a ntiphase-condensate states have been 
discussed in Ref. d415l) . In contrast to superconductors with 
zero-momentum condensate, the state has a full Fermi surface, 
which is visible mainly as arcs near the nodal direction after 
matrix elements have been taken into account. 

While the idea of vanishing Josephson coupling by sym- 
metry is appealing, a few issues are open: (i) The antiphase- 
condensate state alone does not explain the full hierarchy of 
temperature scales seen experimentally. For instance, in the 
simplest model of weakly coupled planes there is still a sin- 
gle ordering temperature (which approaches the Kosterlitz- 
Thouless temperature in the limit of vanishing inter-plane cou- 
pling). A plausible assumption, also consistent with the global 
cuprate phase diagram, would be that bulk superconductivity 
at low T emerges from a zero-momentum condensate, i.e., 
the transition into the Meissner state involves a different con- 
densate. If this co-exists with the condensate at finite Q p , 
then charge order at wavevector Q p is induced, which is not a 
higher harmonic of the charge order at 2Q p . This is testable 
experimentally. In addition, defects in the stripe order and 
glass-like behavior are likely required to explain the compli- 
cated hierarchy of scales in LBCO d53l) . (ii) T he full Fermi 
surface of the antiphase-condensate state d415l) is not easily 
compatible with the d-wave-like gap observed in ARPES and 
STM experiments on LBCO- 1/8 dl36l: l247h . (iii) It is not 
obvious that the antiphase-condensate state is favored over a 
more conventional superconducting state (f53l) . Whereas vari- 
ational Monte Carlo calculations for the t-J model (42) found 
an interval of t'/t where the antiphase-condensate state is 
sligh t ly favored , results from renormalized mean-field theory 
d328l : 13291 : [330) indicate that the antiphase-condensate state, 
although competitive, is always higher in energy than a stripe 
state with zero-momentum condensate. 

An alternative model for the resistivity drop in LBCO 
invokes spin-liquid physics emerging from nearly straight 
Fermi-surface segments near the antinodal points d444l) . but 
its relation to the established stripe order in LBCO is open. 



regime. The observed quantum oscillations are instead in- 
terpreted as evidence for the presence of a Fermi surface of 
quasiparticles, at least in the regime of the large magnetic 
fields of order 50 T applied in the experiment. 

A follow-up question is concerned with momentum-space 
area of the oscillation orbit which is a direct measure of 
the enclosed Fermi volume in two dimensions. A standard 
analysis of the oscillation frequencies gives a Fermi volume 
somewhat smaller than the nominal hole doping: The oscil- 
lation with frequency (530 T) -1 in YBCO-6.5 translates into 
a Fermi pocket of size 0.075 per Cu, while the doping level 
is about 10%. A similar mismatch is found for YBCO-124. 
(Note that oscillations with a frequency corresponding to the 
electron concentration, i.e., to a large Fermi surface, are not 
observed.) Parallel magnetotransport measurements revealed 
a negative Hall resistance in both materials, which was then 
used to argue in favor of electron (instead of hole) pockets 
(234). A later de-Haas-van-Alphen experiment in YBCO-6.5 
(1448b identified a second oscillation period of (1650 T)~\ as- 



sociated with an additional smaller pocket. 

One line of interpretation is in terms of a conventional 
metallic state where (perhaps field-induced) translational 
symmetry breaking and associated Fermi-surface reconstruc- 
tion (via band backfolding) induces small Fermi surface pock- 
ets. A detailed analysis of such pockets in various candidate 
ordered states shows that (tt, it) antiferromagnetism is insuf- 
ficient, but incommensurate (or long-period commensurate) 
SDW and/or CDW order could in principle lead to hole and 
electron pockets which would cau se o s cilla tions consistent 
with the experimental observations (448; 449). 

While this interpretation may support the notion that static 
stripe-like order is crucial for the phenomenology of under- 
doped high-T c cuprates, care is required. One complication is 
that the experiments are performed in the mixed state. This 
situation was analyzed by Stephen d450h in a regime where 
quasiparticle scattering on vortex lines can be treated pertur- 
batively. Considering that the effective Fermi energy for the 
relevant orbits may be small (i.e. of the same order of mag- 
nitude as the superconducting gap), it is unclear at present 
whether Stephen's analysis is applicable (l45ll) . 



G. Stripes in other experiments 

In this section, we turn to other experimental findings in 
cuprates, which have been interpreted as to provide indirect 
evidence for translational symmetry breaking. 

1 . Quantum oscillations and Fermi-surface reconstruction 

Recently, quantum oscillations have been observed for 
the first time in underdoped cuprates, both in de-Haas-van- 
Alphen and Shubnikov-de-Haas measurements. The com- 
pounds used for these experiments are YBCO-6.5 d446k l448h 
and YBCO-124 dii^), both having exceptionally low disorder 
due to ordered oxygen dopants. The findings came as a sur- 
prise, as some theoretical scenarios assumed the complete ab- 
sence of coherent electronic quasiparticles in the underdoped 



2. Stripes and Raman scattering 

Inelastic light scattering has revealed rich information 
on v ariou s correlated electron systems, for a review see 
Ref. (|452|) . Being a q = probe, it couples to both SDW and 
CDW modes at finite q only via two-particle processes. Nev- 
ertheless, Raman scattering was used to reveal both charge or- 
dering at low temperature and collective CDW motion at high 
temperatures in the ladder system Sri2Cu2404i d453l) . 

For the cuprates, a number of Raman scattering studies have 
been published on the 214 materials (l454t |455|) . In LSCO 
the difference between Bi 9 and B2 9 spectra has bee n use d to 
argue in favor of tendencies toward stripe formation d455l) : At 
low temperature an anomalous low-energy peak shows up in 
the B2 9 channel for x — 0.02 and in the Bi g channel for x = 
0.10 - note that the selection rule of the two geometries are 
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equivalent up to a 45° in-plane rotation. This interpretation is 
supported by a theoretical calculation of the Raman response 
in the presence of soft CDW collective modes d456l) . One 
should, however, note that details of the Raman response are 
not understood d452l) . 

3. Doping evolution of the chemical potential 

Core-level photoemission experiments have been utilized 
to determine the doping dependence of the chemical poten- 
tial p in cuprates. Particularly interesting is the behavior in 
214 materials: While electron-doped Nd2- a; Ce 2 ;Cu04 dis- 
plays a roughly constant chemical-potential slope dp/dx, \x 
is very weakly doping dependent in hole-doped LSCO for 
x < 0.15, while the slope for larger x is similar to the one 
on the electron-doped side d457t 1458b . The anomalous small- 
doping behavior has been attributed to stripe physics: In a pic- 
ture with doping-independent stripe filling pi (which also re- 
sults in a stripe incommensurability e oc x), the chemical po- 
tential p(x) will be constant, whereas a deviation from e oc x 
will induce a slope in p(x). Hence, the behavior of the stripe 
incommensurability as shown in the Yamada plot, Fig. [5] nat- 
urally ties in with the behavior of p(x). Of course, this picture 
is too simplistic, as a constant p(x) implies phase separation, 
and a fixed pi (x) (which could only arise from locking due to 
an incompressibility of the underlying state) implies insulat- 
ing stripes. Nevertheless, a weakly doping-dependent pi for 
x < 1/8 appears consistent with the experimental observa- 
tions. This is also borne out f rom mean-field studies of metal- 
lic stripes, see e.g. Ref. (1133b . 

A few caveats of this interpretation should be noted: (i) No 
static charge stripes have been reported in LSCO (without Nd 
or Eu co-doping). Whether fluctuating stripes would be con- 
sistent with the experimental data is unclear, (ii) The chemical 
potential in CCOC has been found to display a much larger 
slope over the entire doping range as compared to LSCO 
(459). This has been attributed to the periodicity of the charge 
order in CCOC being weakly doping dependent. 

Numerical studies of t— J models found a rather strong in- 
fluence of the longe r-range hopping terms, t' and t", on the 
slope of p(x) d460l) . This might explain the experimentally 
detec ted differences in dp/dx between the cuprate families 
(14591) . Obviously, a more detailed study of the chemical po- 
tential in other cuprates, including those with static charge or- 
der, would be desirable. 



4. Possible signatures of fluctuating stripes 

If the tendency to stripe order is common to many under- 
doped cuprates, but static long-range order is restricted to cer- 
tain 214 compounds, then an obvious question is: "How to 
detect fluctuating stripes?" 47 A rough estimate of a rele- 



Incidentally, this was the title of the 2003 review article by Kivelson et al. 

(HI). 



vant fluctuation frequency, using the STM charge correlation 
length £ c 10 . . . 20 do and a characteristic UV cut-off scale 
of lOOmeV, results in mode frequencies u>f of order of a few 
meV or, equivalently, one THz. 

The observation strategies for fluctuating orders are differ- 
ent for the spin and charge sectors, due to the limitations in ex- 
perimental probes and due to the different effects of quenched 
disorder on both sectors. In the following, we only discuss the 
charge sector; in the spin sector, the low-energy incommen- 
surate fluctuations seen by inelastic neutron scattering are a 
well-established precursor of order. 

A direct observation of stripe fluctuations would be via low- 
energy collective modes in the charge sector, by measuring the 
dynamic charge susceptibility x" c (?> W X or tne dielectric func- 
tion e(q,ui), at wavevectors q ~ Q c and u> ~ uit. Alterna- 
tively, it was proposed that a superconductor with fluctuating 
stripe should display a shear photon mode d46ll) . This mode 
shows up in the dielectric response as well, now at q ~ 1 /£ c 
and oj ~ ujf. In both cases, no present-day experimental tech- 
nique has the required energy resolution, as electron energy- 
loss spectroscopy (EELS) is currently limited to resolutions 
> 0.1 eV. 

Less direct signatures of fluctuating charge order may be 
found in the phonon and spin excitation spectra, in the opti- 
cal and Raman response, as well as in single-particle spectra. 
However, the interpretation of these probes is rarely unam- 
biguous, with examples given throughout this article. 

If an experimental parameter (like pressure or magnetic 
field) is available to continuously tune the system through a 
quantum phase transition associated with stripe order, then the 
system close to quantum criticality can be expected to display 
fluctuating stripes. However, given the large intrinsic energy 
scales of cuprates, such tuning is only possible if the material 
at zero pressure/field is already close to criticality. With the 
exception of field-induced SDW order in LSCO and under- 
doped YBCO (Sec. |III.E.3t , no conclusive experiments in this 
direction have been reported. 

Taking into account impurity pinning, "fluctuating" stripes 
do no longer fluctuate. This simplifies matters, i.e., allows for 
a detection using static probes, but at the same time compli- 
cates matters, because the distinction between ordering ten- 
dencies and impurity effects becomes subtle, see Sees. III.GI 
andHlLC] 



H. Weak-coupling vs. strong-coupling description of 
incommensurate order 

As became clear from the theoretical approaches sketched 
in Sec. |IV| symmetry -breaking order can emerge in concep- 
tually different ways. In metals, weak interactions of low- 
energy quasiparticles can lead to Fermi-surface instabilities. 
In the presence of strong interactions, the possibilities are 
richer and dominated by collective effects instead of Fermi- 
surface properties. 

From a symmetry point of view, weak-coupling and strong- 
coupling limits of density-wave order may be continuously 
connected; exceptions are insulating or non-Fermi liquid or- 
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dered states which cannot be obtained from weak coupling. 
An obvious requirement for a meaningful weak-coupling 
treatment is the existence of well-defined fermionic quasi- 
particles in the relevant temperature range. A feature of 
weak-coupling approaches to cuprate models is that ordering 
is genetically spin-driven, and charge order is parasitic to a 
collinear SDW state. 

A number of quantitative, but possibly significant, dif- 
ferences between weak-coupling and strong-coupling ap- 
proaches can be identified from the available theories (fill) . 
The energy dependence of the dynamic spin susceptibility \" 
provides one criterion: In weak coupling, the intensity will 
be broadly distributed over a wide energy range up to the 
Fermi energy Ep, whereas x" °f the local-moment antifer- 
romagnet is dominated by energies of order of or less than 
the exchange energy J. Furthermore, the strong-coupling 
Xs is typically much more structured in momentum space. 
This is also reflected in the size of the ordered moment in 
an SDW state, which is of order unity and of order T sp /Ep 
(in units of /is) in the strong and weak-coupling limits, re- 
spectively. The physics in the vicinity of the ordering tran- 
sition may provide a second criterion: In weak coupling, the 
order is carried by low-energy quasiparticles, implying spec- 
tral weight transfer over small energy scales only. Moreover, 
the fluctuation regime of a weak-coupling transition is usually 
very narrow. In both respects, the opposite is true for strong 
coupling. A third criterion is provided by the response to 
quenched disorder in a nearly ordered situation: In the weak- 
coupling limit, disorder will primarily scatter quasiparticles 
and thus broaden Fermi-surface-related features, whereas pin- 
ning of low-energy collective modes is dominant in the strong- 
coupling case. Thus, the low-energy parts of both x" ar, d Xc 
will be suppressed (enhanced) by quenched disorder in weak 
(strong) coupling. 

Before coming to experimental data, it should be em- 
phasized that experiments may well be in an intermediate- 
coupling regime, where both approaches describe at least 
part of the data. Moreover, falsifying a weak-coupling RPA 
calculation is much simpler than falsifying the statement 
that "strong-coupling physics is involved": RPA can provide 
numbers, e.g., for spin-fluctuation weights which can be di- 
rectly compared to experiments, whereas strong-coupling ap- 
proaches are diverse and often more phenomenological. 

Having said this, the collected experimental data for stripe 
order in 214 cuprates speak in favor of a strong-coupling de- 
scription: (i) The ordering temperatures obey T c h > T sp , i.e., 
magnetic order is not a prerequisite for charge order, (ii) The 
doping dependence of the CDW wavevector is opposite to 
what is expected from a nesting scenario, (iii) The ordered 
moments in the SDW phases are not small. The energy scale 
of magnetic excitations appears to be set by J. A large frac- 
tion of the spectral weight is found at energies below J and 



There are a few interesting cases where insulating or non-Fermi liquid be- 
havior can be obtained from weak interactions: Perfect Fermi-surface nest- 
ing, Luttinger liquids in strictly one-dimensional sy stem s, or the nematic 
Fermi fluid in a continuum system described in Ref. (385|). 



near wavevector (it, ir). (iv) The fluctuation regime of mag- 
netism is generically wide, (v) Substituting Zn for Cu pins 
stripes in LSCO. (vi) Although cuprate stripes are not insu- 
lating, a well-defined full Fermi surface is not a property of 
stripe compounds: Stripe order emerges from the pseudogap 
regime, not from a well-developed Fermi liquid. I note that 
the existence of low-energy (nodal) quasiparticles near the 
momentum-space diagonals is not necessarily in contradic- 
tion with a strong-coupling perspective: For instance, valence- 
bond stripes display nodal quasiparticles even for sizeable 
modulation amplitudes because of the approximate <i-wave 
form factor of the charge order (1361) . 

The STM observations in BSCCO-2212, BSCCO-2201, 
and CCOC are less clear-cut: Low-energy LDOS modula- 
tions have been reported which display a doping dependence 
of the CDW wavevector compatible with a nesting scenario. 
As these materials are good superconductors (with the excep- 
tion of underdoped CCOC), this may not be surprising: It is 
known that the onset of superconductivity in the cuprates ren- 
ders quasiparticles much more coherent as compared to the 
normal state, opening an avenue for Fermi-surface driven or- 
dering. However, the tunnel-asymmetry maps show a rather 
robust period-4 signal which is particularly strong a elevated 
energies. 

At present, it is unclear whether the data constitute a true 
contradiction. A possibility is that strong-coupling and nest- 
ing effects co-operate in producing ordered phases. 



I. Consequences of nematic order 

Static nematic order of d x i_ y i type may arise, e.g., from a 
Pomeranchuk instability of the Fermi surface or as a precursor 
to stripe order, see Secs. IH.Cl and lIV.Gl Nematic order breaks 
the rotation symmetry in the CUO2 plane from C4 down to 
C2. This has a number of consequences, among them (i) lo- 
cally anisotropic single-particle properties, visible in STM, 
(ii) globally anisotropic single-particle properties, e.g. a dis- 
torted Fermi surface, visible in ARPES, (iii) anisotropic spin 
fluctuation spectra, and (iv) anisotropic transport. With the ex- 
ception of (i), these signatures require an experimental system 
in a single-domain nematic state, as expected in de-twinned 
YBCO. 

About anisotropic STM spectra: It has been proposed 
(fTTl) to employ spatial derivatives of LDOS maps to detect 
nematic order. The simplest quantities are Q xx (r, E) = 
(flg - d 2 y ) P {r,E) and Q xy (r,E) = 2d x d y p(f,E), which 
then should be integrated over some energy interval to re- 
move noise and QPI features to only retain long-wavelength 
information. A practical problem might be that gap modula- 
tions inherent to BSCCO (27) hamper the procedure. To my 
knowledge, a detailed analysis along these lines has not been 
performed. However, the STM data of Kohsaka et al. d37l) 
show stripe-like patterns in the tunnel-asymmetry map which 
obviously break the C4 rotation symmetry locally. 

I now turn to features of single-domain broken rotation 
symmetry. Assuming a Pomeranchuk instability without un- 
derlying charge order, issues (ii) and (iii) have been investi- 
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gated the oretically on the basis of a RVB-type slave-boson 
approach (13791 : [462), with an eye towards the experiments on 
underdoped YBCO. The starting point is a t-J model with 
a small built-in hopping anisotropy of (St/t)o = 5%. The 
mean-field solution then leads to an effective, i.e. correlation- 
enhanced, hopping anisotropy of up to (8t/t) e g ~ 20% at 
low temperature and doping. The renormalized band struc- 
ture is then used to calculate the spin-fluctuation spectrum us- 
ing standard RPA, with a downward-renormalized interaction 
adjusted such that the system at doping 0.07 is close to the 
magnetic QCP. 

In Ref. (13791) . the authors presented a detailed study of the 
RPA spin fluctuation spectrum, as function of temperature 
both ab ove a nd below the pairing temperature (denoted Trvb 
in Ref. (13791) ). The theory reproduces a number of feature s 
found in neutron scattering on de-twinned YBCO 
i.e., a "resonance peak" and anisotropic incommensurate cor- 
relations developing below the pairing temperature. A number 
of differences should also be noted: The anisotropy of the in- 
tensity distribution obtained from RPA is somewhat smaller 
than obtained experimentally in YBCO-6.6, the RPA energy 
dependence of the anisotropy is too weak, and the RPA reso- 
nance peak is somewhat too sharp. The subsequent Ref. d462l) 
studied a non-superconducting state at low doping. The ob- 
tained momentum-space profile of the low-energy spin fluctu- 
ations shows some agreement with the neutron scattering data 
from YBCO-6.45 of Hinkov et al. (19 lb - However, the spin 
correlations from RPA are not incommensurate, in contrast 
to the experimental data, which may be due to the neglect of 
pairing. (Pairing, however, leads to a distinct resonance peak 
in RPA which is not observed experimentally.) 

Th e neu tron scattering data on less underdoped YBCO- 
6.85 d 157b have also been modeled us ing RPA for a fixed 
anisotropic hopping with St/t = 6% d406t |407|) . Although 
the calculations differ in details of the band structure and RPA 
interaction functions, reasonable agreement w ith ex periment 
was found in both cases. In addition, Ref. (|407|) also ac- 
counted for a sub-dominant s-wave pairing component which 
is generically present in YBCO - this was found to lead to a 
90° rotation of the anisotropy pattern as function of energy. 

Possible consequences of a Pomeranchuk scenario have 
also been investigated for 214 compounds d463t 1464b . Slave- 
boson plus RPA calculations reproduce an approximate hour- 
glass spin excitation spectrum. However, the experimentally 
observed doping evolution of the incommensurability is not 
easily recovered in this approach. 

A distinct viewpoint onto nematic phases, e.g. in YBCO, 
is based on underlying fluctuating stripes. Then, spin in- 
commensurabilities are primarily driven by fluctuating charge 
modulations instead of Fermi surface distortions. As shown 
in Sec. IV.C.2I such calculations, albeit based on phenomeno- 
logical input, can describe salient features of the experimental 
data on YBCO as well. 

An important property of all nematic models based on 
Fermi-surface distortions only is that the effective Fermi sur- 
face anisotropies, which are required to fit the neutron data, 
are large, with the Fermi surface topology chan ged compared 
to the undistorted case, see Fig. 4 of Ref. d379b . This should 



be easily detectable in ARPES experiments and hence pro- 
vides a clear-cut distinction to a scenario of a nematic phase 
originating from fluctuating stripes, where the Fermi surface 
distortion for plausible parameters is much smaller (14161) . 



J. Consequences of loop-current order 

Two types of loop-current order have been proposed to oc- 
cur in the pseudogap regime of the cuprates, namely the Cu-O 
loop currents within a unit cell with Q = of Varma J3TI l32h . 
and the d-density wave state with Q = (it, n) of Chakravarty et 
al. d3Qb - In both cases, the resultant orbital antiferromagnetic 
order should be visible in elastic neutron scattering, and the 
relevant experiments were summarized in Sec. IIII.AI Interest- 
ingly, the phase transitions can be relate d to thos e of vertex 
models in classical statistical mechanics (465; 466). We shall 
discuss a few theoretical aspects and additional properties of 
these states in the following. 

The Cu-0 loop-current order of Varma has been originally 
derived from a mean-field theory of the three-band Hubbard 
model. It has a number of interesting properties, which have 
been reviewed in Ref. (1466b . For instance, the fermionic spec- 
trum (at the mean-field level) displays a d-wave-like gap tied 
to the Fermi level, with the quasiparticle energies given by 

Ef=e £ ±D(%) for E^fi, (19) 

where eg is the bare dispersion, D(k) oc cos 2 (2</>)/[l + 
(ej/e c ) 2 ] the mean-field order parameter including form fac- 
tor, <fi is the angle of k, and e c is a band cutoff energy. Thus, 
the T = spectrum consists of four Fermi points along the 
Brillouin zone diagonals, which have been shown to broaden 
into arcs at finite temperatures d467b . The orbital moments 
of the loop-current state should be oriented perpendicular to 
the CUO2 planes. However, the experimentally detected mag- 
netic order at Q = i s cha racterized by moment directions 
canted by roughly 45° (Il37b . Spin-orbit coupling in the low- 
symmetry YBCO structure has been proposed as a source of 
canting d468b . However, the fact that moment directions are 
very similar in HgBa2Cu04+a, with this type of spin-orbit 
coupling being absent due to the tetragonal symmetry, sug- 
gest that currents involve oxygen orbitals outside the CUO2 
planes dl39h . 

The transition into the Varma loop-current state is be- 
lieved 4 9 to be d escribed by a variant of the Ashkin-Teller 
model (466t l47lb . with weak therm odynamic signatures at the 
finite temperature transition (1472b . The quantum critical fluc- 
tuations associated with the breakup of the loop-current order 
have been shown (l466h to be of the scale-invariant form hy- 
pothesized to lead to a Marginal Fermi Liquid d473tl474b . In 
particular, the fluctuations have been proposed to mediate d- 



Not all symmetries of the Ashkin-Teller model are shared by the original 
loop-current model (at the microscopic level). 
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wave pairing d470h as well as to cause a linear- i n-T r esistivity 
(466) and a linear-in-T single-particle lifetime (|469|) . 

Currently, the microscopic conditions for the appearance of 
the Varma loop-current state (e.g. its doping dependence) are 
not well investigated. The reliability of the original mean-field 
approach i s unclear, considering that numerical investigations 
(14751 1476|) of three-band models did not provide evidence for 
loop-current order of sufficient strength to explain the pseudo- 
gap physics (unless stron g hybridizations with apical oxygen 
atoms are included (14761) ). 

Let me now come to d-density wave order (30). This type 
of order appears frequently in RVB-type mean-field theories 
of Hubbard and t—J model (50), although the proposal in 
Ref. (30) was of phenomenological nature and did not rely 
on RVB physics. Within a one-band model, <i-density wave 
order can be characterized by an order parameter of the type 
J2k ^k( c \ + Q c k<y) wave vector Q = (w, tt) and form 
factor dz — cos k x — cos k y . (Also incommens urate variants 
of d-density wave order have been discussed (478) in con- 
nection to the quantum oscillation experiments (448).) For 
Q = (tt, tt), the fermionic spectrum displays again a d-wave 
gap, 

E ± = _fe L+Q ± J V ^ (2()) 

fc 2 2 

with oc dt being the order parameter. Here, the resul- 
tant gap between the two quasiparticle bands is not tied to the 
Fermi level, i.e., in the presence of doping and particle-hole 
asymmetry the state displays Fermi pockets (of both electron 
and hole type) and a dip in the density of states at some fi- 
nite energy away from the Fermi level. This gap behavior has 
been one reason for criticism towards the <i-density wave state 
being an explanation for the pseudogap regime: The gap re- 
ported in STM measurements appea rs to be centered at the 
Fermi level also above T c {mExlllQil). Also > ARPES ex- 
periments did not detect clear-cut evidence for electron pock- 
ets near (tt, 0), expected in a d-density wave state. It has 
been proposed (1477b that long-range correlated disorder in a 
(i-density-wave state may lead to ARPES spectra with Fermi 
arcs, but no pockets near (tt, 0). 

The transition into the <i-density wave state has been argued 
to be the same as that of the 6-vertex model d465h . As above, 
the thermodynamic singularity at the finite-temperature tran- 
sition is weak. To my knowledge, the quantum critical prop- 
erties have not been worked out in detail. 



VI. IMPLICATIONS FOR THE CUPRATE PHASE DIAGRAM 

Understanding the physics of cuprates is still a major chal- 
lenge in condensed matter physics. Although the supercon- 
ducting state appears reasonably well described by BCS-type 
(i-wave pairing, issues of ongoing debate are the nature of the 
pairing mechanism, the non-Fermi liquid normal-state prop- 
erties, and the pseudogap regime at small doping above T c 
(14791) . Ordering phenomena beyond superconductivity and as- 
sociated quantum phase transitions play a central role in this 



debate. 

This final section is therefore devoted to a critical discus- 
sion of broader aspects of lattice symmetry breaking in the 
cuprates. Among other things, it was proposed (i) that stripes 
or stripe fluctuations are a central ingredient to the cuprate 
pairing mechanism, (ii) that stripe quantum criticality is re- 
sponsible for the non-Fermi liquid behavior around optimal 
doping, and (iii) that stripe order is the cause of the pseudo- 
gap in underdoped samples. Similar proposals have also been 
made in the context of other ordering phenomena. 



A. Lattice symmetry breaking: Universality? 

Given the experimental evidence for lattice symmetry 
breaking in a variety of cuprates, a crucial question is that of 
universality: Which of the described features are special to a 
particular family of compounds, and which may be common 
to all high-T c cuprates? A problem is that not all probes are 
available for all cuprate families (e.g. due to surface problems 
or the lack of large single crystals), making a direct compari- 
son difficult. 

Static stripe order with large spatial correlation length is 
only present in the single-layer 214 cuprates. The STM data 
of Ref. ( 37), showing a period-4 valence-bond glass in under- 
doped BSCCO and CCOC, suggest that these materials have 
a tendency towards stripe order as well - this moderate con- 
clusion remains true even if the static order seen in STM is 
a surface effect only. The idea of universal stripe physics 
is indirectly supported by the hour-glass magnetic excitation 
spectrum, observed in 214 compounds as well as in YBCO 
and BSCCO, which is consistent with a concept of fluctuat- 
ing or spatially disordered stripes, see Sec. IV.C.2I However, 
in some cuprates with particularly high T c , e.g. multilayer Bi 
and Hg compounds, no stripe signatures have been detected to 
my knowledge. 

The experimental data may be summarized in the follow- 
ing hypothesis: The tendency toward bond/stripe order ap- 
pears common to doped cuprates, but competes with <i-wave 
superconductivity. Two trends are suggested by the data: (i) 
The LTT phase of 214 compounds is most effective in host- 
ing stripes, (ii) An increasing number of CUO2 layers per unit 
cell shifts the balance between stripes and superconductivity 
towards homogeneous superconductivity. 

Trend (i) is straightforward to explain: The LTT distor- 
tion pattern induces an electronic in-plane anisotropy which 
is favorable for stripe order. Hence, stripes are stabilized by 
electron-phonon coupling. 

About trend (ii): Why would multilayer cuprates b e less 
stripy? The phase diagram of multilayer Hg compounds (l480b 
indicates that, with increasing number of CUO2 layers, super- 
conductivity is stabilized together with commensurate antifer- 
romagnetism. While the presence of different hole doping lev- 
els in the different layers certainly plays a role here, the data 
prompts a speculation: An increasing number of CUO2 lay- 
ers shifts the spin-sector competition between commensurate 
antiferromagnetism and valence-bond order towards antifer- 
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romagnetism, possibly due to an increasing effective inter- 
layer coupling. Assuming that valence-bond order is the driv- 
ing force for stripe formation, then stripes are suppressed to- 
gether with valence bonds. 

Notably, no signatures of stripe or nematic physics have 
been identified in any of the electron-doped compounds. 
All neutron scattering experiments show that magnetic or- 
der or magnetic fluctuations are peaked at the antiferromag- 
netic wavevecto r Qaf, r ather than away from it, see e.g. 
Refs. d408ll48lll48al483l) . Similarly, no evidence for charge 
inhomogeneities has been reported, but here little experimen- 
tal data is available. In line with the above speculation, the 
increased stability of commensurate antiferromagnetism for 
electron doping may be connected to the absence of stripes. 

Signatures of loop-current order with Q = have been de- 
tected by now in YBCO and HgBa 2 Cu04 +< 5, but not else- 
where. Therefore, it is too early for a judgement regarding 
universality here. 

B. Phase transitions and quantum criticality 

All ordered phases discussed in this article admit zero- 
temperature phase transitions J484I) which could be associated 
with interesting quantum critical behavior. In the cuprates, 
quantum criticality has been discussed widely, on both the ex- 
perimental and theoretical side. 

In the normal state, the most puzzling piece of data is the 
linear temperature depe ndence of t he normal-state resistivity 
around optimal doping (485; 486). This and other experi- 
mental results appear well described in t he fr amework of the 
marginal Fermi liquid phenomenology d473b |474|) which is 
based on the assumed behavior of the electronic self-energy 

E(k F ,co) A (win— (21) 

where oj c is a ultraviolet cutoff energy. The linear-in-T re- 
sistivity has been attributed to quantum criticality early on, 
although reliable transport calculations are scarce. Other 
normal-state observations also point towards a quantum criti- 
cal point near optimal doping, associated with the disappear- 
ance of l ocal-mom ent magnetism (12941) . a change in the Fermi 
surface 235; 487), and a distinct crossover behavior in 
the resistivity d488l) . In contrast, signatures of quantum criti- 
cality inside the superconducting phase have not been conclu- 
sively identified. 

The next two subsections contain a summary on what 
is known theoretically about quantum phase transitions into 
states with broken lattice symmetry. The quantum critical 
behavior depends on whether or not a transition takes place 
in the presence of background superconductivity. Hence, we 
shall distinguish critical theories in the d-wave superconduct- 
ing state and in the metallic normal state, with the former (lat- 



ter) being appropriate for temperatures below (above) the su- 
perconducting T c . We shall focus our discussion to criticality 
in two space dimensions. 



1. Nematic transition 

The first transition to be considered is between a disordered 
state, without broken lattice symmetries, and a nematic state, 
with focus on the d-wave nematic. The continuum and lattice 
situations need to be distinguished. The continuum case, with 
the orde r paramete r being a director field, has been studied 
in Refs. d385l ; l489h . In the following, we concentrate on the 
square-lattice case, where the order parameter is a real scalar, 
i.e. of Ising type. 

Importantly, in the presence of a full Fermi surface, there 
is Landau damping of the order parameter, resulting in a dy- 
namic exponent of z = 3. Hence, the normal-state phase tran- 
sition is above its upper-critical dimension. Critical behav- 
ior can be calculated by a perturbative expansion about the 
Gaussian fixed point, as is standard in the Landau-Ginzburg- 
Wilso n (LGW) (or H ertz-Millis) approach to metallic criti- 
cality d490l; |49 lb |492|) . However, the effect of critical fluc- 
tuations on the fermions is strong due to Q = 0. Physically, 
the Fermi surface becomes soft at the transition (384). The 
dynamical Fermi-surface fluctuations at criticality in d = 2 
have been analyzed recently (I493I) : the electronic self-energy 
scales as w 2 / 3 , thus destroying the Fermi liquid at all wavevec- 
tors except for the momentum-space diagonals (i.e. at "cold 
spots")- 51 The quantum-critical transport scattering rate was 
found to be linear in T except at the cold spots, leading to 
a resistivity varying as p(T) oc T 3 / 2 in the clean limit. In 
contr ast, in an impurity-dominated regime, p(T) — po oc T 
d49l . 

Within higher-dimensional bosonization it has been demon- 
strated that the fermionic correlation functions display "local" 
behavior in the non-Fermi liquid quantum critical regime, i.e., 
spatial correlations remain short-ranged d489l) . Note that this 
result has been derived in the continuum limit. While it has 
been suggested that the quantum critical properties of contin- 
uum and lattice cases are similar, this is not entirely correct 
in d = 2: In the continuum situation, a second type of critical 
mode with z = 2 appears due to the combination of rotational 
invariance, d-wave nature of the order parameter, and Landau 
damping. This leads to a complicated interplay of two time 
scales near the transition, resulting in non-trivial logarithmic 
corrections d495h . 

In the d-wave superconducting state, fermionic excitations 
only exist at the nodal points along the momentum space di- 
agonals. In the present Q = case, the order-parameter fluctu- 
ations acquire a relevant coupling to the nodal fermions. The 
critical theory has been focus of recent work and turns out to 



The competition between antiferromagnetism and VBS order is well 
known and studied in undoped Mott insulators on the square lattice 1 3). 



The Landau damping of the nematic order parameter is not qualitatively 
chang ed by the oj 2 ' 3 self energy as compared to the Fermi-liquid case 

EH). 
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be subtle. In Ref. d47k |496|) . an expansion in e = 3 — d was 
employed, with the result that no sta ble c ritical fixed point 
was found (for cases D and E in 14961) ). This suggested 
a fluctuation-induced first-order transition. However, a recent 
approach based on a 1/Nf expansion, where Nf is the num- 
ber of fermion flavors, uncovere d that a continuous transition 
is possible in d — 2 d497t |498|) . Interestingly, the velocity 
anisotropy of the nodal quasiparticles becomes large near crit- 
icality, in contrast to the theories of QPT to d + is or d + id 
states, where the ordered state preserves the C4 lattice symme- 
try and t he qu antum critical theory takes a Lorentz-invariant 
form d47c 14961) , The strong velocity anisotropy may be impor- 
tant for understanding ARPES experiments and quasiparticle 
interference as measured by STM (432). 



2. CDW and SDW transitions 

Next I turn to transitions into a state with broken transla- 
tional symmetry, which could either display a charge mod- 
ulation only, or both CDW and SDW. The transition into a 
uni-directional density-wave state can occur as direct transi- 
tion from a disordered state, or from a state with nematic or- 
der; the latter can then can be understood as a precursor of 
density-wave order. 

The order parameter <f> of the density wave is a complex 
scalar (CDW) or vector (SDW) and carries a lattice momen- 
tum Q 7^ 0. Commensurate lattice pinning reduces the U(l) 
symmetry of the complex phase of to Zjf. Multiple copies 
of <j) are needed for the inequivalent directions of Q, e.g., a 
cuprate CDW or SDW transition in a tetragonal environment 
requires two fields for horizontal and vertical modulations. 

As above, the presence of a full Fermi surface leads to Lan- 
dau damping, here with z = 2, i.e., the damping arises from 
hot spots (lines) on the Fermi surface connected by Q in 2d 
(3d). In 2d, the normal-state phase transitions are at their 
upper-critical dimension. This can be treated by the standard 
LGW approach, but requires a re-summation of perturbation 
theory, with the result of mean-field behavior supplemented 
by logarithmic corrections. 52 A calculation of the electronic 
self-energy shows that the quasiparticle picture breaks down 
in d < 3, but only along the hot lines. A reliable transport 
theory becomes difficult even for the simplest antiferromag- 
netic transition: In c ontra st to d = 3 where a Boltzmann de- 
scription is possible (500), in d = 2 more elaborate methods 
are required. Ref. d50ll) studied transport near a 2d antiferro- 
magnetic critical point, taking into account vertex corrections, 
but neglecting impurity effects. The result for the resistivity 
shows approximately T-linear behavior over an inter media te 
range of temperatures. However, from the 3d results (500) a 
complex interplay of magnetic and impurity scattering can be 



It has been argued that the LGW approach breaks down for the 2d metallic 
antiferromagnet, as the low-energy modes of the Fermi liquid induce an 
infinite number of marginal operators. As a result, a continuous transition 
with non-trivial exponents emerges (499). 



expected also in 2d which has not been studied to date. 

Concrete applications to cuprates of criticality associated 
with an incommensu rate metal l ic CDW have b een worked out 
by the Rome group d39ll: l39l IBol IBol IBol) . In particular, 
singular scattering near CDW-induced hot spots of the Fermi 
surface has been invoked to explain the non-Fermi liquid char- 
acteristics in both the single-particle and transport properties. 
However, a full transport calculation was not presented. From 
an experimental perspective, a strain-controlled CDW QCP 
has been suggested, based on EXAFS measurements of the 
bond-length distributions in various cuprates (509). 

A non-trivial interplay of nematic and density-wave fluctu- 
ations occurs for a transition from a nematic to a CDW (or 
"smectic") state in the continuum (i.e. without underlying lat- 
tice) in d = 2 d507l) . The reason is the Goldstone (director) 
mode of the nematic phase which causes non-Fermi liquid be- 
havior inside the nematic phase d385l) and str ongl y influences 
the critical properties of the CDW transition 

In the d-wave superconducting state, the fate of the criti- 
cal theory depends on whether or not the ordering wavevector 
Q connects two nodal points. This is not the case without 
fine tuning, and such fine tuning seems to be absent experi- 
mentally. Then the critical theory is that of an insulator, with 
dynamical exponent z = l. While non-trivial exponents will 
occur in order-parameter correlations, the fermions are by- 
standers only. However, this is not the full truth: As discussed 
in Ref. (15061) for the case of a SDW transition in a d-wave su- 
perconductor on the square lattice, the coupling between the 
nodal fermions and an Ising nematic field, constructed from 
the two SDW fields, is irrelevant, but with a tiny scaling di- 
mension. Hence, although fermions and order-parameter fluc- 
tuations formally decouple at the QCP, there will be strong 
damping of the nodal quasiparticles, with a nearly T-linear 
scattering rate, from nematic fluctuations. If, instead Q is fine- 
tuned to connect nodal points, then non-trivial critical behav- 
ior the coupled system of order parameter and fermions can 
be expected, in analogy to Refs. d497l:l498l) . 

Finally, I mention an interesting route to exotic phases aris- 
ing f rom the tendency toward SDW order. As discussed in 
Ref. d442l) . the order parameter for a collinear SDW can be 
written as <fi sa — e l0 n a . This description involves a Z2 de- 
generacy associated with a simultaneous change 9 — * 9 + it 
and n a — > — n a , which may be implemented into a field the- 
ory via a Z2 gauge field. Consequently, it is conceivable to 
have a deconfined phase with "excitation fractionalization": 
This is a phase without broken symmetries, but with topo- 
logical order, where the 9 and n a degrees of freedom form 
separate excitations (442). Similar ideas of "stripe fractional- 
ization" were put forward in the context of fluctuating stripes 
in Ref. (15081) . Experimental signatures of such fractionalized 
phases (which are very different from spin-charge separated 
states in RVB-like theories) have not been identified to date. 



3. Cuprate quantum criticality? 

Various cuprate experiments suggest the existence of a 
quantum critical point near optimal doping, with quantum 
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critical signatures above T c , i.e., in the so-called strange-metal 
regime. The theoretically most challenging observation is 
the linear-in-T resistivity which extends up to high temper- 
atures, namely 600 K in YBCO and almost 1000K in LSCO 
d485tl486h . In the following, I try to summarize the theoretical 
status on this issue. 

In general, quantum critical points can be interacting or 
mean-field-like, depending on whether the theory is below 
or above its upper critical dimensions. Interacting critical 
fixed points display single-parameter sca ling in thermody- 
namics and strong hyperscaling properties d484l) . The hypoth- 
esis of single-parameter scaling in the cuprates was explored 
in Ref. d510h . where a number of scaling laws for thermody- 
namic data were proposed. Tests require detailed measure- 
ments of the doping dependence of the chemical potential and 
the electronic compressibility, which are not available to date. 

Proposals which relate cuprate quantum criticality to tran- 
sitions into states with lattice symmetry breaking face a num- 
ber o f objections: (i) Order-parameter theories of Hertz -Millis 
type (490; 491; 14921) have a dynamical exponent z = 2 or 3, 
and are not below their upper critical dimension in d = 2. Con- 
sequently, hyperscaling is violated, and it is unclear whether 
a robust linear p(T) can be expected. This^applies e.g. to 
the charge-density wave criticality of Ref. d504h . (ii) An or- 
der parameter carrying a finite momentum Q (like stripes) will 
primarily affect low-energy fermions which can be connected 
by Q in momentum space. However, ARPES experiments ap- 
pear to be characterized by the absence of quasiparticles at 
all wavevectors in the strange-metal regime, (iii) Quantum 
criticality associated with a conventional order parameter is 
unlikely to yield quantum critical behavior up to 600 K. 

Let me discuss the last objection in somewhat more detail: 
Usually, quantum critical behavior with well-defined power 
laws requires that the order-parameter correlation length is 
large compared to microscopic scales, and that no other en- 
ergy scales intervene. For simple quantum magnets, charac- 
terized by an exchange scale J, it has been sh own that the 
quantum critical regime can extend up to J/2 (15 1 lh . How- 
ever, it appears unlikely that long-ranged correlations exist 
in optimally doped cuprates at temperatures of several hun- 
dred Kelvin. On the one hand, such correlations, if existing 
in the spin or charge channels, should have been detected 
experimentally, but e.g. neutron scattering has found the 
magnetic correlation length at optimal doping to be around 
two lattice constants at low T. On the other hand, numeri- 
cal investigations of relevant microscopic models using, e.g., 
quantum Monte-Carlo techniques, have not found apprecia- 
ble correlations at optimal doping and elevated temperatures. 
These arguments could be invalid if the order parameter is 
difficult to detect, like nematic order (under the influence of 
quenched disorder) or the circulating current patterns pro- 
posed by Varma (I3ltl32l) : here also objections (i) and (ii) may 
not apply. It remains to explain why phonons, existing in the 
temperature range up to 600 K, do not affect the resistivity. 

It is fair to say that, at present, the robustness of the linear- 
in-T resistivity is a puzzle. Even under the hypothesis of an 
interacting critical fixed point, a linear p(T) in the quantum 
critical regime does not automatically follow. In fact, explicit 



candidates for a linear p(T) are scarce. Perhaps with the ex- 
ceptions of early gauge-theory descriptions of Mott physics 
d512l) and a sc enario of a doped disordered spin-liquid Mott 
insulator d513l) . no robust theoretical explanation seems avail- 
able. This leaves open the options of (A) some unknown form 
of strong quantum criticality or (B) the resistivity not being 
truly linear, in which case the search for a quantum critical 
mechanism could be pointless. 

Int erestingly, recent numerical studies of the 2d Hubbard 
SI 

and t-J models (HH), using cluster extensions of 
DMFT, have provided some hints toward criticalit y aro und 
optimal doping. However, the results of Refs. d348t 13491) are 
at best indicative and provide little phenomenological under- 
standing on the physics of the critical point, moreover they 
disagree with respect to the nature of the phase below opti- 
mal doping. Clearly, more investigations in this direction are 
required. 



C. Pseudogap 

The suppression of low-energy fluctuations in underdoped 
cuprates significantly above T c , dubbed "pseudogap", is cen- 
tral to the cuprate phenomenology (5). The pseudogap tem- 
perature, T*, monotonically decreases with doping, in striking 
difference to the superconducting T c . T* has been suggested 
to extrapolate to the scale J of the magnetic exchange in the 
limit of zero doping, and to vanish either around optimal dop- 
ing or at the overdoped end of the superconducting dome. To 
my knowledge, pseudogap signatures have been unambigu- 
ously identified in all hole-doped cuprates, while in electron- 
doped materials the issue is controversial. 

The list of proposed explanations for the pseudogap is long 
and ranges from genuine Mott-gap physics over preformed, 
phase-fluctuating Cooper pairs to ordered states competing 
with superconductivity. While an extensive discussion of 
pseudogap physics is beyond the scope of this review, I will 
briefly summarize a few important aspects. 

Genuine Mott physics is difficult to describe on a phe- 
nomenological level. In single-site DMFT d342l) . a Mott gap 
occurs in the large-C/ insulating phase of the single-band Hub- 
bard model, however, this phase suffers from an artificial spin 
degeneracy. With spatial correlations included, short-range 
singlet for mation ma y be responsible for a partial gap forma- 
tion J345L I346L 13471) . In the framework of Hubbard-model 
field theories, attempts have been made to identify the Mott 
gap with the dynamics of a charge-2e boson which connects 
the low-energy sector to the upper Hubbard band J5141) . The 
verification of theoretical ideas in this direction is open. 

Signatures of preformed pairs above T c (15 151) have been 
identified in a numbe r of expe riments, most not ably Nerns t 



identified in a number ot experiments, most notably JNernst 
effect measureme nts (5161 : 15171) . photoemission d248t EH), 
and STM studies d519 ). Theoretical calculations, relati ng th e 
Nernst signal to phase fluctuating superconductivity d520h . 
give a plausible description of most of the data. This inter- 
pretation is supported by the observation of fluctuating dia- 
magnetism which often varies in proportion to the Nernst co- 
efficient d52ll) . In general, however, care is required, as there 
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are three sources of a sizeable Nernst signal (15171) : quasipar- 
ticles with a small Fermi energy, vortices (i.e. phase fluctu- 
ations), and short-lived Cooper pairs (i.e. amplitude fluctua- 
tions), and it remains to be seen which is most important. 53 In 
any case, the characteristic onset temperature of pairing fluc- 
tuations, as determined from Nernst effect or STM, is signif- 
icantly below the established pseudogap temperature in un- 
derdoped cuprates. This casts doubts on the assumption of 
preformed pairs being the exclusive source of the pseudogap. 

A third class of proposals links the pseudogap to an or- 
dering phenomenon in the particle-hole channel which is as- 
sumed to compete with superconductivity. This would also 
offer a natural explanation for the suppression of T c in the un- 
derdoped regime. Among the concrete proposals for compet- 
ing phases are spin and charge density waves, e.g., stripes, and 
various forms of circulating-current orders. An common ob- 
jection against these proposals is that the pseudogap line at T* 
does not appear to be associated with a thermodynamic phase 
transition. However, this can be circumvented either by invok- 
ing quenched disorder which tends to smear the transition (see 
Sec. III. Gi l or by postulating a special form of phase transition 
with weak thermodynamic singularities, e.g. of Kosterlitz- 
Thouless or Ashkin-Teller type d472l) . In fact, a weak, but 
rather sharp, signature in the uniform susceptibility has re- 
cently been detected in YBCO samp les o f different doping, 
tracking the pseudogap temperature (Il45l) . This may repre- 
sent a distinct thermodynamic phase transition, and is possi- 
bly connected to the broken time reversal as detected in polar 
Kerr effect measurements dl4ll) . 

In the context of the present article's topic, I will briefly 
discuss the hypothesis that stripes are the cause of the pseudo- 
gap, which has bee n voic ed on the basis of both phenomeno- 
logical theory d504t 15051) and experimental data d207l) . In my 
view, this hypothesis is problematic on several grounds: (i) 
Order in the particle-hole channel with a finite wavevector Q 
will cause distinct signatures in the quasiparticle band struc- 
ture, in particular gaps at momenta separated by Q. Those 
have not been observed, (ii) The pseudogap appears to be a 
universal phenomenon in hole-doped cuprates, with very sim- 
ilar properties in the different families. This suggests that the 
pseudogap has a common origin in all hole-doped cuprates. 
In contrast, stripes are strongest in single-layer compounds of 
the 214 family, but are weak or absent e.g. in materials with 
more than two CUO2 layers per unit cell. 

On the theory side, microscopic calculations using cluster 
extensions of DMFT have established the existence of a pseu- 
dogap in the 2d Hubbard model at small d oping above the 
superconducting instability d345t l346t 13471) . Unfortunately, 
the numerical data provide limited insight into the origin of 
gap formation: The pseudogap occurs in the absence of long- 
range order and is apparently related to strong short-range cor- 
relations (which may eventually become long-ranged at low 



Nernst effect measurements in LESCO and LNSCO have observed a dis- 
tinct enhancement of the Nernst signal far above T c , which was interpreted 
in terms of a stripe-induced reconstruction of the Fermi surface 4235h - 



temperatures). Whether this pseudogap should be attributed 
to genuine Mott physics is unclea r; alt ernatively, both antifer- 
romagnetic d345l) and bond-order d347l) correlations have been 
made responsible for the pseudogap. 

In summary, experimental and theoretical results consis- 
tently show that forms of non-superconducting order are en- 
hanced at small doping, i.e., in the pseudogap regime. This in- 
cludes spin and charge density waves, nematic order, and pos- 
sibly loop-current order. Then, such order can either be (A) 
the cause of the pseudogap, in which case the phenomenon 
should be common to all cuprates, or (B) a secondary effect of 
some other phenomenon causing the pseudogap. For stripes, 
experimental evidence points towards scenario (B). For other 
forms of order, more experiments are required to check or ver- 
ify their universal occurrence. 



D. Pairing mechanism 

While it appears well established that the superconduct- 
ing state in the cuprates is characterized by Cooper pairs 
of d-wave symmetry (38), the pairing mechanism is contro- 
versial. Weak-coupling calcu l ation s in the Hubbard-model 
framework, using RPA (522; 15231) or the more sophisti- 
cated functional renormalization group (fRG, see Sec. |IV.DT2l ) 
d358l) . show that rf-wave pairing can in principle be mediated 
by antiferromagnetic fluctuations. Electron-phonon coupling, 
various experimental signatures for which have been identi- 
fied, has also been suggested as pairing force - this, however, 
is not easily compatible with d-wave symmetry. The distinct 
maximum in T c as function of doping has triggered alterna- 
tive proposals, with strong pairing driven by the presence of a 
quantum critical point near optimal doping. 

Stripe physics as as a source of pairing was suggested by 
different g roups, either via quantum critical CDW fluct uations 
d502l ; l503l) . or via topological stripe fluctuations (508), or via 
the interp lay o f co-existing Luttinger-liquid and hole-pair ex- 
citations d524l) . At present, it is fair to say that none of these 
suggestions has been worked out into a testable theory of d- 
wave superconductivity. 

A diff erent but relate d proposal is that of Kivelson and co- 
workers d525tl526c 15271) who argue that the presence of stripes 
in a superconducting system enhances the transition temper- 
ature T c . This idea starts from the assumption of a pairing 
scenario where the superconducting gap and stiffness are anti- 
correlated as a function of doping, as is the case in cuprates 
below optimal doping. Then, one can show that an inhomo- 
geneous structure, combining regions of large gap with those 
of large stiffness, leads to an increased T c as compared to the 
homogeneous system, at the expense of having a smaller su- 
perfluid stiffness. This concept of "optimal inhomogeneity" 
was proposed to be relevant for understanding the high transi- 
tion temperature in cuprates. 

While at present it is difficult to judge these proposals, they 
meet the same problem as was noted above: The empirical 
anticorrelation between "stripyness" and high T c points to- 
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wards stripes being a competitor to superconductivity. The 
possibility that stripe fluctuations at elevated energies mediate 
pairing cannot be ruled out, however. 



VII. CONCLUSIONS 

Ordered phases in strongly correlated electronic systems 
offer a fascinating variety of phenomena. In this article, I 
aimed to give a balanced and critical account on a particu- 
larly intensively studied class of ordered phases, namely those 
associated with lattice symmetry breaking in cuprate high- 
temperature superconductors. I tried to cover exciting experi- 
mental developments, which include the observation of charge 
order via resonant soft x-ray scattering and scanning tunnel- 
ing microscopy and the identification of a seemingly universal 
spin excitation spectrum. I also discussed theoretical works, 
in particular those dealing with the microscopic origins of 
symmetry breaking as well as those providing a more phe- 
nomenological modeling of experimental data. Considering 
the wealth of published papers and the existence of previous 
review articles, the main emphasis was on recent works which 
appeared during the last five years. 

A quick physics summary may be given as follows: The 
tendency towards states with modulated spin and charge den- 
sities (stripes) appears to be common to hole-doped cuprates 
- a view which has only been established recently - although 
the strength of the phenomenon varies from family to family. 
While general scenarios of stripe formation have been devel- 
oped, the underlying microscopic physics is not completely 
understood. In particular, oxygen orbitals appear to play an 
important role, implying that the three-band Hubbard model 
has to be considered to gain a quantitative picture. Although 
tendencies to stripes are strongest in the pseudogap regime, 
stripes may well be a result rather than a cause of the pseu- 
dogap. The same applies to other experimentally identified 
symmetry-broken states, e.g., nematic and loop-current order. 
Indications for these ordering phenomena being crucial ingre- 
dients to the cuprate pairing mechanism are weak at present. 

Independent of their actual role for cuprate superconduc- 
tivity, stripes and nematics appear to be common to a vari- 
ety of correlated oxides, ranging from cuprates to nickelates 
J22I : I23I) . manganites 1241) . and perhaps also ruthenates J46h . 
While this suggests a common driving mechanism, differ- 
ences are apparent: While ruthenates are good metals, stripe 
phases in nickelates and manganites are insulating and display 
rather robust charge order, rendering them more "classical". 
Thus, cuprates are indeed special, as they offer a unique com- 
bination of the proximity to a Mott insulator and very strong 



The resistivity drop in LBCO-1/8, occurring at T Bp and interpreted as 
fluctuating superconductivity (73; 266), together with the fact that pairing 
sets in at a lower T for all other doping levels of LBCO, has been invoked 
as evidence for a positive correlation between stripes and pairing. However, 
a concise picture for the data of Refs. 1 73; 266) has not yet emerged, see 
Sec. lVR2l 



quantum effects. Understanding their puzzles remains a chal- 
lenge. 
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